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Modified gravity theories have received increased attention lately due to combined motivation 
coming from high-energy physics, cosmology and astrophysics. Among numerous alternatives 
to Einstein's theory of gravity, theories which include higher order curvature invariants, and 
specifically the particular class of /(-R) theories, have a long history. In the last five years there has 
been a new stimulus for their study, leading to a number of interesting results. We review here /(-R) 
theories of gravity in an attempt to comprehensively present their most important aspects and 
cover the largest possible portion of the relevant literature. All known formalisms are presented — 
metric. Palatini and mctric-afSne — and the following topics are discussed: motivation; actions, 
field equations and theoretical aspects; equivalence with other theories; cosmological aspects and 
constraints; viability criteria; astrophysical applications. 
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I. INTRODUCTION 
A. Historical 

As we are approaching the closing of a century after the 
introduction of General Relativity (GR) in 1915, ques- 
tions related to its limitations are becoming more and 
more pertinent. However, before coming to the contem- 
porary reasons for challenging a theory as successful as 
Einstein's theory, it is worth mentioning that it took only 
four years from its introduction for people to start ques- 
tioning its unique status among gravitation theories. In- 
deed, it was just 1919 when Weyl, and 1923 when Edding- 
ton (the very man that three years earlier had provided 
the first experimental verification of GR by measuring 
light bending during a solar eclipse) started considering 
modifications of the t heory by including higher orde r in- 
variants in its action (|Eddingtonl . Il923tlwevll . Il919l ). 
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These early attempts were triggered mainly by scien- 
tific curiosity and a will to question and, therefore, un- 
derstand the then newly proposed theory. It is quite 
straightforward to realize that complicating the action 
and, consequently, the field equations with no apparent 
theoretical or experimental motivation is not very ap- 
pealing. However, the motivation was soon to come. 

Beginning in the 1960's, there appeared indications 
that complicating the gravitational action might indeed 
have its merits. GR is not renormalizable and, therefore, 
can not be conventionally quantized. In 1962, Utiyama 
and De Witt showed that renormalization at one-loop de- 
mands that the Einstein-Hilbert action be supplemented 
by hi gher order curvature terms ( Utivama and DeWittl . 
[1962!). Later on, Stelle showed that higher or der ac- 
tions are indeed renormalizable (but not unitary) (jStelld . 
Il977f ). More recent results show that when quantum 
corrections or string theory are taken into account, 
the effective low energy gravitatio nal action admits 
higher order curvature inva r iants (iBirrell and Daviesl . 
119821 : iBuchbinder et aLl . [19921 : IVilkoviskvlll992[ ). 

Such considerations stimulated the interest of the 
scientific community in higher-order theories of gravity, 
i.e., modifications of the Einstein-Hilbert action in order 
to include higher-order c urvature i nvaria nts with respect 
to the Ricci scalar [see (ISchmidtl . [200I for a historical 



review and a list of references to early work]. However, 
the relevance of such terms in the action was considered 
to be restricted to very strong gravity regimes and 
they were expected to be strongly suppressed by small 
couplings, as one would expect when simple effective 
field theory considerations are taken into account. 
Therefore, corrections to GR were considered to be 
important only at scales close to the Planck scale and, 
consequently, in the early universe or near black hole sin- 
gularities — and indeed there are relevant studies, such 
as the well-k i iown curvature-driven inflation scenario 
( Starobinskvl Il980l ) and attempts to avoid cosmo- 



1992. 1993. 1995: 


Brandenbereer et al, 1993 


Mukhanov and Brande 


nbereeil 119921 Shahid-Saless 


19901: Trodden et al. 


I993I). However, it was not 



expected that such corrections could affect the gravita- 
tional phenomenology at low energies, and consequently 
large scales such as, for instance, the late universe. 



matter refers to an unkown form of matter, which has 
the clustering properties of ordinary matter but has not 
yet been detected in the laboratory. The term dark en- 
ergy is reserved for an unknown form of energy which 
not only has not been detected directly, but also does 
not cluster as ordinary matter does. Mor e rigorously , 
one could use the various energy conditions (IWaldl . ll984D 
to distinguish dark matter and dark energy: Ordinary 
matter and dark matter satisfy the Strong Energy Gondi- 
tion, whereas Dark Energy does not. Additionally, dark 
energy seems to resemble in high detail a cosmological 
constant. Due to its dominance over matter (ordinary 
and dark) at present times, the expansion of the universe 
seems to be an accelerated one, contrary to past expec- 
tations.^ 

Note that this late time speed-up comes to be added to 
an early time accel erated epoch as predicted by th e infla- 
tionary paradigm (|Guthl . Il98ll: iKolb and Turneil . Il992t 
iLinde . 1990) . The inflationary epoch is needed to ad- 
dress the so-called horizon, flatness and monopole prob- 
lems (KoVo and Turneil Il992l: [Lindd . Il990l: [Misned . Il968t 
[Weinberg, 197^, as well as to provide the mechanism 
that generates primordial inhomogeneities ac ting as seeds 
for th e formation of large scale structures ( Mukhaiiovl . 
I2OO3D . Recall also that, in between these two periods 
of acceleration, there should be a period of decelerated 
expansion, so that the more conventional cosmological 
eras of radiation domination and matter domination can 
take place. Indeed, there are stringent observational 
bounds on the abundances of light elements, such as deu- 
terium, helium and lithium, which require that Big Bang 
Nucleosynthesis (BBN), the production of nuclei other 
than hydrogen, takes place during radiation domin ation 
(|Burles et all I2OOII : iGaTroll and Kaplinghatl l2002t ). On 
the other hand, a matter-dominated era is required for 
structure formation to take place. 

Puzzling observations do not stop here. Dark mat- 
ter does not only make its appearance in cosmological 
data but also in astrophysical observations. The "missing 
mass" q uestion had alr eady been posed in 1933 for galaxy 
clus ters (IZwickv . 1933 ) and i n 1959 for individual galax- 
ies (|Kahn and Woltieil . ll959l) and a satisfactory final an- 



swer has bee n pending ever sin ce (Bosma, 1978; Ellii 



20021: iMoord. [2001 : Persic et all ,1996.: Rubin and Ford 
19701 : iRubin efali [l98oir 



One, therefore, has to admit that our current picture 



B. Contemporary Motivation 



More recently, new evidence coming from astrophysics 
and cosmology has revealed a quite unexpected picture of 
the universe. Our latest datasets coming from different 
sources, such as the Gosmic Microwave Background Radi- 
ation (GMBR) and supernovae surveys, seem to indicate 
that the energy budget of the universe is the following: 
4% ordinary baryonic mat t er, 20% dark matter and 76% 



dark erterm (lAstier et al . 2006: Ei senstein et l2005t 
iRiess et all [20041: ISpergel aU 2^)- The term dark 



^ Recall that, from GR in the absence of the cosmological con- 
stant and under the standard cosmological assumptions (spatial 
homogeneity and isotropy etc.), one obtains the second Fried- 
mann equation 

a A-kG , 

- = — (P-H3P), (1) 

a i 

where a is the scale factor, G is the gravitational constant and p 
and P are the energy density and the pressure of the cosmological 
fluid, respectively. Therefore, if the Strong Energy Condition 
p + 3P > is satisfied, there can be no acceleration (gravity is 
attractive) . 
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of the evolution and the matter/energy content of the 
universe is at least surprising and definitely calls for an 
explanation. The simplest model which adequately fits 
the data creating this picture is the so called concordance 
model or ACDM (A-Cold Dark Matter), supplemented by 
some inflationary scenario, usually based on some scalar 
field called inflaton. Besides not explaining the origin 
of the inflaton or the nature of dark matter by itself, the 
ACDM model is burden ed with t he wel l known cosmolo g 
ical constant problems ( Carrolll boOlat IWeinberd . Il98£ 



the magnitude problem, according to which the observed 
value of the cosmological constant is extravagantly small 
to be attributed to the vacuum energy of matter fields, 
and the coincidence problem, which can be summed up 
in the question: since there is just an extremely short 
period of time in the evolution of the universe in which 
the energy density of the cosmological constant is com- 
parable with that of matter, why is this happening today 
that we are present to observe it? 

These problems make the ACDM model more of an 
empirical fit to the data whose theoretical motivation 
can be regarded as quite poor. Consequently, there 
have been several attempts to either directly motivate 
the presence of a cosmological constant or to propose 
dynamical alternatives to dark energy. Unfortunately, 
none of these attempts are problem-free. For instance, 
the so-called anthropi c reasoning for the m agnitude of 
A ()Barrow and Tipleil . Il986t ICarteil . Il974[ ). even when 
placed into the firmer grounds through the idea of the 
"anthropic or string landscape" (,Susskind . .2003) , still 
makes many physicists feel uncomfortable due to its 
probabilistic nature. On the other hand, simple sce- 
narios for dynamica l dar k energy, such as quintessence 



(iBahcall et al. 



1998; 



. 199 9: Cald well et a/.l, 119981: iCarrof 

Ostriker and Steinhardt, 1995^ 'Peebles and Rat 
Ratra and Peebles, 1988; Wang et al, 2000; We tterich 



1988; 



1988t ) do not seem to be as well motivated theoretically 
as one would desire.^ 

Another perspective towards resolving the issues de- 
scribed above, which might appear as more radical to 
some, is the following: gravity is by far the dominant in- 
teraction at cosmological scales and, therefore, it is the 
force governing the evolution of the universe. Could it be 
that our description of the gravitational interaction at the 
relevant scales is not sufficiently adequate and stands at 
the root of all or some of these problems? Should we con- 
sider modifying our theory of gravitation and if so, would 
this help in avoiding dark components and answering the 
cosmological and astrophysical riddles? 

It is rather pointless to argue whether such a perspec- 



^ We are referring here to the fact that, not only the mass of the 
scalar turns out to be many orders of magnitude smaller than any 
of the masses of the scalar fields usually encountered in particle 
physics, but also to the inability to motivate the absence of any 
coupling of the scalar field to matter (ther e is no mechanism or 
symmetry preventing this) l lCarro 



five would be better or worse than any of the other so- 
lutions already proposed. It is definitely a different way 
to address the same problems and, as long as these prob- 
lems do not find a plausible, well accepted and simple, 
solution, it is worth pursuing all alternatives. Addition- 
ally, questioning the gravitational theory itself definitely 
has its merits: it helps us to obtain a deeper understand- 
ing of the relevant issues and of the gravitational interac- 
tion, it has high chances to lead to new physics and it has 
worked in the past. Recall that the precession of Mer- 
cury's orbit was at first attributed to some unobserved 
("dark") planet orbiting inside Mercury's orbit, but it 
actually took the passage from Newtonian gravity to GR 
to be explained. 



C. f{R) theories as toy theories 

Even if one decides that modifying gravity is the way to 
go, this is not an easy task. To begin with, there are nu- 
merous ways to deviate from GR. Setting aside the early 
attempts to generalize Einstein's t heory , most of which 
have been shown to be non-viable (IWil]| . ll98lD . and the 
most well known alternative to GR , scalar-te nsor t heory 
(iBergmann . 1968 ; Brans and Dick c, 1961; D ickd , Il962l : 
iFaraonil l2004aHNordtvedtl , Il970l: IWagoneii I1970[) . there 
are still numerous proposals for modified gravity in con- 
temporary literature. Typical examples a r e DGP (Dvali- 
Gabadadze-Porrati) gravity to vali et al. I, l2000D . brane- 
world gr avity ( Maartens, ^QM) , TeVeS (Tensor- Vector- 
Scalar) (iBekens tein, 2004 ) and Einstein- Aether theory 
( Jacobson and Mattinglvr[200l[ ). The subject of this re- 
view is a different class of theories, f{R) theories of grav- 
ity. These theories come about by a straightforward gen- 
eralization of the Lagrangian in the Einstein-Hilbert ac- 
tion. 



'S'eh — — 



(2) 



where K = 87rG', G is the gravitational constant, g is 
the determinant of the metric and R is the Ricci scalar 
(c = h = 1), to become a general function of R, i.e., 



(3) 



Before going further into the discussion of the details 
and the history of such actions — this will happen in 
the forthcoming section — some remarks are in order. 
We have already mentioned the motivation coming from 
high-energy physics for adding higher order invariants 
to the gravitational action, as well as a general motiva- 
tion coming from cosmology and astrophysics for seeking 
generalizations of GR. There are, however, still two ques- 
tions that might be troubling the reader. The first one 
is: Why specifically f{R) actions and not more general 
ones, which include other higher order invariants, such 
as R^,^R^"'1 
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The answer to this question is twofold. First of all, 
there is simplicity: f{R) actions are sufficiently gen- 
eral to encapsulate some of the basic characteristics of 
higher-order gravity, but at the same time they are sim- 
ple enough to be easy to handle. For instance, viewing / 
as a series expansion, i.e., 



i?3 



, (4) 



where the ai and f3i coefficients have the appropriate 
dimensions, we see that the action includes a number 
of phenomenologically interesting terms. In brief, f{R) 
theories make excellent candidates for toy-theories — tools 
from which one gains some insight in such gravity mod- 
ifications. Second, there are serious reasons to believe 
that f{R) theories are unique among higher-order grav- 
ity theories, in the sense that they seem to be the only 
ones which can av oid the long kno wn and fatal Ostro- 
gradski instability (jWoodardl . l2007l) . 

The second question calling for an answer is related to 
a possible loophole that one may have already spotted in 
the motivation presented: How can high-energy modifi- 
cations of the gravitational action have anything to do 
with late-time cosmological phenomenology? Wouldn't 
effective field theory considerations require that the co- 
efficients in eq. (|H) be such, as to make any corrections to 
the standard Einstein-Hilbert term important only near 
the Planck scale? 

Conservatively thinking, the answer would be posi- 
tive. However, one also has to stress two other serious 
factors: first, there is a large ambiguity on how grav- 
ity really works at small scales or high energies. Indeed 
there are certain results already in the literature claiming 
that terms responsible for late time gravitational phe- 
nomenology might be predicted by some more funda- 
mental theory, such as strin g theory [see, for instance, 
(jNoiiri and Odintsovl [2003b[ )]. On the other hand, one 
should not forget that the observationally measured value 
of the cosmological constant corresponds to some energy 
scale. Effective field theory or any other high-energy the- 
ory consideration has thus far failed to predict or explain 
it. Yet, it stands as an experimental fact and putting 
the number in the right context can be crucial in ex- 
plaining its value. Therefore, in any phenomenological 
approach, its seems inevitable that some parameter will 
appear to be unnaturally small at first (the mass of a 
scalar, a coefficient of some expansions, etc. according to 
the approach). The real question is whether this initial 
"unnaturalness" still has room to be explained. 

In other words, in all sincerity, the motivation for in- 
frared modifications of gravity in general and f{R) grav- 
ity in particular is, to some extent, hand-waving. How- 
ever, the importance of the issues leading to this motiva- 
tion and our inability to find other, more straightforward 
and maybe better motivated, successful ways to address 
them combined with the significant room for speculation 
which our quantum gravity candidates leave, have trig- 
gered an increase of interest in modified gravity that is 



probably reasonable. 

To conclude, when all of the above is taken into ac- 
count, f{R) gravity should neither be over- nor under- 
estimated. It is an interesting and relatively simple alter- 
native to GR, from the study of which some useful con- 
clusions have been derived already. However, it is still 
a toy-theory, as already mentioned; an easy-to-handle 
deviation from Einstein's theory mostly to be used in or- 
der to understand the principles and limitations of mod- 
ified gravity. Similar considerations apply to modifying 
gravity in general: we are probably far from concluding 
whether it is the answer to our problems at the moment. 
However, in some sense, such an approach is bound to be 
fruitful since, even if it only leads to the conclusion that 
GR is the only correct theory of gravitation, it will still 
have helped us to both understand GR better and secure 
our faith in it. 



II. ACTIONS AND FIELD EQUATIONS 

As c an be found in many textbooks — see, for ex- 
ample (jMisner et all \197A IWaldl . I1984D — there are 
actually two variational principles that one can apply 
to the Einstein-Hilbert action in order to derive Ein- 
stein's equations: the standard metric variation and a 
less standard variation dubbed Palatini variation [even 
th ough it was Einstein and not Palatini who introduced 
it ( Ferraris et aLl . [T982[ )]. In the latter the metric and the 
connection are assumed to be independent variables and 
one varies the action with respect to both of them (we 
will see how this variation leads to Einstein's equations 
shortly), under the important assumption that the mat- 
ter action does not depend on the connection. The choice 
of the variational principle is usually referred to as a for- 
malism, so one can use the terms metric (or second order) 
formalism and Palatini (or first order) formalism. How- 
ever, even though both variational principles lead to the 
same field equation for an action whose Lagrangian is lin- 
ear in R, this is no longer true for a more general action. 
Therefore, it is intuitive that there will be two version of 
f{R) gravity, according to which variational principle or 
formalism is used. Indeed this is the case: f{R) gravity 
in the metric formalism is called metric f{R) gravity and 
f{R) gravity in the Palatini fo rmalism is called Palatini 
fiR) gravity (iBuchdahll . [l970l) . 

Finally, there is actually even a third ver- 
sion of f(R) gravit y : m etric-affine f{R) gravity 
(|Sotiriou and Liberatil . l2007allbl ). This comes about 



if one uses the Palatini variation but abandons the 
assumption that the matter action is independent of the 
connection. Clearly, metric affine f{R) gravity is the 
most general of these theories and reduces to metric or 
Palatini f{R) gravity if further assumptions are made. 
In this section we will present the actions and field 
equations of all three versions of f{R) gravity and point 
out their difference. We will also clarify the physical 
meaning behind the assumptions that discriminate 



5 



them. 

For an introduction to metric f{R) gravity see 
also ( Nojiri and Odintsovl l2007al ). for a shorter 
review of metric and Pa la tini f(R) gravity see 
( Capozziello and Francaviglial l2008l ) and for an ex- 
tensive analysis of all versions of f (R) grav i ty and other 
alternative theories of gravity see (jSotirio ul. l2007bD . 



A. Metric formalism 



Beginning from the action ^ and adding a matter 
term 5a/, the total action for f{R) gravity takes the form 



Sri 



-9/(«) + 5j/(s„^.*), (5) 



where ip collectively denotes the matter fields. Variation 
with respect to the metric gives, after some manipula- 
tions and modulo surface terms 

(6) 



where, as usual, 



T — 



-2 5S 



M 



-9 Sg^-' ' 



(7) 



a prime denotes differentiation with respect to the ar- 
gument, is the covariant derivative associated with 
the Levi-Civita connection of the metric, and □ = 
V^ Vu. Metric .f{R ) gravity was first rigorously studied 



in (|Buchdahi[l970l) .3 

It has to be stressed that there is a mathematical jump 
in deriving eq. ([6]) from the action ([5]) having to do with 
the surface terms that appear in the variation: as in the 
case of the Einstein-Hilbert action, the surface terms do 
not vanish just by fixing the metric on the boundary. 
For the Einstein-Hilbert action, however, these terms 
gather into a total variation of a quantity. Therefore, 
it is possible to add a total divergence to the action in 
order to "heal" it and arrive to a well-defined variational 
principle (this is the well known Gibbons-Hawking- Yo rk 
surface term (jGibbons and Hawking Il977t lYorkl Il972[ ) ) . 
Unfortunately, the surface terms in the variation of the 
action ([3]) do not consist of a total variation of some 
quantity (the reader is urged to calculate the variation in 
order to verify this fact) and it is not possible to "heal" 
the action by just subtracting some surface term before 
performing the variation. 

The way out comes from the fact that the action in- 
cludes higher order derivatives of the metric and, there- 
fore, it should be possible to fix more degrees of freedom 
on the boundary than those of the metric itself. There is 



^ Specific attention to higher-dimensional f(R) gravity was paid i n 
llGunther et all |2002| . |2003| . |2005|: ISaidov and Zhukl. I2006l.l2007t) . 



no unique prescription for such a fixing in the literature so 
far. Note also that the choice of fixing is not void of phys- 
ical meaning, since it will be relevant for the Hamiltonian 
formulation of the theory. However, the field equations 
(|6]) would be unaffected by the fixing chosen and from a 
purely classical perspective, such as the one follo wed here , 
the fie ld equations are all that one needs [see ( Sotirioul . 
l2007bD for a more detailed discussion on these issues!. 

Setting aside the complications of the variation we can 
now focus on the field equations (|6]). These are obviously 
fourth order partial differential equations in the metric, 
since R already includes second derivatives of the latter. 
For an action which is linear in R, the fourth order terms 
— the last two on the left hand side — vanish and the 
theory reduces to GR. 

Notice also that the trace of eq. ^ 



f{R)R~2f{R) + 3nf = kT, 



(8) 



where T = g^'^Tf^^, relates R with T differentially and 
not algebraically as in GR, where R = —kT. This is al- 
ready an indication that the field equations of f{R) theo- 
ries will admit a larger variety of solutions than Einstein's 
theory. As an example, we mention here that the Jebsen- 
Birkhoff's theorem, stating that the Schwarzschild solu- 
tion is the unique spherically symmetric vacuum solution, 
no longer holds in metric f{R) gravity. Without going 
into details, let us stress that T = no longer implies 
that i? = 0, or is even constant. 

Eq. ([5]) will turn out to be very useful in studying var- 
ious aspects of f{R) gravity, notably its stability and 
weak-field limit. For the moment, let us use it to make 
some remarks about maximally symmetric solutions. Re- 
call that maximally symmetric solutions lead to a con- 
stant Ricci scalar. For R = constant and T^^ = 0, eq. dS]) 
reduces to 



/'(i?)i? - 2/(i?) - 0, 



(9) 



which, for a given /, is an algebraic equation in i?. If 
i? = is a root of this equation and one takes this root, 
then eq. © reduces to R^^i, = and the maximally sym- 
metric solution is Minkowski spacetime. On the other 
hand, if the root of eq. ^ is R = C, where C is a con- 
stant, then eq. ^ reduces to Rf^^ — g^uC/4: and the 
maximally symmetric solution is de Sitter or anti-de Sit- 
ter space depending on the sign of C, just as in GR with 
a cosmological constant. 

Another issue that should be stressed is that of en- 
ergy conservation. In metric f{R) gravity the matter 
is minimally coupled to the metric. One can, therefore, 
use the usual arguments based on the invariance of the 
action under diffeomorphisms of the spacetime manifold 
[coordinate transformations — >■ x'^ — x'^ -I- fol- 
lowed by a puUback, with the field vanishing on the 
boundary of the spacet ime region c onsidered, leave the 
physics unchanged, see ( Waldl . [l98^ ] to show that T^^ is 
divergence-free. The same can be done at the level of the 
field equations: a "brute force" calculation reveals that 



6 



the left hand side of eq. ([6]) is divergence-free (g eneralized 
Bianch i identity) implying that V^T^'' = (|Koivistol . 
[20063)4 

Finally, let us note that it is possible to write the field 
equations in the form of Einstein equations with an ef- 
fective stress-energy tensor composed of curvature terms 
moved to the right hand side. This approach is question- 
able in principle (the theory is not Einstein's theory and 
it is artificial to force upon it an interpretation in terms 
of Einstein equations) but, in practice, it has been proved 
to be useful in scalar-tensor gravity. Specifically, eq. dU 
can be written as 



the metric and an independent connection (Palatini vari- 
ation) . The action is formally the same but now the Rie- 
mann tensor and the Ricci tensor are constructed with 
the independent connection. Note that the metric is not 
needed to obtain the latter from the former. For clarity 
of notation, we denote the Ricci tensor constructed with 
this independent connection as TZf^i, and the correspond- 
ing Ricci scalar^ is TZ = g'^'"R-^v The action now takes 
the form 



Spai ^ y d^^V^fC^) + SA^{g^,,y,'4')■ 



f'{R) '^'^^ 2f'{R) 
f'{R) 



[f{R)-Rf{R)] 



G - ^ (t + T^^'f 

^ f'{R) v^^ f^" 



•If) 



(10) 



(11) 



where the quantity Ge// = G/f'{R) can be regarded as 
the effective gravitational coupling strength in analogy 
to what is done in scalar-tensor gravity — positivity of 
Gef f (equivalent to the requirement that the graviton is 
not a ghost) imposes that f'{R) > 0. Moreover, 



llf{R)-Rf'iR) 



-g^,n/'(i?) 



5p. + V^V,/'(i?) 



(12) 



is an effective stress-energy tensor which does not have 
the canonical form quadratic in the first derivatives of 
the field f'{R), but contains terms linear in the second 
derivatives. The effective energy density derived from 
it is not positive-definite and none of the energy condi- 
tions holds. Again, this situation is analogous to that 
occurring in scalar-tensor gravity. The effective stress- 
energy tensor (fT2|) can be put in the form of a perfect 
fluid energy-momentum tensor, which will turn out to be 
useful in Sec. ITVl 



B. Palatini formalism 

We have already mentioned that the Einstein equa- 
tions can be derived using, instead of the standard met- 
ric variation of the Einstein-Hilbert action, the Palatini 
formalism, i. e. , an independent variation with respect to 



GR will come about, as we will see shortly, when f{TZ) = 
TZ. Note that the matter action Sm is assumed to de- 
pend only on the metric and the matter fields and not on 
the independent connection. This assumption is crucial 
for the derivation of Einstein's equations from the linear 
version of the action (|13p and is the main feature of the 
Palatini formalism. 

It has already been mentioned that this assumption has 
consequences for the p h ysical m ea ning of the independent 
connection (jSotirioul . l2006blldl : ISotiriou and Liberatil . 
l2007bi) . Let us elaborate on this: recall that an affine 
connection usually defines parallel transport and the co- 
variant derivative. On the other hand, the matter action 
S M is supposed to be a generally covariant scalar which 
includes derivatives of the matter fields. Therefore, these 
derivatives ought to be covariant derivatives for a general 
matter field. Exceptions exist, such as a scalar field, for 
which a covariant and a partial derivative coincide, and 
the electromagnetic field, for which one can write a co- 
variant action without the use of the covariant derivative 
[it is the exterior derivative that is actually nee ded, see 
next section and I Sotiriou and Liberatil . l200"7bl )]. How- 
ever, Sm should include all possible fields. Therefore, 
assuming that Sm is independent of the connection can 
imply one of two things (Sotiriou, 2006d): either we are 
restricting ourselves to specific fields, or we are implic- 
itly assuming that it is the Levi-Civita connection of the 
metric that actually defines parallel transport. Since the 
first option is implausibly limiting for a gravitational the- 
ory, we are left with the conclusion that the independent 
connection r-^^^, does not define parallel transport or the 
covariant derivative and the geometry is actually pseudo- 
Riemannian. The covariant derivative is actually defined 
by the Levi-Civita connection of the metric {^^i,}- 

This also implies that Palatini /(i?) gravity is a met- 
ric theo ry in the s ense that it satisfies the metric pos- 
tulates Il98l[ ). Let us clarify this: matter is mini- 
mally coupled to the metric and not coupled to any other 
fields. Once again, as in GR or metric f{R) gravity, 
one could use diffeomorphism invariance to show that 



^ The term "f{R) gravity" is used generically for a theory in which 
* Energy-momentum com plexes in the sphe r ically symmetric case the action is some function of some Ricci scalar, not necessarily 



-h/nergy-momentum com plexes m the sphe r ically symmetric case the 
have been computed in l lMultamaki et all . |200^ . R. 
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the stress energy tensor is conserved by the covariant 
derivative defined with the Levi-Civita connection of the 
metric, i.e., V^T^" = (but V^T^" ^ 0). This can 
also be shown by using the field equations, which we will 
present shortly, in order to calculate the divergence of T^^ 
with respect to the Levi- Civita connection of the metric 
and show that it vanishes ( Barraco et al ]. ll999HKoivistol 
l2006al ).^ Clearly then, Palatini f{R) g ravity is a metric 
theory according to the definition of (jWilt Il98l[ ) (not 
to be confused with the term "metric" in "metric f{R) 
gravity" , which simply refers to the fact that one only 
varies the action with respect to the metric). Conven- 
tionally thinking, as a consequence of the covariant con- 
servation of the matter energy-momentum tensor, test 
particles should follow geodesies of the metric in Palatini 
f{R) gravity. This can be seen by considering a dust 
fiuid with = pu^Uu and projecting the conserva- 
tion equation \/^T^p = onto the fluid four- velocity . 
Similarly, theories that satisfy the metric postulates are 
supp osed to sati sfy the Einstein Equivalence Principle as 
well (jWilll . lTgSlI ). Unfortunately, things are more compli- 
cated here and, therefore, we set this issue aside for the 
moment. We will return to it and attempt to fully clarify 
it in Sees ■ I VI. B] and fVI . C . 2 1 For now, let us proceed with 
our discussion of the field equations. 

Varying the action (|13|) independently with respect to 
the metric and the connection, respectively, and using 
the formula 



yields 



(14) 



(15) 



(16) 



where T^^ is defined in the usual way as in eq. ([7]), 
denotes the covariant derivative defined with the inde- 
pendent connection F'*'^^, and {^v) and [^v] denote sym- 
metrization or anti-symmetrization over the indices fi and 
ly, respectively. Taking the trace of eq. P^ . it can be 
easily shown that 



0, 



(17) 



which implies that we can bring the field equations into 
the more economical form 

/'(7^)7^(^,) - i/(7^).g,,, = kGT,,,, (18) 



0, 



(19) 



It is now easy to see how the Palatini formalism leads 
to GR when f{TZ) = TZ; in this case f'{TZ) = 1 and 
eq. (fT9|) becomes the definition of the Levi-Civita connec- 
tion for the initially independent connection F^^^^. Then, 
TZ^i, — Rfj,y, TZ — R and eq. (|T8l) yields Einstein's equa- 
tions. Thi s reproduces the result that can be found in 
textbooks (|Misner et all Il973l : IWaldl . Il984[ ). Note that 
in the Palatini formalism for GR, the fact that the con- 
nection turns out to be the Levi-Civita one is a dynamical 
feature instead of an a priori assumption. 

It is now evident that generalizing the action to be a 
general function of TZ in the Palatini formalism is just 
as natural as it is to generalize the Einstein-Hilbert 
action in the metric formalism.^ Remarkably, even 
though the two formalisms give the same results for 
linear actions, they lead to difFer ei it results for more 
genera l actions dB uchda hl,, 1970l: Burton and Mann . 
'1998a', b; Exirifard and Sheikh- Jabbaril l2008t iQuerella . 
11998; ,Shahid-Saless. 1987lv 

Finally, let us present some useful manipulations of the 
field equations. Taking the trace of eq. (|18p yields 



/'(7^)7^ - 2/(7^) = kT. 



(20) 



As in the metric case, this equation will prove very useful 
later on. For a given /, it is an algebraic equation in 
TZ. For all cases in which T = 0, including vacuum and 
electrovacuum, TZ will therefore be a constant and a root 
of the equation 



f{TZ)TZ - 2/(7^) = 0. 



(21) 



We will not consider cases for which this equation has 
no roots since it can be shown that the fi eld equations 
are then inconsistent ( Ferraris et all Il992l ). Therefore, 
choices of / that lead to this behaviour should simply 
be avoided. Eq. (PT|) can also be identically satisfied if 
f{TZ) cx TZ"^. This very partic ular choice for / lead s to 
a conformally invariant theory ()Ferraris et aZ.I . [19921 ). As 
is apparent from eq. ([20|) . if f{Ti) oc TZ^ then only con- 
formally invariant matter, for which T = identically, 
can be coupled to gravity. Matter is not generically con- 
formally invariant though, and so this particular choice 
of / is not suitable for a low energy theory o f gravity . 
We wi ll, therefore, not consider it further [see ( Sotirioul , 
l2006bD for a discussion! . 

Next, we consider eq. (IT^ . Let us define a metric 
conformal to g^^i, as 



hf,,. = f'{TZ)g^u- 



(22) 



' Energy supertensors and pseudotensora i n Palatini f{K) grav- 
ity were studied in llBarraco et alV Il999l : iBorowiec et al\ . 1 19941 . 
|199^: IFerraris ali Il992h and alternative energy definitions 
were given in l lDeser and Tekinl . |2002| . l2003al lbl [SoOTI) . 



See, however, llSotirioul. l2007bD for further analysis of the f(R) 
action and how it can be derived from first principles in the two 
formalisms. 
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It can easily be shown that® 



(23) 



Then, eq. becomes the definition of the Levi- 

Civita connection of h^^, and can be solved algebraically 
to give 



1 



or, equivalently, in terms of g^^, 
1 1 



(24) 



2/'(7^)■ 

-9,(/'(7^)g^,) 



5p(/'(7^).9..) + a.(/'(7^)5^.) 



(25) 



Given that eq. (1^01) relates 7?, algebraically with T, and 
since we have an explicit expression for F^^^, in terms 
of TZ and g'^'^, we can in principle eliminate the inde- 
pendent connection from the field equations and express 
them only in terms of the metric and the matter fields. 
Actually, the fact that we can algebraically express T'*'^^ 
in terms of the latter two already indicates that this con- 
nections act as some sort of auxiliary field. We will ex- 
plore this further in Sec. IIIII For the moment, let us 
take into account how the Ricci tensor transforms under 
conformal transformations and write 



1 



2 (/'(7^))2 
1 



(v^/'(7^))(v./'(7^)) 



fin) y^^^" ~ 2-^^''° ) ^^^^ 

Contraction with g'^'^ yields 

n^R + ^^^,^^(v^/'(7^))(v^/'(7^)) 

Note the difference between TZ and the Ricci scalar of h^i, 
due to the fact that g^^, is used here for the contraction 
of TZ^u- 

Replacing eqs. (PE)) and (P7|) in eq. p^ . and after some 
easy manipulations, one obtains 



1 

3 1 
~2F 



(V,/')(V./')--5p.(V/')' 



(28) 



Notice that, assuming that we know the root of eq. (|20l) . 
TZ = TZ{T), we have completely eliminated the indepen- 
dent connection from this equation. Therefore, we have 
successfully reduced the number of field equations to one 
and at the same time both sides of eq. p8| depend only 
on the metric and the matter fields. In a sense, the the- 
ory has been brought to the form of GR with a modified 
source. 

We can now straightforwardly deduce the following: 

• When f{TZ) = TZ, the theory reduces to GR, as 
discussed previously. 



• For matter fields with T — 0, due to eq. (|2T|) . TZ and 
consequently f{TZ) and f'{TZ) are constants and the 
theory reduces to GR with a cosmological constant 
and a modified coupling constant G/f. If we de- 
note the value of TZ when T = as TZo, then the 
value of the cosmological constant is 



TZo 



/(7^o) 

/'(7^o) 



TZo 

4 ' 



(29) 



where we have used eq. (1211) . Besides vacuum, T — 
also for electromagnetic fields, radiation, and any 
other conformally invariant type of matter. 

• In the general case T ^ 0, the modified source on 
the right hand side of eq. (1^51) includes derivatives 
of the stress-energy tensor, unlike in GR. These are 
implicit in the last two terms of eq. (1^51) , since /' is 
in practice a function of T, given that^ /' = f'{TZ) 
and TZ = TZ{T). 

The serious implications of this last observation will 
become clear in Sec. I VI. C. II 



C. Metric-affine formalism 

As we already pointed out, the Palatini formalism 
of f{R) gravity relies on the crucial assumption that 
the matter action does not depend on the independent 
connection. We also argued that this assumption rele- 
gates this connection to the role of some sort of aux- 
iliary field and the connection carrying the usual ge- 
ometrical meaning — parallel transport and definition 
of the covariant derivative — remains the Levi-Civita 
connection of the metric. All of these statements will 
be supported further in the forthcoming sections, but 
for the moment let us consider what would be the out- 
come if we decided to be faithful to the geometrical in- 
terpretation of the independent connection F"^^;^: this 



This calculation holds in four dimensions. When the num- 
ber of dimensions D is different from 4 then, instead of us- 
ing eq. II22II . the conformal metric /i^^ should be introduced as 
= [/'{7?.)]2/(o-2)g^^ in order for eq. ((23J to still hold. 



^ Note that, apart from special cases such as a perfect fluid, T^^ 
and consequently T already include first derivatives of the matter 
fields, given that the matter action has such a dependence. This 
implies that the right hand side of eq. II28I I vfill include at least 
second derivatives of the matter fields, and possibly up to third 
derivatives. 



would imply that we would define the covariant deriva- 
tives of the matter fields with this connection and, there- 
fore, we would have Sm = 'S'A/(5;^iy, F^^^, Vj)- The ac- 
tion of this theory, dubbed metric-affine /(i?) gravity 
( Sotiriou and Liberati I20d7bl) . would then be [note the 
difference with respect to the action (fT3|) ] 



(30) 



1. Preliminaries 

Before going further and deriving field equations from 
this action certain issues need to be clarified. First, since 
now the matter action depends on the connection, we 
should define a quantity representing the variation of Sm 
with respect to the connection, which mimics the defi- 
nition of the stress-energy tensor. W e call this quantity 
the h ypermomentum and is defined as ( Hehl and Kerlina 
I1978D 



5S 



M 



-g 



(31) 



Additionally, since the connection is now promoted to 
the role of a completely independent field, it is interesting 
to consider not placing any restrictions to it. Therefore, 
besides dropping the assumption that the connection is 
related to the metric, we will also drop the assumption 
that the connection is symmetric. It is useful to define 
the following quantities: the non-metricity tensor 



(32) 



which measures the failure of the connection to covari- 
antly conserve the metric, the trace of the non-metricity 
tensor with respect to its last two (symmetric) indices, 
which is called the Weyl vector. 



and the Cartan torsion tensor 



(33) 



(34) 



which is the antisymmetric part of the connection. 

By allowing a non-vanishing Cartan torsion tensor we 
are allowing the theory to naturally include torsion. Even 
though this brings complications, it has been consid- 
ered by some to be an advantage for a gravity the- 
ory since some matter fields, such as Dirac fields, can 
be coupled to gravit y in a way which might be con- 
sidered more natural (jHehl et aU [l995): one might ex- 
pect that at some intermediate or high energy regime, 
the spin of particles might interact with the geometry 
(in the same sense that macroscopic angular momen- 
tum interacts with geometry) and torsion can naturally 
arise. Theories with torsion have a long history, probably 



star ting with the Einstein-Cartan(-Sciama-Kibble) the- 



Cartanl. Il922[ | 1923L Il924t iHehl et al. l. ll976HKibbl3 . 
Sciama . 1964) . In this theory, as well as in other 



theories with an independent connection, some part of 
the connection is still related to the metric {e.g., the 
non-metricity is set to zero). In our case, the connec- 
tion is left completely unconstrained and is to be deter- 
mined by the field equations. Metric-affine gravity with 
the linea r version of the action (|30l) was initially pro- 
posed in (jHehl and Kerlintdl 19781 ) and the generalization 
to .f(l^ ) actions was considered in ( Sotiriou and Liberatil . 
l2007allbl ). 

Unfortunately, leaving the connection completely un- 
constrained comes with a complication. Let us consider 
the projective transformation 



(35) 



where is an arbitrary covariant vector field. One can 
easily show that the Ricci tensor will correspondingly 
transform like 



(36) 



However, given that the metric is symmetric, this implies 
that the curvature scalar does not change 



7^ ^- 7^, 



(37) 



i.e., TZ is invariant under projective transformations. 
Hence the Einstein-Hilbert action or any other action 
built from a function of TZ, such as the one used here, 
is projective invariant in metric-affine gravity. However, 
the matter action is not generically projective invariant 
and this would be the cause of an inconsistency in the 
field equations. 

One could try to avoid this problem by generalizing 
the gravitational action in order to break projective in- 
variance. This can be done in several ways, such as 
allowing for the metric to be non-symmetric as well, 
adding higher order curvature invaria nts or terms in- 
cluding the Cartan torsion t ensor [see ( Sotirioul . l2007bt 
ISotiriou and Liberatil . [2007bD for a more detailed discus- 
sion]. However, if one wants to stay within the framework 
of f{R) gravity, which is our subject here, then there is 
only one way to cure this problem: to somehow constrain 
the connection. In fact, it is evident from eq. pSI) that, 
if the connection were symmetric, projective invariance 
would be broken. However, one does not have to take 
such a drastic measure. 

To understand this issue further, we should re-examine 
the meaning of projective invariance. This is very similar 
to gauge invariance in electromagnetism (EM). It tells us 
that the corresponding field, in this case the connections 
F^^j,, can be determined from the field equations up to a 
projective transformation [eq. psp]. Breaking this invari- 
ance can therefore come by fixing some degrees of free- 
dom of the field, similarly to gauge fixing. The number 
of degrees of freedom which we need to fix is obviously 
the number of the components of the four-vector used 
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for the transformation, i.e., simply four. In practice, this 
means that wc should start by assuming that the con- 
nection is not the most general which one can construct, 
but satisfies some constraints. 

Since the degrees of freedom that we need to fix 
are four and seem to be related to the non-symmetric 
part of the connection, the most obvious prescription 
is to demand that = S„.^ b e equa l to zero, which 



was first suggested in (|Sandberd . flQTSi) for a linear ac- 
tion and shown to work als o for an /(T?,) action in 
( Sotiriou and Liberati l2007bD .^° Note that this does 



not mean that T should vanish, but merely that 
r^gf = T^^ . Imposing this constraint can easily be done 
by adding a Lagrange multiplier i?^. The additional term 
in the action will be 



d'^Xy/^B'^Sf,. 



(38) 



The action (|30p with the addition of the term in eq. ((38)) 
is, therefore, the action of the most general metric-affine 
f{R) theory of gravity. 



Taking the trace of eq. (|43p over the indices /i and A and 
using eq. (|44|) yields 



3 " 

Therefore, the final form of the field equations is 



-VA(^/^/'(W 



5, 



0. 



(45) 



(46) 



-2/'(7^).g^'^5,/ = k{A^^- - -Af^'-S^^h), (47) 



(48) 



Next, we examine the role of A^'^ . By splitting eq. p7)) 
into a symmetric and an antisymmetric part and per- 
formi ng contractions and manipula tions it can be shown 
that (|Sotiriou and Liberatil . l2Q07bl ) 



A 



0^S„ 



0. 



(49) 



2. Field Equations 

We are now ready to vary the action and obtain field 
equations. Due to space limitations, we will not present 
the various steps of the variation here. Instead we merely 
give the formula 



-1'"- "T 2r'^j^;^](5r'^^^, (39) 

which is useful to those wanting to repeat the variation 
as an exercise, and we also stress our definitions for the 
covariant derivative 



-pa AL> 



(40) 



and for the Ricci tensor of an independent connection 



(41) 



per 
rA fiu 



CTZ^ /^A • 



The outcome of varying independently with respect to 
the metric, the connection and the Lagrange multiplier 
is, respectively. 



^/'(7^)g'"')+V„ 



(42) 



+2/'(7^) {g'^'-S^f ~ g^PSpfS'^x + 5^"^.a'') 
= k(A/"-B1^5''1a), 



(43) 
(44) 



The proposal of l|HehI and Kerliii3 . Il978h to fix part of the non- 
metricity, namely the Weyl vector Q^, in ord er to break pro- 
jective invariance works only when /(7^) = l|Sotiriod . l2007bl : 
ISotiriou and Liberatil. l2007bl) . 



This straightforwardly implies two things: a) Any torsion 
is introduced by matter fields for which A^^''^'^^ is non- 
vanishing; b) torsion is not propagating, since it is given 
algebraically in terms of the matter fields through A_^''^'^^ . 
It can, therefore, only be detected in the presence of such 
matter fields. In the absence of the latter, spacetime will 
have no torsion. 

In a similar fashion, one can use the symmetrized ver- 
sion of eq. (|47l) to show that the symmetric part of the 
hypermomentum A-^^"^ is algebraically related to the 
non-metricity Qp^^x- Therefore, matter fields with non- 
vanishing A_^*-'"^'' will introduce non-metricity. However, 
in this case things are slightly more complicated because 
part of the non-metricity is also due to the functional 
form o f the Lagrangian itself [see ( Sotiriou and Liberatil . 
l2007bD ]. 

We will not perform a detailed study of different mat- 
ter fields and their role in metric-affine gravity. We 
refer the reader to the mo re exhaustive analysis of 
( Sotiriou and Liberatil . [2007bl ) for details and we restrict 
ourselves to the following remarks: Obviously, there are 
certain types of matter fields for which A^;^^"^ — 0. Char- 
acteristic examples are 

• A scalar field, since in this case the covariant 
derivative can be replaced with a partial deriva- 
tive. Therefore, the connection does not enter the 
matter action. 

• The electromagnetic field (and gauge fields in gen- 
eral), since the electromagnetic field tensor is 
defined in a covariant manner using the exterior 
derivative. This definition remains unaffected when 
torsion is incl uded [this can be related to ga uge in- 
variance, see ([Sotiriou and Liberatil l2007bl) for a 
discussion! . 
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On the contrary, particles with spin, such as Dirac fields, 
generically have a non-vanishing hypermomentum and 
will, therefore, introduce torsion. A more complicated 
case is that of a perfect fluid with vanishing vorticity. If 
we set torsion aside, or if we consider a fluid describing 
particles that would initially not introduce any torsion 
then, as for a usual perfect fluid in GR, the matter ac- 
tion can be written in terms of three scalars: t he energy 
density, the pressu re, and the velocity potential ( Schaket 
119961: IStonel . l2000l) . Therefore such a fluid wfll lead to a 
vanishing ^)^^ ■ However, complications arise when tor- 
sion is taken into account: Even though it can be argued 
that the spins of the individual particles composing the 
fluids will be randomly oriented, and therefore the expec- 
tation value for the spin should add up to z ero, fluctua- 
tions around this value will affect sp acetime ( Hehl et al\ . 
119761 : ISotiriou and Liberati l2007bli . Of course, such ef- 
fects will be largely suppressed, especially in situations 
in which the energy density is small, such as late time 
cosmology. 

It should be evident by now that, due to eq. (P^ . the 
field equations of metric f{R) gravity reduce to cqs. ([15]) 
and (|16p and, ultimately, to the field equations of Palatini 



f{R) gravity and ([TOl) . for all cases in which 
0. Consequently, in vacuo, where also T^^ = 0, they 
will reduce to the Einstein equations with an effective 
cosmological constant given by eq. (l29t . as discussed at 
the end of Sec. III. Bl for Palatini f{R) gravity. 

In conclusion, metric-affine f{R) gravity appears to 
be the most general case of /(i?) gravity. It includes 
enriched phenomenology, such as matter-induced non- 
metricity and torsion. It is worth stressing that torsion 
comes quite naturally, since it is actually introduced by 
particles with spin (excluding gauge fields). Remarkably, 
the theory reduces to GR in vacuo or for conformally in- 
variant types of matter, such as the electromagnetic field, 
and departs from GR in the same way that Palatini f{R) 
gravity does for most matter fields that are usually stud- 
ied as sources of gravity. However, at the same time, it 
exhibits new phenomenology in less studied cases, such 
as in the presence of Dirac fields, which include torsion 
and non-metricity. Finally, let us repeat once more that 
Palatini f{R) gravity, despite appearances, is really a 
metric theory according to the definition of (jWill . 1981!) 
(and the geometry is a priori pseudo-Riemannian).^^ On 
the contrary, metric-affine f{R) gravity is not a metric 
theory (hence the name). Consequently, it should also 
be clear that T^^'^ is not divergence-free with respect to 
the covariant derivative defined with the Levi-Civita con- 
nection (nor with actually). However, the physical 
meaning of this last statement is questionable and de- 
serves further analysis, since in metric-affine gravity T^j^ 



As mentioned in Sec. III.BI although the metric postulates are 
manifestly satisfied, there are ambiguities regarding the physical 
interpretation of this property and its relation with the Einstein 
Equivalence Principle (see Sec. IVI.C.ll l. 



does not really carry the usual meaning of a stress-energy 
tensor (for instance, it does not reduce to the special rel- 
ativistic tensor at an appropriate limit and at the same 
time there is also another quantity, the hypermomentum, 
which describes matter characteristics). 



III. EQUIVALENCE WITH BRANS-DICKE THEORY 
AND CLASSIFICATION OF THEORIES 

In the same way that one can make variable redefini- 
tions in classical mechanics in order to bring an equation 
describing a system to a more attractive, or easy to han- 
dle, form (and in a very similar way to changing coordi- 
nate systems), one can also perform field redefinitions in 
a field theory, in order to rewrite the action or the field 
equations. 

There is no unique prescription for redefining the fields 
of a theory. One can introduce auxiliary fields, perform 
renormalizations or conformal transformations, or even 
simply redefine fields to one's convenience. 

It is important to mention that, at least within a clas- 
sical perspective such as the one followed here, two the- 
ories are considered to be dynamically equivalent if, un- 
der a suitable redefinition of the gravitational and matter 
fields, one can make their field equations coincide. The 
same statement can be made at the level of the action. 
Dynamically equivalent theories give exactly the same 
results when describing a dynamical system which falls 
within the purview of these theories. There are clear ad- 
vantages in exploring the dynamical equivalence between 
theories: we can use results already derived for one the- 
ory in the study of another, equivalent, theory. 

The term "dynamical equivalence" can be considered 
misleading in classical gravity. Within a classical per- 
spective, a theory is fully described by a set of field 
equations. When we are referring to gravitation theories, 
these equations describe the dynamics of gravitating sys- 
tems. Therefore, two dynamically equivalent theories are 
actually just different representations of the same theory 
(which also makes it clear that all allowed representations 
can be used on an equal footing). 

The issue of distinguishing between truly different the- 
ories and different representations of the same theory (or 
dynamically equivalent theories) is an intricate one. It 
has serious implications and has been the cause of many 
misconceptions in the past, especially when conformal 
transformations are used in order to redefine the fields 
{e.g., the Jordan and Einstein frames in scalar-tensor the- 
ory). It goes beyond the scope of this review to present 
a detailed analysis of this issue. We refer the reader to 
the literature, and specifically to ( Sotiriou et al\ . |2007| ) 
and references therein for a detailed discussion. Here, we 
simply mention that, given that they are handled care- 
fully, field redefinitions and different representations of 
the same theory are perfectly legitimate and constitute 
very useful tools for understanding gravitational theories. 

In what follows, we review the equivalence between 
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metric and Palatini f{R) gravity with specific theories 
within the Brans-Dicke class with a potential. It is shown 
that these versions of f{R) gravity are nothing but dif- 
ferent representations of Brans-Dicke theory with Brans- 
Dicke parameter ujq — and ujq — —3/2, respectively. We 
comment on this equivalence and on whether preference 
to a specific representation should be an issue. Finally, 
we use this equivalence to perform a classification of f{R) 
gravity. 



A. Metric formalism 

It has been noticed quite early that metric quadratic 
gravity can be cast into the form of a Brans-Dicke theory 
and it did not take long for these results to be extended 
to more general act ions which are functions of the Ricci 
scalar of the metric ('Barrow', '1988'; 'Barrow and Cotsakisl 
119881: iTevssandier an d Tourrcnc, 1983; Wands, 1994) Jsee 
also (|FlanaganL l2'003l ) and (jCecottil . Il987l : IWandd . Il99l 
for the extension to theories of the type f{R, O^R) with 
fc > 1 of interest in supergravity] . This equivalence has 
been re-examined recent ly due to the increased i nterest 
in metri c f(R) gravity (|Chibal . 120031 : iFlanaganl l2004al: 
ISotiriouL[2006bh . Let us present this equivalence in some 
detail. 

We will work at the level of the action but the same 
approach can be used to work directly at the level of the 
field equations. We begin with metric f{R) gravity. For 
the convenience of the reader, we rewrite here the action 



©: 



1 

2k 



(50) 



One can introduce a new field x and write the dynami- 
cally equivalent action 



Srt 



1 

2k 



'9[fix) + f'{x){R-x)] 

+Sm{9,iv,^)- 
Variation with respect to x leads to the equation 

/"(x)(i?-x)-o. 



(51) 



(52) 



Therefore, x = ^ if /"(x) 7^ Oi which reproduces the 
action Redefining the field x by = /'(x) ^nd 

setting 

n0)=x(0)0-/(x(0)), (53) 
the action takes the form 
1 



S^et = — / y(,^)] +5m(5m-,V'). (54) 



The action is sometimes cal led "R-regular" by mathematic al 
physicists if f"{R) ^ [e.g., jMagnano and Sokolowskil[l994l 'l]. 



This is the Jordan frame representation of the action 
of a Brans-Dicke theory with Brans-Dicke parameter 
LjQ — 0. An LjQ — Brans-Dicke theory Jsome- 
times called "massive dilaton g ravity" (IW andsl . Il994[) ] 
was originally proposed by (|0'IIanloni . il972br ) in or- 
der to generate a Yukawa term in the Newtonian limit 
and has b e en o ccasionally considere d in the literature 
(lAndersonl Il971, : Barber. 2003 : Dabrowski al 



1972al : lO'Hanlon and Tupper 



1972) 



2007 



Davidsqnl 120051: iDeseil. Il970l: iFuiiil . Il982l: lO'Hanlon 



It 



should be 

= fix) is 
like all 



stressed that the scalar degree of freedom (f> — 
quite different from a matter field; for example 
nonminimally cou pled scalars, it can violate all of the 
energy conditions ( Faraonil l2004a|) . 

The field equations corresponding to the action ([M)) 
are 

K 1 

+ ^{^^N.(p-9^..a(j)), (55) 

R = V'{(j)). (56) 

These field equations could have been derived directly 
from eq. ([6]) using the same field redefinitions that were 
mentioned above for the action. By taking the trace of 
eq. ([55t in order to replace R in eq. ([551) . one gets 



3n(f> + 2V{(l)) - (j) 



dV 



kT. 



(57) 



This last equation determines the dynamics of for given 
matter sources. 

The condition /" for the scalar-tensor theory to be 
equivalent to the original f{R) gravity theory can be seen 
as the condition that the change of variable cj) — f'{R) 
needed to express the theory as a Brans-Dicke one ([5i)) 
be invertible, i.e., d(j)/dR = /" ^ 0. This is a suffi- 
cient but not necessary condition for invertibility: it is 
only necessar y that f'{R) be continuous and one-to-one 
()01mol |2007[ ). By looking at eq. it is seen that 

/" implies (j) = f'{R) and the equivalence of the ac- 
tions ^ and ([^ . When /" is not defined, or it vanishes, 
the equality (j) — f'{R) and the equivalence between the 
two theories can not be guaranteed (although this it is 
not a priori excluded by /" = 0). 

Finally, let us mention that, as usual in Brans-Dicke 
theory and more general scalar-tensor theories, one can 
perform a conformal transformation and rewrite the ac- 
tion in what is called the Einstein frame (as opposed 
to the Jordan frame). Specifically, by performing the con- 
formal transformation 

gi-iu g^Lu = f'{R) g^iy = 4>gi-iu (58) 

and the scalar field redefinition ((> = f'{R) (p with 



2wo + 3 



2k 



(59) 
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a scalar-tensor theory is mapped into the Einstein frame 
in which the "new" scalar field (j) couples minimally to 
the Ricci curvature and has canonical kinetic energy, as 
described by the gravitational action 



(60) 



For the loq — equivalent of metric f{R) gravity we have 
= /'(i?) = eV¥^, (61) 



Even though the gravitational part of this action is for- 
mally the same as that of the action this action is 
not a Brans-Dicke one with luq — 0'- TZ is not the Ricci 
scalar of the metric g^i^. However, we have already seen 
that the field equation (ITS)) can be solved algebraically 
for the independent connection yielding eq. (1251) . This 
implies that we can replace the connection in the ac- 
tion without affecting the dynamics of the theory (the 
independent connection is practically an auxiliary field). 
Alternatively, we can directly use eq. (I?7)) . which relates 
R and TZ. Therefore, the action (p5|) can be rewritten, 
modulo surface terms, as 



Rf'iR) - f{R) 



2K{r{R)Y 

where R = R{(j)), and the complete action is 



(62) 



^-\d"~4>d^~^-U{4>) 



^Sm-,^'). (63) 



A direct transformation to the Einstein frame, without 
the intermediate p assage from the Jordan frame, has 
been discovered in ( Barrow and Cotsakisl Il988l : IWhittl . 
[1984). 

We stress once more that the actions ([5]) , (|54p , and ((63)) 
are nothing but different representations of the same 
theory. Additionally, there is nothing exceptional 
about the Jordan or the Einstein frame of the Brans- 
Dicke representation, an d one can actually find infinitely 
many conformal frames (|Flanaganl . r2004aHSotiriou tt all 
l2007t) . 



B. Palatini formalism 

Palatini /(i?) gravity can also be ca st in the form of a 
Brans-Dicke theory w i th a po tential ( Flanaganl l2004bl : 
lOlmol l2005bl: ISotirioul . l2006bl) . As a matter of fact, be- 
ginning from the Palatini /(i?) action, which we repeat 
here for the reader's convenience 



Spal 



2k. 



(64) 



and following exactly the same steps as before, i.e., in- 
troducing a scalar field x which we later redefine in terms 
of (p, yields 

Spal = ^ / d^xV^ [4>n - V {(!>)] + SM{g^.u, ^). (65) 



This has been an issue of debate and confusi on, see for example 
the references in l lFaraoni and Nadeaul . [2007^ . 



1 

2k 



d Xy/—g (f)R 



3 

20 



+SM{g^,v,'4^)■ 



(66) 



This is the action of a Brans-Dicke theory with Brans- 
Dicke parameter wg = —3/2. The corresponding field 
equations obtained from the action ()66|) through varia- 
tion with respect to the metric and the scalar are 



1 



(v^ - 5Mi^n0) 



V 
20 



□0 



:{R~V') 



20 



V^0V, 



(67) 
(68) 



Once again, we can use the trace of eq. ()67|) in order to 
eliminate R in eq. ([S5)) and relate directly to the matter 
sources. The outcome is 



2V ~(j)V' = K T. 



(69) 



Finally, one can also perform the conformal transfor- 
mation (1551) in order to rewrite the action (p5|) in the 
Einstein frame. The result is 



(70) 

where J7(0) = y(0)/(2/t0^). Note that here we have not 
used any redefinition for the scalar. 

To conclude, we have established that Palatini f{R) 
gravity can be cast into the form of an cjq = —3/2 Brans- 
Dicke theory with a potential. 



C. Classification 

The scope of this section is to present a classification 
of the different versions of f{R) gravity. However, before 
doing so, some remarks are in order. 

Let us, first of all, use the Brans-Dicke representation 
of both metric and Palatini f{R) gravity to comment 
on the dynamics of these theories. This representation 
makes it transparent that metric f{R) gravity has just 



14 



one extra scalar degree of freedom with respect to GR. 
The absence of a kinetic term for the scalar in the ac- 
tion (f54| or in eq. (|56|) should not mislead us to think 
that this degree of freedom does not carry dynamics. 
As can be seen by eq. ((57)) . (/) is dynamically related to 
the matter fields and, therefore, it is a dynamical de- 
gree of freedom. Of course, one should also not fail to 
mention that eq. (j56l) does constrain the dynamics of (j). 
In this sense metric f{R) gravity and ujq — Brans- 
Dicke theory differs from the general Brans-Dicke theo- 
ries and constitutes a special case. On the other hand, 
in the ujf) = —3/2 case which corresponds to Palatini 
f{R) gravity, the scalar appears to have dynamics in 
the action (|66|) or in eq. (|68l). However, once again this is 
misleading since, as is clear from eq. ((69|). </> is in this case 
algebraically related to the matter and, therefore, carries 
no dynamics of its own [indeed the field eqs. ((67|) and ((69)) 
could be combined to give eq. (j28p . eliminating (p com- 
pletely]. As a remark, let us state that the equivalence 
between Palatini /(i?) gravity and wq — —3/2 Brans- 
Dicke theory and the clarifications just made highlight 
two issues already mentioned: the fact that Palatini f{R) 
gravity is a metric theory according to the definition of 
( Willi . Tl98ll ). and the fact that the independent connec- 



tion is actually some sort of auxiliary field. 

The fact that the dynamics of (j) are not transparent at 
the level of the action in both cases should not come as 
a big surprise: (f> is coupled to the derivatives of the met- 
ric (through the coupling with R) and, therefore, partial 
integrations to "free" 5<j} or 5g^'^ during the variation are 
bound to generate dynamical terms even if they are not 



J 



initially present in the action. The wq = —3/2 case is 
even more intricate because the dynamical terms gener- 
ated through this procedure exactly cancel the existing 
one in the action. 

We already saw an example of how different represen- 
tations of the theory can highlight some of its charac- 
teristics and be very useful for our understanding of it. 
The equivalence between f{R) gravity and Brans-Dicke 
theory will turn out to be very useful in the forthcoming 
sections. 

Until now we have not discussed any possible equiva- 
lence between Brans-Dicke theory and metric-affine f{R) 
gravity. However, it is quite straightforward to see that 
there cannot be any. Metric- afhne /(i?) gravity is not a 
metric theory and, consequently, it can not be cast into 
the form of one, such as Brans-Dicke theory. For the 
sake of clarity, let us state that one could still start from 
the action ([50]) and follow the steps of the previous sec- 
tion to bring its gravitational part into the form of the 
action (j66p . However, the matter action would have an 
explicit dependence from the connection. Additionally, 
one would not be able to use eq. (I?7l) to eliminate TZ in 
favour of R since this only holds in Palatini f{R) gravity. 

In conclusion, metric-afHne f{R) gravity is the most 
general case of f[R) gravity. Imposing further assump- 
tions can lead to both metric or Palatini f{R) gravity, 
which can be cast into the form of wq = and = —3/2 
Brans-Dicke theories with a potential. In both cases, re- 
stricting the functional form of the action leads to GR. 
These results are summarized in the schematic diagram 

of Fig. m 




f{R) GRAVITY 
F"^^^ and independent 
METRIC-AFFINE f{R) 

Sm = SM{gfj.v,ip) 

PALATINI f{R) METRIC f{R) 

f{R) = R fiR) = R 

f"{R) + "^^"^^^ GR r\R) ^ 




BRANS-DICKE, 



BRANS-DICKE, = 



FIG. 1 Classification of f{R) theories of gravity and equivalent Brans-Dicke theories. The flowchart shows the list of assump- 
tions that are needed to arrive to the various versions of ,f{R) gravity and GR beginning from the the general f{R) action. It 
also includes the equivalent Brans-Dicke classes. Taken from (jSotirioul . l2006bl l. 



D. Why f{R) gravity then? tation, one might be led to ask two questions: first, why 

should we consider the f{R) representation and not just 

Since f{R) gravity in both the metric and the Palatini 
formalisms can acquire a Brans-Dicke theory represen- 
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work with the Brans-Dicke one, and second, why, since 
we know a lot about Brans-Dicke theory, should we re- 
gard f{R) gravity as unexplored or interesting? 

The answer to the first question is quite straightfor- 
ward. There is actually no reason to prefer either of the 
two representations — at least as far as classical grav- 
ity is concerned. There can be applications where the 
f{R) representation can be more convenient and applica- 
tions where the Brans-Dicke representation is more con- 
venient. One should probably mention that habit affects 
our taste and, therefore, an f{R) representation seems 
more appealing to relativists due to its more apparent 
geometrical nature, whereas the Brans-Dicke represen- 
tation seems more appealing to particle physicists. This 
issue can have theoretical implications. To give an exam- 
ple: if f{R) gravity is considered as a step towards a more 
complicated theory, which generalisation would be more 
straightfo rward will depend on the chosen representation 
[see also (jSotiriou et all |2007() for a discussion] . 

Whether f{R) theories of gravity are unexplored and 
interesting or just an already-studied subcase of Brans- 
Dicke theory, is a more practical question that certainly 
deserves a direct answer. It is indeed true that scalar- 
tensor theories and, more precisely, Brans-Dicke the- 
ory are well-studied theories which have been extensively 
used in many applications, including cosmology. How- 
ever, the specific choices ujq = 0,-3/2 for the Brans- 
Dicke parameter are quite exceptional, as already men- 
tioned in the previous section. It is also worthwhile point- 
ing out the following: a) As far as the luo — case is 
concerned, one can probably speculate that it is the ap- 
parent absence of the kinetic term for the scalar in the 
action which did not seem appealing and prevented the 
study of this theory, b) The ujq = —3/2 case leads to a 
conform ally invariant the ory in the absence of the poten- 
tial [see (ISotirioul . [2006bl) and references therein], which 
constituted the initial form of Brans-Dicke theory, and 
hence it was considered non-viable (a coupling with non- 
conformally invariant matter is not feasible). However, in 
the presence of a potential, the theory no longer has this 
feature. Additionally, most calculations which are done 
for a general value of loq in the literature actually exclude 
loq = —3/2, mainly because, merely for simplicity pur- 
poses, they are done in such a way that the combination 
2ujo + 3 appears in a denominator (see also Sec. IV. A| . 

In any case, the conclusion is that the theories in the 
Brans-Dicke class that correspond to metric and Palatini 
f{R) gravity had not yet been explored before the recent 
re- introduction of f{R) gravity and, as will also become 
clear later, several of their special characteristics when 
compared with more standard Brans-Dicke theories were 
revealed through studies of f{R) gravity. 



IV. COSMOLOGICAL EVOLUTION AND CONSTRAINTS 

We now turn our attention to cosmology, which moti- 
vated the recent surge of interest in f{R) gravity in order 



to explain the current cosmic acceleration without the 
need for dark energy. Before reviewing how f{R) gravity 
might provide a solution to the more recent cosmological 
riddles, let us stress that the following criteria must be 
satisfied in order for an f{R) model to be theoretically 
consistent and compatible with cosmological observations 
and experiments. The model must: 

• have the correct cosmological dynamics; 

• exhibit the correct behaviour of gravitational per- 
turbations; 

• generate cosmological perturbations compatible 
with the cosmological constraints from the cosmic 
microwave background, large scale structure. Big 
Bang Nucleosynthesis, and gravity waves. 

These are independent requirements to be studied sep- 
arately, and they must all be satisfied. 



A. Background evolution 

In cosmology, the identification of our universe with a 
Friedmann-Lemaitre-Robertson- Walker (FLRW) space- 
time is largely based on the high degree of isotropy 
measured in the cosmic microwave background; this 
identification relies on a formal re sult known as the 
Ehlers-Geren-Sachs (EGS) theorem (|Ehlers et aUllQeSl ) 
which is a kinematical characterization of FLRW spaces 
stating that, if a congruence of timelike freely falling 
observers see an isotropic radiation field, then (as- 
suming that isotropy holds about every spatial point) 
the spacetime is spatially homogeneous and isotropic 
and, therefore, a FLRW one. This applies to a uni- 
verse filled wit h any perfect fluid that is geodesic 
and barotropic ( Clarkson and Barrettl . Il999l : lEllis et all . 
il983al ib). Moreover, an "almost-EGS theorem" holds: 
spacetimes that are close to satisfying the EGS condi- 
tions are close t o FLR W universes in an appropriate sense 
([Stoeger et all Il995[ ). One would expect that the EGS 
theorem be extended to f{R) gravity; indeed, its validity 
for the (metric) theory 



S 



1 

2k 



fiR^uR^'^] + Sm (71) 



was proved in (iMaartens and Tavloii 11994 
iTavlor and Maartensl Il995l ) and the gene ralization to 
arbit rary metric f{R) gravity was given by faipp l et al\ . 
Il996l ). The validity of the EGS theorem can also be seen 
through the equivalence between f{R) and Brans-Dicke 
theory: the the orem was extended to scalar-tensor 
theories in fClark son et~ai\ . 2001, 20"ol. Since metric 
and Palatini f{R) gravities are equivalent to w = 
and Wo = —3/2 Brans-D icke theories re s pectively, it 
seems that the results of ( Clarkson et al 
can be considered as straightforward generalizations 
of the EGS theorem in both versions of f{R) gravity 
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as well. However, in the case of Palatini f{R) gravity 
there is still some doubt regarding this is sue due to 
complications in averaging ( Flanagam l2004b( ) . 



1. Metric f{R) gravity 

Considering the discussion above, it is valid to use the 
FLRW line element 



-df + a^{t) 



dr^ 



1 — kr'^ 



(72) 

as a local description of spacetime at cosmological scales, 
where (<, r, 0, (p) are comoving coordinates. We remind 
the reader that k — —1,0,1 according to whether the 
universe is hyperspherical, spatially flat, or hyperbolic 
and that a{t) is called the scale factor. Part of the stan- 
dard approach, which we follow here as well, is to use 
a perfect fluid description for matter with stress-energy 
tensor 



T'^" = (p + P) u'^u" + Pg 



(73) 



where u'^ denotes the four- velocity of an observer comov- 
ing with the fluid and p and P are the energy density and 
pressure of the fluid, respectively. 

Note that the value of k is an external parameter. As 
in many other works in the literature, in what follows we 
choose fc = 0, i.e., we focus on a spatially flat universe. 
This choice in made in order to simplify the equations and 
should be viewed sceptically. It is sometimes claimed in 
the literature that such a choice is favoured by the data. 
However, this is not entirely correct. Even though the 
data [e.g. (jSpergel et a/.l . l2007l )] indicate that the current 
value of rik is very close to zero, it should be stressed that 
this does not really reveal the value of k itself. Since 



(74) 



the current value of flk is sensitive to the current value 
of a{t), i.e. to the amount of expansion the universe has 
undergone after the Big Bang. A significant amount of 
expansion can easily drive flk very close to zero. The 
success of the inflationary paradigm is exactly that it 
explains the flatness problem — how did the universe 
become so flat — in a dynamical way, allowing us to 
avoid fine-tuning the parameter k (the value fc = is 
statistically exceptional) . 

The above having been said, choosing fc = for sim- 
plicity is not a dramatic departure from generality when 
it comes to late time cosmology. If it is viewed as an 
approximation and not as a choice of an initial condi- 
tion, then one can say that, since Vlk as inferred from 
observations is very close to zero at current times, the 
terms related to k will be subdominant in the Friedmann 
or generalised Friedmann equations and, therefore, one 
could choose to discard them by setting k = 0, with- 
out great loss of accuracy. In any case, results derived 



under the assumption that k = should be considered 
preliminary until the influence of the spatial curvature 
is precisely determined, since there are indications that 
even a very small value of may have an effect on them 
[see, for instance ( Clarkson et a/.l . [2007() ]. 

Returning to our discussion, inserting the flat FLRW 
metric in the field equations ([6|) and assuming that the 
stress-energy tensor is that of eq. (1751) yields 



H' = 



2H + 



K 

K 

'7 



P +{kff"' + 2HRJ" 
1 



-Rf' + lif-Rf) 



(75) 



(76) 



With some hindsight, we assume that /' > in or- 
der to have a positive effective gravitational coupling 
and f" > to avoid the Dolgov-Kawasaki instability 
( Dolgov and Kawasakil [2003al : iFaraonl l2006al) discussed 
in Sec. EH 

A significant part of the motivation for f{R) gravity 
is that it can lead to accelerated expansion without the 
need for dark energy (or an inflaton field). An easy way 
to see this is to define an effective energy density and 
pressure of the geometry as 



Peff 



eff 



Rf - f 3HRf" 



V 



/' 



+ 2HRf" + Rf" + \{f - Rf) 
T 



(77) 
, (78) 



where Pe// has to be non- negative in a spatially fiat 
FLRW spacetime, as follows from the inspection of 
eq. ([75]) in the limit p 0. Then, in vacuo eqs. ([75)1 
and ([75]) can take the form of the standard Friedmann 
equation 



a 
a 



6 



[Peff + 3Pe//] 



(79) 

(80) 



Hence, in vacuo the curvature correction can be viewed 
as an effective fluid. ^"^ 

The effective equation of state parameter Wef / of mod- 
ifled gravity can be expressed as 



"eff 



Peff _ R\r + 2HRf" + Rf" + i(/ - Rf) 

Peff 

(81) 



^^^^ - 3HRf" 



Note the following subtlety though: should we have included 
matter it would enter the Friedmann equations with a modified 
coupling k/J'. In general this effective fluid representation is 
used only for demonstrative purposes and should not be overes- 
timated or misinterpreted. 
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Since the denominator on the right hand side of eq. ((8T|) 
is strictly positive, the sign of Wgff is determined by its 
numerator. In general, for a metric f{R) model to mimic 
the de Sitter equation of state Weff = —1, it must be 



/"' _ RH -R 



(82) 



Let us also give two simple examples that can be found 
in the literature for demonstrative purposes and without 
considering their viability: First, one can consider the 
function / to be of the form f{R) oc i?". It is quite 
straightforward to calculate Weff as a function of n if 
the scale factor is assumed to be a generic power law 
a{t) = ao{t/to)°' (a general a(t) would l ead to a time 
dependent We//) (|Capozziello et all . I2003D . The result is 



6n^ - 7n - 1 



"eff 



Qrfl — 9n + 3 
for n 7^ 1 , and a is given is terms of n as 
-2n^ + 3n - 1 



(83) 



(84) 



A suitable choice of n can lead to a desired value for 
Weff - For instance, n = 2 yields Weff = — 1 and a — oo, 
as expected, considering that quadratic corrections to the 
Einstein-Hilbert Lagr angian were used in the well known 
Starobinsky inflation ( Starobinskvl . Il980l) . 

The second example which we will refer to is a model 
of the form f{R) = i?— /x^^^+^Vfi", wh ere ^ is a suitably 
chosen parameter ([Carroll et all l2004 l) . In this case, and 
once again if the scale factor is assumed to be a generic 
power law, Wpff can again be written as a function of n 
(ICarroll et an . \2()d4 : 



Weff 



2(n + 2) 



3(2n + l)(n+ 1) 



(85) 



The m ost typical model w ithin this class is that with 
n = 1 ( Carroll et all [200^ . in which case Weff — —2/3. 
Note that in this class of models, a positive n implies 
the presence of a term inversely proportional to R in the 
action, contrary to the situation for the i?" models. 

In terms of the quantity (piR) = f'{R) one can rewrite 
eq. ((HI) as 



Weff 

and 



-1+2- 



(0 - H^'^ 
R(l)~ f ~ QH4) 



(0-^0) 
30i72 



K (0 



Peff + Peff 



H<j) 



d 

di 




(86) 



(87) 



An exact de Sitter solution corresponds to (p ~ f"{R)R = 
0, or to = Ca{t) = Cqq e^"*, where C ^ is an integra- 
tion constant. However, the second solution for (j>{t) is 



not acceptable because it leads to f"{R)R = Caoe^"*, 
which is absurd because the left hand side is time- 
independent (for a de Sitter solution), while the right 
hand side depends on time. 

One could impose energy conditions for the effective 
stress-energy tensor (|12p of f{R) gravity. However, 
this is not very meaningful from the physical point of 
view since it is well known that effective stress-energy 
terms originating from the geometry by rewriting the 
field equations of alternative gravities as effective Ein- 
stein equations do, in g eneral, violate all the energy 
conditions [e.g., ([Faraoni . 2004a)]. Also, the concept 
of gravitational energy density is, anyway, ill-defined 
in GR and in all metric theories of gravity as a con- 
sequence of the Equivalence Principle. Moreover, the 
violation of the energy conditions makes it possible to 
have H > and bounc i ng uii iverses ([Carloni et aLl . [20061: 
iNovello and BerdiaffaL [200I . 

The field equations are clearly of fourth order in a{t). 
When matter is absent (a situation of interest in early 
time inflation or in a very late universe completely domi- 
nated by f{R) corrections), a{t) only appears in the com- 
bination H = d/a. Since the Hubble parameter _ff is a 
cosmological observable, it is convenient to adopt it as 
the (only) dynamical variable; then the field equations 
([75)1 and (fZS)) are of third order in H. This elimination 
of a is not possible when fc 7^ 0, or when a fluid with 
density p — p{a) is included in the picture. 

Regarding the dynamical field content of the the- 
ory, the fact that quadratic corrections to the 
Einstein-Hilbert action introduce a massive scalar field 
was n oted in (iBuchbinder et al 



1992; StcUc, 1978, 



1977[: [Strominged . [1984 lutivama and DeWitt> . .1962 : 
VilkoviskvLll992D : this applies to an y f(R) gravity t heory 
in the metric formal i sm [see, e.g., ( Ferraris et al\ . [l988l : 
iHindawi et aZI . [19961: [Olmol l2007l) ]. The metric tensor 
contains, in principle, various degrees of freedom: spin 2 
modes, and vector and scalar modes, which can all be 
massless or massive. In GR we find only the massless 
graviton but, when the action is allowed to include terms 
that depend on i?, Rf^^R^" , R^j.vpaR^"'"^ , other modes 
show up. In f{R) gravity, a massive scalar mode appears, 
which is evident in the equivalence with scalar-tensor the- 
ory (see Sec. IIII|) . As discussed in Sec. Illl.Ci the scalar 
field (f) ~ R is dynamical in the metric formalism and 
non-dynamical in the Palatini formalism. 



15 In (ISantos et all l2007f l. the Null Energy Condition and the 
Strong Energy Condition for metric f{R) gravity have been de- 
rived by using the Raychaudhuri equation and imposing that 
gravity be attractive, whereas for the Weak Energy Condition 
and the Dominant Energy Condition an effec tive stress-energy 
tenso r which includes the matter was used. In llPerez Bergliaffal . 
I2OO6I ). a different approach was followed, in which the standard 
energy conditions on matter where used in an attempt to con- 
strain f{R) gravity. 
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2. Palatini f{R) gravity 

As already mentioned, some concerns have been ex- 
pressed on whether the homogeneity approximation can 
justify the use of the FLRW metri c as a cosmologi - 
cal soluti on in Palatini f (R) gravity ( Flanaganl l2004bl ) 
see also (|Li et al\ . \200m ]. Therefore, even though it is 
standard practice in the literature to assume a FLRW 
background and a perfect fluid description for mat- 
ter when studying cosmology in Palati ni f(R) gravity 
\e.q. (Allcmandi et ai, 2004, 2005a; A marzguioui et ai . 
2006; M eng and Wand. 120031 I2004bl 120051: ISotiriml 



2006ai ia: iVolhckl l2003lV |. and we are going to review this 



approach here, the reader should approach it with some 
reasonable skepticism until this issue is clarified further. 

Under the assumptions that the spacetime is indeed 
described at cosmological scales by the FLRW metric, 
eq. ((7^ . that the stress-energy tensor of matter is that of 
eq. (1731) . and that fc = 0, easy manipulations reveal that 
the field eqs. ([T^ and yield the f ollowing modified 
Fricdmann equation [see for instance ( Meng and Wanj . 
12003; Sotir iou. 2 006^) ]: 



2/'] 6 /' 6/'' 



where the overdot denotes differentiation with respect to 
coordinate time. Note that when / is linear, /' = 1 
and, therefore, /' = 0. Taking into account eq. ([20| . one 
can easily show that in this case eq. (I55|) reduces to the 
standard Fricdmann equation. 

We will avoid representing the extra terms in eq. ([55)) 
with respect to the standard Fricdmann equation as a 
an effective stress energy density and pressure since, as 
it is not that difficult to see, the former equation does 
not really carry more dynamics than the latter. Indeed, 
assume as usual that the cosmological fluid is composed 
by pressureless dust (F„i = 0) and radiation (P^ = Pr/3) 
and p = Pm + Pp and P ^ Pm + Pr where pm, Pp and P,„, 
Pr denote the energy density and the pressure of dust and 
radiation, respectively. Due to eq. ([20)) and the fact that 
for radiation T = 0, it is quite straightforward to derive 
an algebraic relation between TZ and the energy density 
of the dust. Combining thi s with energy conservation, 
one obtains ( Sotirio ul. l2006eD 



n 



3H{nf' - 2/) 

nr - f ■ 



(89) 



This equation can be used to replace TZ in eq. ([55)) . yield- 
ing 



1 2Kp + nf-f 
6/' f. _ 3 f"(nf'-2f) 

2 f'{TZf"-f') 



(90) 



Considering now that, due to eq. ([^ . TZ is just an alge- 
braic function of pm, it is easy to realize that eq. ()90p is 



actually just the usual Fricdmann equation with a modi- 
fied source. The functional form of / will determine how 
the dynamics will be affected by this modification. 

It seems, therefore, quite intuitive that by tampering 
with the function / one can affect the cosmological dy- 
namics in a prescribed way. Indeed it has been shown 
that for f{TZ) = TZ — /{2>TZ) one approaches a de Sit- 
ter expansion as the density goes to zero ( Vollic^, ,2003.) . 
In order to match observations of the expansion history. 



one needs to choose a ~ 10 



-67 



eV)2 



10 



'53 , 



Ad- 



ditionally, in regimes for which np ^ a, eq. ()90|) reduces 
to high precision to the standard Fricdmann equation. 
The above can very easily be verified by replacing this 
particular choice of / in eq. ([90)) . We refer the reader to 
the literature for more details. 

One could, of course, consider more general functions 
of TZ. Of particular interest would be having posi- 
tive powers of TZ higher than the first power added in 
the action (since one could think of the Lagrangian as 
a series expansion). Indeed this has been consider ed 
([Meng and Wan 3, l2003l l2004bl 120051 : ISotiriod . [2006all3 ) . 
However, it can be shown that such terms do not re- 
ally lead to interesting phenomenology as in metric f{R) 
gravity: for instance they cannot drive inflation as, un- 
like in the scenario proposed by ( Starobinskvl . Il980[ ) in 
the metric formalism, here there are no extra dynam- 
ics and inflation cannot end gracefully ( Meng and Wangj . 
l2004bl : ISotiri"OTjl . l2006al) . As a matter of fact, it is more 
likely that positive powers of TZ will lead to no interesting 
cosmological phenomenology unless their coeffi cients are 
large enough to make the models non-viable ( Sotirioul . 
l2006all . 



B. Cosmological eras 

As stated in the Introduction, the recent fiurry of 
theoretical activity on f{R) models derives from the 
need to explain the present acceleration of the uni- 
verse discovered with supe r novae of type la (lAstier et al. 



2006'; ''Harris et 

Knop et al. , 2003; Perlmuttcr et al. 
19981 IT999I 120041: ISchmldTera^.l . 



'2004"; 'FiliDDcnko and Ricss', '199S; 

, 1998; Ricss et ql 
Il998i: iTonrv et al 



2003) ). We have seen in the previous section how f{R) 



gravity can achieve cosmic acceleration and an effective 
equation of state parameter We// ^ ^1; on the other 
hand, it was already known from i?^-inflationary scenar- 
ios of the early universe that this is possible, so we are 
actually witnessing a resurrection of this theoretical pos- 
sibility in models of the late universe — this parallels 
the use of scalar fields to drive early inflation or late- 
time acceleration in quintessence models. There are also 
attempts to unify early in flation and late time acceler- 
ation in modified g: r avity ( Bamba and Odintsovi l2008l 
iNoiiri and Odintsovi . l2007dll2008aibl ldldD. However, any 
model attempting to explain the cosmic speed-up at late 
times should not spoil the successes of the standard cos- 
mological model which requires a definite sequence of eras 
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to follow each other, including: 

1. early inflation 

2. a radiation era during which Big Bang Nucleosyn- 
thesis occurs; 

3. a matter era; 

4. the present accelerated epoch, and 

5. a future era. 



Big Bang Nucleosyn th esis i s well constrained 



see 



(IBrookfield et al. 
2005at lEvans et al 



20071: 



Clifton and Barrow 



, .Kneller and Steigman 

2004; iLambiase and Scarpettal . 120061 : iNakamura et al 



20061 ) for such constraints on f{R) models. The matter 
era must last long enough to allow the primordial 
density perturbations generated during inflation to grow 
and become the structures observed in thd universe 
today. The future era is usually found to be a de Sitter 
attractor solution, or to be truncated at a finite time by 
a Big Rip singularity. 

Furthermore, there must be smooth transitions be- 
tween consecutive eras, which may not happen in all 
f{R) models. In particular, the exit from the ra- 
diation era has been studied and claimed to origi- 
nate problems for many forms of f{R) in the met- 



2(n-Hl) , 



n > 



200e; 



2006) 



20081: 



ric formalism, inclu d ing / = R — n , 

jA mendola et a/.l. l2007allbl: 'Brookfield ef^aZI. 
ICapoz ziollo e t al\ . l2006cl : [Noiir i and Odi ntsovl 
[but not in the Pa latini formalism (jCarvalho et al\ . 
iFav et all l2007b( )]. However, the usual model f[R) = 
i? — /x'^/i? with "bad" beha viour was studied u sing singu- 
lar perturbation methods ( Evans et al\ . 120071) . definitely 
finding a matter era which is also sufficiently long. 

Moreover, one can always find choices of the func- 
tion f{R) with the correct cosmological dynamics in 
the following way: one can prescribe the desired 
form of the scale factor a{t) and integrate a dif- 
ferential equation f or f{R) that pr o duces the de- 



sired scale factor (ICapozziello et al. 



2005bl l2006c : 



de la Cruz-Dombriz and Dobadol 2006 : Faulkner et al. 
20071: iFav et o.li l2007 al'b': ' Hu and Sawickl 120073111^ : 
Multa,maki and Viliai ri2006a: Nojiri and Odintsovl . l200e , 
2007bllcl : ISong et all [2Q07D . In general, this "designer 



f{R) gravity" produces forms of the function /(i?) that 
are rather contrived. Moreover, the prescribed evolution 
of the scale factor a{t) does not determine uniquely the 
forn i of f{R) but, at bes t , only a class of f { R) mo d- 
els (iMultamaki and Vilial . boOGal: ISokolowskil l2007bll3 : 
IStarobinskvl l2007^ . Therefore, the observational data 
providing information on the history of a{t) are not suf- 
ficient to reconstruct /(i?): one needs additional infor- 
mation, which may come from cosmological density per- 
turbations. There remains a caveat on being careful to 
terminate the radiation era and allowing a matter era 
that is sufficiently long for scalar perturbations to grow. 



While sometimes it is possible to find exact solutions 
to the cosmological equations, the general behaviour of 
the solutions can only be assessed with a phase space 
analy sis, which constitutes a powerful t o ol in cosmol- 
ogy (|Colevl . 120031: IWainwright and Ellisl . Il997t ). In a 
spatially flat FLRW universe the dynamical variable is 
the Hubble parameter H, and a convenient choice of 
phase space variables in this case is {H,R). Then, for 
any form of the function f{R), the phase space is a 
two-dimensional curved manifold embedded in the three- 
dimensional space (^H, R, Rj with de Sitter spaces as 

fixed points ( de Souza and Faraoni [20071 ): the structure 
of the phase space is simplified with respect to that of 
general scalar-tensor cosmology ( Faraoni l2005b[ ). 

Studies of the phase space of f{R) cosmology 
(not limited to the spatially flat FLRW case) were 
common in the p re- 199 8 hterature on j?^ -in flation 
(lAmendola et al'. . 19921: Capozziello et Ml, Il993t 
iMuUer et all \l99(i IStarobinskvl . Il980l) . Thepresence or 



absence of chaos in metric 
([Barrow and Cotsakisl [l98i: 



:ravity was studied in 
Such studies with 



dynamical system methods have become widespread 
with the recent surge of interest in f{R) gravity to 
explain the present cosmic acceleration. Of course, de- 



choices of the function f(R) ([Abdelwahab et all 2008 


Amendola et al.. 2007ajb c': 'Amendola and Tsuiikawa 


2008 


; 


Carloni et all 2008a: Carloni and Dunsby 


2007 


Carloni et all 2005 


. 20071: Carroh et all [2005 


Chfton 


. 2006a. 20071: 


Clifton and Barrowl [2005b 


Easson 


, 2004; Fav et all 1 




2008: Leach et all '2006. 2007: [Li and 


Barrowl [20071 


Noiiri and Odintsovl. ,2004bi: iSami et al.. 


20051). 



C. Dynamics of cosmological perturbations 

Obtaining the correct dynamics of the background 
cosmological model is not sufficient for the theory to 
be viable: in fact, the FLRW metric can be obtained 
as a solution of the field equations of most gravitation 
theories, and it is practically impossible to discrimi- 
nate between f{R) gravity and dark energy theories 
(or between different f{R) models) by using only 
the unperturbed FLRW cosmological model, i.e., by 
using only probes that are sensitive to the expansion 
history of the universe. By contrast, the growth of 
cosmological perturbations is sensitive to the theory of 
gravity adopted and constitutes a possible avenue to 
discriminate between dark energy and modified gravity. 
Changing the theory of gravity affects the dynamics of 
cosmological perturbations and, among other things, 
the imprints that these leave in the cosmic microwave 
background (which currently provide the m ost sensitive 
cosmological probe) and in galaxy surveys ([Knox et al 
2006; Koivisto, 2006 b; [Kovama and Maartensl [2006 : 
Li and Barrowl [20071: iLi and Chu[. [2006[: [Sealfon et al 
200A [Shirata et al\ . [20051 [20071 Iskordis et al 
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20061: ISong et all 12007 



Tsuiikawa 



Zhang et al\ . 




■ .Stabenau and Jainl 
IWhite and Kochanekl 



20061 iLi et all l2007bl : lUddin et all 



2001 



to constrain f{R) 
background data 



This originated various efforts 
Jjravity with cosmic microwave 



( Amendola and Tsuiikawal 2008 : 



Applebv and Battvel l2007t ICarloni et al^ ] l2008b : 



Hu and Sawicki', '2007b; Li and Barrow', '2007"; 'Li et al. 
2007a; Pog osian and Silvcstri, 2008; St arobinsky. ,2007 
Tsuiilawalj200g: ffsuiikawa eit'all . boOSHWei and Zhang 



20081) 



Most of these works are restricted to specific choices 
of the function f{R), but a few general results have 
also been obtained. The growth and evolution of lo- 
cal scalar perturbations, which depends on the theory of 
gravity employed, was studied in metric f{R) gravity the- 
ories wh ich reproduce GR a, t high curvatures in various 
pape r s (ICarroU et al. [■ 120061: |de la Cruz-Dombriz et all 
120081 : ISong et all l2007^ by assuming a scale factor evo- 
lution typical of a ACDM model. Vector and tensor 
modes are unaffected by f{R) corrections. It is found 
that f"{R) > is requi r ed for the stability of scalar per- 
turbations (|Song et all [2003), which matches the anal- 
ysis of Sec. (jV.B.2p in a locally de Sitter background. 
The corrections to the Einstein-Hilbert action produce 
qualitative differences with respect to Einstein gravity: 
they lower the large angle anisotropy of the cosmic mi- 
crowave background and may help explaining the ob- 
served low quadrupole; and they produce different cor- 
relations betwe en the cosmic micr owave background and 
galaxy surveys ( Song et oZI . |2007( ). Further studies chal- 



lenge the viability of .f(R) gravity in comparison with the 
ACDM model: in ("Bea n et aH . l2007l) it is found that large 
scale density fluctuations are suppressed in comparison 
to small scales by an amount incompatible with the ob- 
servational data. This makes it impossible to fit simulta- 
neously large scale data from the cosmic microwave back- 
ground and small scale data from galaxy surveys. Also, 
a qua si-static approximation used in a previous analysis 
(| Zhan g. 2007) is found to be invalid. 

In (jde la Cruz-Dombriz et all . l2008l ). the growth of 
matter density perturbations is studied in the longitu- 
dinal gauge using a fourth order equation for the den- 
sity contrast Sp/ p, which reduces to a second order one 
for sub-horizon modes. The quasi-static approximation, 
which does not hold for general forms of the function 
f{R), is however found to be valid for those forms of this 
function that describe successfully the present cosmic ac- 
celeration and pass the Solar System tests in the weak- 
field limit. It is interesting that the relation between the 
gravitational potentials in the metric which arc responsi- 
ble for gravitational lensing, and the matter overdensities 
depends on the theory of gravity; a study of this relation 
in f{R) gravity (as well as in other g ravitational theories) 
is contained in ( Zhang et all [20071 ) . 

Cosmological density perturbations in the Palatini 



formalism have been studied in (lAmarzguioui et al. 


20061: Carroll et all 20061: Koivisto. 2006b 


- 


2007 


Koivisto and Kurki-Suonio, 120061 iLeel 2007 


2008 



l2007fl. Two differ- 

ent formalisms developed in ([Hwang; and Noh . [2002 
[Koivisto and Kurki-Suoniol 120061 ) and (|Lue e^'aLl . 1200 



were compared for the model f{R) — R — p 



2(n+l) I Jin 

and it was found that the two models agree for scenarios 
that are "close" (in parameter space) to the standard con- 
cordance model, but give different results for models that 
differ significantly from the ACDM model. Although this 
is not something to worry about in practice (all models 
aiming at explaining the observational data are "close" 
to the standard concordance model), it signals the need 
to test the validity of perturbation analyses for theories 
that do differ significantly from GR in some aspects. 



V. OTHER STANDARD VIABILITY CRITERIA 

In addition to having the correct cosmological dynam- 
ics and the correct evolution of cosmological perturba- 
tions, the following criteria must be satisfied in order for 
an f{R) model to be theoretically consistent and com- 
patible with experiment. The model must: 

• have the correct weak-field limit at both the New- 
tonian and post-Newtonian levels, i.e., one that is 
compatible with the available Solar System exper- 
iments; 

• be stable at the classical and semiclassical level (the 
checks performed include the study of a matter in- 
stability, of gravitational instabilities for de Sitter 
space, and of a semiclassical instability with respect 
to black hole nucleation); 

• not contain ghost fields; 

• admit a well-posed Cauchy problem; 

These independent requirements are discussed sepa- 
rately in the following. 



A. Weak-field limit 

It is obvious that a viable theory of gravity must have 
the correct Newtonian and post-Newtonian limits. In- 
deed, since the modified gravity theories of current inter- 
est are explicitly designed to fit the cosmological obser- 
vations. Solar System tests are more stringent than the 
cosmological ones and constitute a real testbed for these 
theories. 



1. The scalar degree of freedom 

It is clear from the equivalence between /(i?) and 
Brans-Dicke gravities discussed in Sec. IIIII that the for- 
mer contains a massive scalar field 4> [see eqs. ([51)1 and 
([M)) ]. While in the metric formalism this scalar is dy- 
namical and represents a genuine degree of freedom, it 
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is non-dynamical in the Palatini case. Let us, therefore, 
consider the role of the scalar field in the metric forma- 
lism as it will turn out to be crucial for the weak-field 
limit. Using the notations of Sec. IIII.AI the action is 
given by eq. ([54)) and the corresponding field equations 
by eq. ([55]). 

Equation (|52l) for x has no dynamical content because 
it only enforces the equality x — R- However, x = i? is 
indeed a dynamical field that satisfies the wave equation, 

3/"(x)nx + 3/"'(x)V„xV"x 

+ x/'(x)-2/(x)-«T. (91) 

When /" 7^ a new effective potential Vt^(x) 7^ ^(x) can 
be introduced, such that 



dW _ nT-xf'ix)+2fix) 
dx 3 fix) 



(92) 



The action can be seen as a Brans-Dicke action with 
Wo = if the field = /'(x) = f'{R) is used instead of x 
as the independent Brans-Dicke field: 



S = 



(93) 



where V {4>) is given by eq. (|53p . 

Now one may think of studying the dynamics and sta- 
bility of the model by looking at the shape and extrema 
of the effective potential V{x) but this would be mis- 
leading because the dynamics of x are not regulated by 
V{x) (indeed, the wave equation (|9T1) does not contain 
V), but are subject the strong constraint x — R, and R 
(or /'(i?)) is ruled by the trace equation 

The following example shows how the use of the poten- 
tial V{x) can be misleading. As is well known, the effec- 
tive mass of a scalar field (corresponding to the second 
derivative of the potential evaluated at the minimum) 
controls the range of the force mediated by this field. 
Thus, when studying the weak-field limit of the theory it 
is important to know the range of the dynamical scalar 
field X = present in the metric formalism in addition 
to the metric field 5^^, as this field can potentially vio- 
late the post-Newtonian constraints obtained from Solar 
System experiments if the scalar field gives observable ef- 
fects at the relevant scales. One way to avoid Solar Sys- 
tem constraints, however, is to have x have a sufficiently 
short range (see Sec. IV. A. 21 for more details). Consider 
the example f{R) — R + aR^, with a a positive constant. 
By naively taking the potential, one obtains 



V{x)^ax' 



2 ^ 



(94) 



with effective mass squared = 2a. Then, the small 
values of a generated by quantum corrections to GR im- 
ply a small mass mi and a long r ange field \ might 
be detectable at Solar Syste m scales ( Chiba et qLI . l2007t 
iJin et aJ. l. l2006HOlmol . l2007D . However, this conclusion is 
incorrect because toi is not the physical mass of x- The 



true effective mass is obtained from the trace equation ([8]) 
ruling the evolution of R which, for f{R) — R + aR^, re- 
duces to 



DR-^^^, 
6a 6a 



(95) 



and the identification of the mass squared of x = -R as 



m 



1 

6a 



(96) 



is straightforward.^^ A small enough value of a now leads 
to a large value of m and a short range for x- The situ- 
ation is, however, more complicated; the chameleon effect 
due to the dependence of the effective m ass on the curva- 
ture may change the ran ge of the scalar ( Faulkner et al\ . 
120071 : IStarobinskvl l200iD . 

For a general f{R) model, the effective mass squared of 
X = i? is obtained in the weak-field limit by considering 
a small, spherically symmetric, perturbation of de Sitter 
space with constant curvature Rq. One finds 



Ro 



(97) 



This equation coincid es with eq. (6) of (jMuller et all 
|1990D. with eg. f26) of (IOlmoLl2007D . and with eq. (17) of 
( Navarro and Van AcoleveiiT 20071 ). It also appears in a 
calculation of the propagato r for fjR) gravity in a locally 
fiat background (eq. (8) of (jNunez and Solganilj . 12004 )). 
The same expression is recovered in a gauge-invariant 
stabi lity analysis of de Sitter space ( Faraoni and Nadeaul . 
l2005t ) reported in Sec. EED below. 

Another possibility is to consider the field 4> = f'{R) 
instead of x = -R, and to define the effective mass of 
(j) by using the Einstein frame scalar-tensor analog of 
f{R) gravity instead of its Jordan frame cousin already 
discussed (Chibai, i2003i ). By performing the conformal 
transformation 



and the scalar field redefinition — f'{R) cj) with 



2a;o + 3 



2k 



(99) 



a scalar-tensor theory is mapped to the Einstein frame 
in which the "new" scalar field </) couples minimally to 



It was already noted by fStcUi, Il978t) that an correction to 
the Einstein-Hilbert Lagrangian generates a Yukawa correction 
to the Newtonian potential — this has to be kept small at macro- 
scopic scales by giving it a short range. 

The deflection of light by the Sun in GR plus quadratic cor- 
rections was studied by calculating the Feynman amplitudes for 
photon scattering, and it was fou nd that, to linearized order, this 
deflection is the same as in GR llAcciolv et aZ.1. 119991) . 
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the Ricci curvature and has canonical kinetic energy, as 
described by the action 



%4> - U{^) 



^g^.,^), (100) 



(note once more the non-minimal coupling of the matter 
in the Einstein frame). For the wq = equivalent of 
metric f{R) gravity we have 



<j, ^ /'(i?) = eV^^ 
Rf'iR) - f{R) 



2n{f{R)Y 



(101) 
(102) 



where R — R{4)). By using dcp/dcj) 
live mass of (/) is defined by 



^ the effec- 



eff 



d02 



1 



4/ 
/' iff 



(103) 



[this equation appears in the footnote on p. 2 of (jChibal . 
l2003t )]. By assuming a de Sitter background with con- 
stant curvature Rq — 12Hq = /o/(6/q), this turns into 



(104) 



In the Einstein frame, it is not the mass to of a particle or 
a field that is measurable, but rather the ratio ifi/ifiu be- 
tween m and the Einstein frame unit of mass rhu, which 
is varying, scaling as to„ = [f'{R)] rriu — rriu, 
where is the constant uni t of mass in the Jordan frame 
(iDickd . 119621: iFaraoni et ai\ . Il999l: iFaraoni and Nadeaul . 
I2007D . Therefore, 



,5,2 



TO 



eff 



m 



eff 



(105) 

In practice, (f> = f'{R) is dimensionless and its value must 
be of order unity in order to obtain the gravitational cou- 
pling strength measured in the Solar System; as a result, 
the Einstein frame metric g^i, and the Jordan frame met- 
ric 5^1/ are almost equal, and the same applies to to„, to„ 
and to rhef / , nief / , respectively. Then, the only rele- 
vant difference between Einstein and Jordan frames is 
the scalar field redefinition 6^6. 



2. Weak-field limit in the metric formalism 

Having discussed the field content of the theory, we 
are now ready to discuss the weak- field limit. Having 
the correct weak-field limit at the Newtonian and post- 
Newtonian levels is essential for theoretical viability. 

From the beginning, works on the weak-field 
(Newtonian and post-Newtonian) limit of f{R) 



gravity led to opposite re sults ap pearing in the 



litera tu re (lAcciolv et al 



2007 



fBarrow and Cli 
2006dL l2007Jd: ' 



"ton 




Baghram et al 



Capozziello et al 



Capozziello and Troisj ^ 



Capozziello and Tsuiikawa 



2006 
2004 
lorio 



Clifton and Barrom 
■Easso4 120041: 
20071: 



2005 



20071: ICembranos , 
2005al |2006t 



Dick . 

'Hu and Sawicki', '2007b; 
Multamaki a nd Vilia, 2006b, 2007 a; 

b; 



Navarro and Van Acolevenl. 1 2005. 2006; Olm g 12005 



Raiaraman. 2003; Ruggiero and lorio, 20 07|: j Shao et al. 
2a06; Soussa and Woodard, 2004; Zhang, l2007|). More 



over, a certain lack of rigour in checking the convergence 
of series used in the expansion around a de Sitter 
background often left doubts even on results tha t, a 
posteriori, turned out to be correct ( Sotirioul . l2006c( l . 



By using the equivalence between f{R) and scalar- 
tensor gravity, Chiba or iginally suggested that all f{R) 
theories are ruled out ( Chibal [2003). This claim was 
based on the fact that metric f{R) gravity is equivalent 
to an ojo = Brans-Dicke th eory, while the observ ational 
constraint is \ujq\ > 40000 (|Bertotti et all 120031 ). This 
is not quite the case and the weak-field limit is more 
subtle than it appears, as the discussion of the previ- 
ous section might have already revealed: The value of 
the parametrized post-Newtonian (PPN) parameter 7, 
on which the observational bounds are directly applica- 
ble, is practically independent of the mass of t he scalar 
only when the latter is small ( Wagoneil . Il970l ). In this 
case, the constraints on 7 can indeed be turned into con- 
straints on ujQ. However, if the mass of this scalar is 
large, it dominates over loq in the expression of 7 and 
drives its value to unity. The physical explanation of 
this fact, as mentioned previously, is that the scalar be- 
comes short-ranged and, therefore, has no effect at Solar 
System scales. Additionally, there is even the possibility 
that the effective mass of the scalar field itself is actually 
scale-dependent. In this case, the scalar may acquire a 
large effective mass at terrestrial and Solar System scales, 
shielding it from experiments performed there while be- 
ing effectively light at cosmological scales. This is the 
cha meleon mechanism, well-known in quintessence mod- 
els (|Khourv and Weltmam . i2004ai ib) . 



Given the above, it is worth examining these issues in 
more detail. Even though early doubts about the validity 
of the dynamical equi valence with scalar-tensor theory in 
the Ne wtonian limit ( Faraonil l2006bl: Kainulainen et all . 
l2007aD have now been dissipated (jFaraonil . l2007b[) . a di- 
rect approach which does not resort to the scalar-tensor 
equivalence i s preferable a. s the former could in principle 
hide things (|01mol . l2005bl ). This was gi ven in the met- 
ric fo rmalism, first in the special case ( Erickcek et all . 
l2006f ) fiR) = R- //i? [whic h is already ruled out 
by th e Ricci scalar i nstability ( Dolgov and Kawasakil 
l2003al: iFaraonj l2006al) ] and in the case f{R) = i?" us- 
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ing: light deflection and othe r Solar System experiments^^ 
([Barrow and Cliftonl. [2006t IClifton and Barrowl l2005al 
[2OO6; Za kharov et aZj . 120061 ). Only later was the case 
of a gene r al function f (R) d iscussed ( Chiba et a/.l . l2007t 
Ijin et cd. I. l2006l:l01moll2007f) . Chiba's result based on the 
scalar-tensor equiyalence eventually turns out to be valid 
subjec t to cert ain assumptions which are not always sat- 
isfied (jChiba et at, 2007; Jin et at, 2006; Olmo, 20Q3) 
— see below. This method, however, does not apply to 
the Palatini version of f{R) gravity. 

In what follows we adhere to, but streamline, the dis- 
cussion of (IChiba et aLl . l2007D with minor modifications, 
in order to compu te the PPN pa rameter 7 for metric f{R) 
gravity [see also ( Olmd . l2007l )]. We consider a spheri- 
cally symmetric, static, non-compact body embedded in 
a background de Sitter universe; the latter can exist in 
an adiabatic approximation in which the evolution of the 
universe is very slow in comparison with local dynamics. 
The condition for the existence of a de Sitter space with 
Rfii^ = ^o.9/ji//4 and constant curvature Rq = 12Hq is 



foRo — 2/0 — 0, 



Ho 




(106) 



The line element is 



ds^ = - [1 + S^lr) - H^r^] dt^ 

+ + 2$(r) -I- H^r^] dr^ + r^dfl^ 



(107) 



in Schwarzschild coordinates, where the post-Newtonian 
potentials 5'(r) and <i>(r) are treated as small pertur- 
bations.^^ The goal is to compute the PPN parameter 
7 — —4'/$ by solving the equations satisfied by these 
potentials. A linearized analysis is performed assuming 



and 



|4'(r)|,|$(r)| « 1, r«Ho 



R{r) = Ro + Ri{r), 



(108) 



(109) 



where the deviation Ri{r) of the Ricci curvature from 
the constant Rq is also a small perturbation.^" 
Three assumptions are made: 



The perihelion precess ion i n modified gravity is studied in 
llBaeh ram et al, 2007; lorid . |2007| : llorio and Ruegierd r2007allbt 
ISchmidii. 120081) . 

Isotropic coordinates are usually employed in the study of the 
weak-field limit of spherically symmetric metrics; however, the 
difference is irrelevant to first order in ^ and <I> jOlmo, 2007). 
The solution derived for the spherically symmetric metric is only 
valid when mr << 1, where m is the effective mass of the scalar. 
If t his assumption is not made, then (for example, according 
to l|jin et al\ . l200d )). it would seem that quantum corrections 
in f{R) = R + aR? with a ~ 10'^'^ GeV'^ are ruled out by 
Solar System constraints, which is not the case because these 
corrections are equivalent to a massive scalar field with short 
range that is not constrained by the available data. 



Assumption 1: f{R) is analytical at Rq. 
Assumption 2: mr << 1, where m is the effective mass 
of the scalar degree of freedom of the theory. In other 
words, this scalar field (the Ricci curvature, which is an 
extra dynamical quantity in the metric formalism) must 
have a range longer than the size of the So lar System— 
if it is much shorter than, say, 0.2 mm ( Hoyle et al\ . 
l200lh . the presence of this scalar is effectively hidden 
from Solar System and terrestrial experiments. In this 
case, this field could not have cosmological effects at late 
times, but could only be important in the very early uni- 
verse at_Jii£h_curvatures, e.g., in Starobinsky-like infla- 
tion (|Starobinskvlll980l) . 

Assumption 3: the pressure P ~ for the energy- 
momentum of the local star-like object. The trace of 
the corresponding energy-momentum tensor reduces to 
Ti ~ -p. 

By expanding f{R) and f'{R) around i?o, the trace 
equation (|153p reduces to 



3/^'ni?i + (/^'i?o-/^)i?i =«;Ti, 



(110) 



where T — Ti, since T is zero in the background. For 
a static, spherically symmetric body, Ri = Ri{r) and 
DRi = V^Ri ^ j^-^ {r'^^)- The reduced trace equa- 
tion (|110p then becomes 



where 



2 _ /o - fpRo 
By using i?o = 12i?o = '^fo/ fo, this reduces to 



m 



(111) 



(112) 



(113) 



This equation is found in vari ous other treatments of per- 
turbations of de Sitter spac e (iFaraon i and Nadeaul. 2005 ; 
Navarro and Van Acolevenl . 2007; Nunez and Solganik , 
20041 : lOlmol . 120071 ). 

Assumption 2 that the scalar i?i is light, which enables 
the f{R) theory to produce significant cosmological ef- 
fects at late times, also allows one to neglect^^ the term 
rn^Ri in cq. (jllip . The Green function of the equation 
V^i?i = —-S-Mt is then Gir) — ~-r^ and the solution is 

^ 3/q ^ '' 47rr 

i?i ~ / d^x' ^^y^G{r - r'), which yields 



Ri 



3fo 

kM 



[mr « 1) . 



(114) 



Although l lChiba et all l2007f) provide Green functions in both 
cases > and m? < 0, the latter corresponds to a spacetime 
instability and is unphysical. This is is irrelevant in the end 
because onl y the case — » is necessary and used in the 
calculation jFaraoni and Lanahan-Tremblavl . |2007| ). 



24 



Now, the condition m r << 1 yields 



1 

3 

and, using H^r << 1, 



/o p 
777 - rto 

JO 



<< 1 



« 1. 



(115) 



(116) 



Let us use now the fuh field equations ([6]); by expand- 
ing f{R) and f'{R) and using /o = QH^fo we get 



- /^'V"V^i?i + = k7>". (117) 

By using again the assumption Hqt << 1, the 
d'Alembertian □ becomes and, for {11,1^) = (0,0), 



(118) 

By computing Rq = 3Hq — V'^\l/(r) and dropping terms 
fo'H^Ri « /qV^^-, etc., we obtain 



f'^V'-^{r)+'-^R,- r^V'R^ = Kp. 



(119) 



Recalling that V^i?i 2:^ — for mr << 1, one obtains 

■J Jr 



/oV'*(r) 



2k p 



(120) 



Eq. (|120p can be integrated from r = to r > rp (where 
To is the radius of the star-like object) to obtain, using 
Gauss' law. 



kM 



dr 67r/^ 487r/^'r2 



(121) 



where M[r) — An dr' [r'Yp{r'). The integration con- 
stant Ci must be set to zero to guarantee regularity of 
the Newtonian potential at r = 0. The potential 5'(r) 
then becomes 



*(r) 



K.M 



kM 



487r/^' 



(122) 



The second term on the right hand side is negligible; in 
fact, 



K, Air 












-kM 







<< 1, 



and 



kM 
67r/o r ' 



(123) 



(124) 



Let us now find the second potential <I>(r) appearing in 
the line element (|107p . By using the field equations © 
with {a,b) = (1,1), 

+ f^RiR\ + /^'Di?! = ^Tl (125) 
with Tl ~ outside the star, and 



R\ — 3-ffg 



d^^f 2 (i$ 
dr^ r dr ' 



(126) 
(127) 



an d neglecting higher order terms, one obtains (eq. (22) 
of (jChiba et all 120071) 1 



dr 



0. (128) 



Now, using eq. ()114p for i?i , one concludes that the third 
term in eq. ()128p is negligible in comparison with the 
fourth term. In fact. 



/o-Ri 

2 






r dr 







« 1. 



(129) 



Then, using again the expression (|114p for dRi/dr and 
eq. (|124p for 4'(r), one obtains 



dr 



12^/^ r' 

which is immediately integrated to 
<I>(r) 



1271/^ r- 



(130) 



(131) 



The post-Newtonian metric (|107p therefore gives the 
PPN parameter 7 as 



7 



$(r) 
l'(r) 



(132) 



This is a gross v iolation of the expe rimental bound 
I7- 1| < 2.3-10-5 (|Bertotti et a^.l . 120031 ) and agrees with 



the calculation of the PPN parameter 7 = found 
by using the equivalence of metric /(f?) gravity with an 
Wo = Brans-Dick e theorv (IChibal. I2003D . 

The results of , (Chiba et all l2007t ) have been re- 
produced by (jOlmo. . ,2007) , who works in isotropic 

coordinates with a slightly different approach. 

( Kainulainen et al\ . l2007al ) have obtained spherically 
symmetric interior solutions matched to the exterior 
solutions of metric f{R) gravity and have confirmed the 
result 7 — 1/2. 

Limits of validity of the previous analysis: One can 
contemplate various circumstances in which the assump- 
tions above are not satisfied and the previous analysis 
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breaks down. It is important to ascertain whether these 
are physically relevant situations. There are three main 
cases to consider. 

Th e case of non-analytic f{R): While ()Chiba et al\ . 
l2007t ) consider functions /(i?) that are analytic at the 
background value i?o of the Ricci curvature, the situ- 
ation in which this fun ction is not analytical has been 
contemplated briefly in (jJin et a l.. 2006). Assuming that 
f{R) has an isolated singularity at i? = i?^, it can be 
expressed as the sum of a Laurent series. 



/(i?) = ^ a„ {R-RsY 



(133) 



n=0 



(|jin et al\ . l2006l ) note that it must he R ^ Rs in the 
dynamics of the universe because a constant curvature 
space with R = Rs can not be a solution of the field 
equations. Therefore, one can approximate the solution 
adiabatically with a de Sitter space with constant cur- 
vature Rq 7^ Rg. The function f{R) is analytical here 
and the previous discussion applies. This is not possi- 
ble if f{R) has an essential singularity, for example, if 

f{R) ^ R-n^ sin (|Jin et aLl . l2006l) . There is, of 

course, no reason other than Occam's razor to exclude 
this possibility. 

Short range scalar field: If the assumption mr << 1 
is not satisfied, the scalar is massive. If its range is suf- 
ficiently short, it is effectively hidden from experiments 
probing deviations from Newton's law and from other 
Newtonian and post-Newtonian experiments in the solar 
neighbourhood. This is the case of quadratic quantum 
corrections to Einstein's gravity, e.g., f{R) ~ R + aR^. 
If the effective mass is to > 10"'^ eV (corresponding to a 
fifth force range less than ~ 0.2 mm, the shortest scale 
currently accessible to weak- field experiments), this cor- 
rection is undetectable and yet it can sti ll have large 
effect s in the early inflationary universe (jStarobinskvl . 
I1980D . However, it can not work as a model for late 
time acceleration. 

Chameleon behavio ur: The chameleon effect 

( Khourv and Weltmanl . l2Q04al lbh. originally discov- 
ered in scalar field models of dark energy, consists of the 
effective mass to of the scalar degree of freedom being 
a function of the curvature (or, better, of the energy 
density of the local environment), so that to can be large 
at Solar System and terrestrial curvatures and densities, 
and small at cosmological curvatures and densities — 
effectively, it is short-ranged in the Solar System and 
it becomes long-ranged at cosmological densities thus 
causing the acceleration of the universe. Th e chameleon 
effect can be applied to met r ic f{R) gravity (ICembranos , 
2006; 'Faulkn er et 'ail. 120071: H avarro and Van Acolevenl . 
2007.: StarobinskvL I2007D. with the result that theories 



of th e kind (|Amendola et all l2007allbl l3: ICarroh et 'all . 
I2004D 



f{R) =R-{1- n)^^ ( 4 



(134) 



are compatible with the observations in the region of 
the p arameter space < n < 0.25 with /i sufficiently 
small ( Faulkner et all l2007t ) . Prec isely, using the Cassin i 
bound on the PPN parameter 7 ( Bertotti et all |2003() . 
the constraint 



n(l — n) 



2(l-r.) -6-5r. 



(135) 



is obtained ( Faulkner et all |2007[ ). Fifth force experi- 
ments give the bounds 



no 



n(l — n) 



10" 



(136) 



Preferred values seem to be to ~ 10 eV~ 10 Hq 
(|Faulkner et aLl . l2007l ). Mote that n > 0, which guaran- 
tees /" > 0, is required for Ricci scalar stability (n = 
reduces the model to GR with a cosmological constant, 
but avoiding the latter was exactly the reason why dark 
energy and modified gravity were introduced in the first 
place). 

These models work to explain the current cosmic ac- 
celeration because, for small curvatures i?, the correc- 
tion in i?" with n < 1 is larger than the Einstein- 
Hilbert term R and comes to dominate the dynam- 
ics. On the negative side, these theories are obser- 
vationally indistinguishable from a cosmological con- 
stant and they have beeii dub bed "vanilla f(R) gravity" 



Amendol a et al. 
2008: Faulkner et al\ 



2007alf9 l3: I Amendola and Tsuiikawal . 

izD- 



20071 ). However, they still have the 



advantage of avoiding a fine-tuning problem in A at the 
price of a much smaller fine-tuning of the parameter /i. 
As for all modified gravity and dark energy models, they 
do not address the cosmological constant problem. 

The weak-field limit of metric f{R) theories which ad- 
mit a global Minkows ki solution a round which to lin- 
earize, was studied by ( Cliftonl . l2008l ). These theories (in- 
cluding, e.g., analytic functions f{R) = X^n^ anR") are 
not motivated by late time cosmology and the Minkowski 
global solution, although p resent, may not be stable 
fclifton and Barrowl r2005al ). which in practice detracts 
from the usefulness of this analysis. Several new post- 
Newtonian potentia ls are found to appear in addition to 
the two usual ones (ICliftonl . [2008h . 



3. Weak-field limit in the Palatini formalism 

Early works on the weak-field limit of Pala- 
tini f{R) gravity often led to contradictory re- 
sults and to several technical problems as well 



(lAUemandi et all 2005bl; lAllemandi and Ruggiero 


2007; 


Barraco and Hamitv. 2000|: Bustelo and Barraco 


2007; 


Dominguez and Barracol 12004 ; Kainulainen et al. 


20071 


r. Meng and Wang. 2004a; Olmol. 2005a 


b|. 2007 


Ruggierol 120071: IRuggiero and loriol 12007: 


Sotiriou 



i2006a) which seem to have been clarified by now. 
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First of all, there seems to have been some confusion 
in the literature about the fact that Palatini f{R) grav- 
ity reduces to GR with a cosmological constant in vac- 
uum and the consequences that this can have on the 
weak-field limit and Solar System tests. It is, of course, 
true (see Sec. III.Bp that in vacuo Palatini f{R) grav- 
ity will have the same solutions of GR plus a cosmo- 
logical constant and, therefore, the Schwarzschild-(anti- 
)de Sitter solution will be the unique vacuum spherically 
symmetric solution (see also Sec. IVI.CJTTj for a discussion 
of t he Jebsen-Birkhoff theor e m). This was interpreted 
in (lAllemandi and Ruggieroi I20Q7I : iRuggiero and loriol 
l2Q07li " as an indication that the only parameter that can 
be constrained is the effective cosmological constant and, 
therefore, models that are cosmologically interesting (for 
which this parameter is very small) trivially satisfy Solar 
System tests. However, even if one sets aside the fact that 
a weak gravity regime is possible inside matter as well, 
such claims cannot be correct: they would completely 
defeat the purpose of performing a parametrized post- 
Newtonian expansion for any theory for which one can 
establish uniqueness of a spherically symmetric solution, 
as in this case we would be able to judge Solar System 
viability just by considering this vacuum solution (which 
would be much simpler). 

Indeed, the existence of a spherically symmetric vac- 
uum solution, irrespective of its uniqueness, does not suf- 
fice to guarantee a good Newtonian limit. For instance, 
the Schwarzschild-de Sitter solution has two free param- 
eters; one of them can be associated with the effective 
cosmological constant in a straightforward manner (using 
the asymptotics). However, it is not clear how the sec- 
ond parameter, which in GR is identified with the mass 
of the object in the Newtonian regime, is related to the 
internal structure of the object in Palatini f{R) gravity. 
The assumption that it represents the mass defined in 
the usual way is not, of course, sufficient. One would 
have to actually match the exterior solution to a solution 
describing the interior of the Sun within the realm of 
the theory in order to express the undetermined param- 
eter in the exterior solution in terms of known physical 
quantities, such as Newton's constant and the Newtonian 
mass. The essence of the derivation of the Newtonian 
limit of the theory consists also in deriving such an ex- 
plicit relation for this quantity and showing that it agrees 
with the Newtonian expression. The parametrized post- 
Newtonian expansion is nothing but an alternative way 
to do that without having to solve the full field equa- 
tions. Therefore, it is clear that more information than 
the form of the vacuum solution is needed in order to 
check whether the theory can satisfy the Solar System 
constraints. 

However, some early attempts towards a Newto- 
nia n and post-Newtonian e xpans ion were also flawed . 
In (|Meng and Wand , l2004a^ and ( Barraco and Hamitvl . 
|200(F for instance, a series expansion around a de Sitter 
background was performed in order to derive the Newto- 



nian limit. Writing 



n = 7^o + 7^l, 



(137) 



where TZq is the Ricci curvature of the background and 
72.1 is a correction, one is tempted to expand in powers 
of TZi/TZq regarding the latter as a small quantity. Since 
one needs the quantities f{TZo + TZi) and f'{TZo + TZi), 
the usual approach is to Taylor-expand around TZ — TZq 
and keep only the l eading order term s in TZi. However, it 
has been shown in ( Sotirioul . l2006cl ) that this can not be 
done for most cosmologically interesting models because 
TZi/TZo is not smaU. 

Take as an example the model ( Vollickl |2003() 



(138) 



and £2 - lO-^^(eV) 



Expanding as 



fin) - fiUo) + f'ino)ni + -f"ino)nj + ... (139) 



and using eq. ()138p yields 



fin) = fiUo) 



nl 2111 ' 



(140) 



where now TZo — ^2- It is then easy to see that the second 
term on the right hand side is of the order of TZ\ , whereas 
the third term is of the order of TZ\/e2- Therefore, in 
order to truncate before the third term, one needs TZi ^ 
nil €2 or 



£2 >72i. 



(141) 

£2 and so this is 



This is not a stringent constraint: 7?.o 
the usual condition for linearization. 

Let us return now to the trace equation ([20|) . For the 
model under consideration, 



7^ = i ( -kT± 



k2 T2 



12£2 



(142) 



According to eq. (|142p . the value of TZ, and conse- 
quently, TZi, is algebraically related to T and, whether 
or not the condition (|141|) is satisfied or not critically 
depends on the value of the energy density. To demon- 
strate this, pick the mean density of the Solar System, 
p ^ 10~^"'^gr/cm'^, which satisfies the weak-field limit cri- 
teria. For this value, \e2/i^T\ ~ lO'^i^ where T ^ ~p. 
The "physical" branch of the solution (|142p is the one 
with positive sign because, given that T < 0, it ensures 
that matter leads to a standard positive curvature in 
strong gravity. Then, 



7^. 



-kT 



3e| 

kT 



(143) 



and 7^l kT r-. up. Thus, £2/7^l lO^^i ^nd it 

is evident that the required condition does not hold for 
some typical densities related to the Newtonian limit. 
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The situation does not improve even with the "unphys- 
ical" branch of eq. (|142p with a negative sign. In fact, in 
this case, TZi ^ e2[3e2/(KT) + \/3] and the correction to 
the background curvature is of the order 62 and not much 
smaUer than that, as it would be required in order to 
trunc ate the expansion (|140p . In ( Barraco and Hamitvl . 
I2OOOD . this fact was overlooked and only linear terms in 
TZi were kept in the expansion of f(TZ) and f'{TZ) around 
72.0. In ( Mens and Wand . l2004aj) . even though this fact 
is noticed in the final stages of the analysis and is actu- 
ally used, the authors do not take it into account properly 
from the ou tset, keeping again only l inear terms [see, e.g., 
eq. (11) of (|Meng and Wantd . l2004al) ]. 

However, the algebraic dependence of 72 on the density 
does not only signal a problem for the approaches just 
mentioned. It actually implies that the outcome of the 
post-Newton ian expansion it self depends on t he density, 
as shown in (|01mol . l2005al lbl: ISotiriou', '2006c'). Consider 



for instance, along the lines of (jSotiriou . .2006a, . c). the 
conformal metric 

h^. = f'in)g,,, (144) 

that was introduced in Sec. III.BI [cf. eq. (|22l) ]. In terms 
of this metric, the field equations can be written in the 
form 



n^. - \nh^, + (/' - 1) 



7? -JLh 
^'^ 2 /' 



= kT, 



(145) 

and TZfiu is the Ricci tensor of the metric h^,^. It is ev- 
ident that, if /' — 1, then h^^ and coincide and 
eq. p45p yields Einstein's equation. However, since 72. 
and consequently /'(72) are functions of the energy den- 
sity, due to eq. (pOj) . the deviation of /' from unity will 
always depend on the energy density and the functional 
form of /. Therefore, one can definitely find some func- 
tion / which, for some range of energy densities, will give 
/' = 1 to high precision. However, for the same function 
/, there will be large deviations from /' = 1 at a dif- 
ferent density range. This dependence of the weak-field 
limit on the energy density is a novel characteristic of 
Palatini f{R) gravity. 

This dependence can be made explicit if the prob- 
lem is app r oached via the equivalent Brans-Dicke the- 
ory ( Olmol [2005al lbh. Note that the usual bounds com- 
ing from Solar System experiments do not apply in the 
cjQ ~ —3/2 case, which is equivalent to Palatini /(i?) 
gravity. This is because the standard treatment of the 
post-Newtonian expansion of Brans-Dicke theory, which 
one uses to arrive to such bounds, is critically based on 
the assumption that wq 7^ —3/2 and the term (2a;o -f 3) 
frequently appears as a denominator. Making this as- 
sumption is not necessary, of course, in order to derive 
a post-Newtonian expansion, but is a convenient choice, 
which allows for this otherwise general treatment. There- 
fore, a differe nt approach, such as the one followed in 
(IOlmd . [2005bl) . was indeed required for the ujq = —3/2 
case. Following the standard assumptions of a post- 
Newtonian expansion around a background specified by a 



cosmolog ical solution (|Wilill98l[ ). the following relations 
were derived for the post-Newtonian limit: 



-v^[hl^-n{T)\ - 



2<j) 
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(146) 
(147) 



where V is the potential of the scalar field 4> and i^{T) = 
log[(/)/0o]- The subscript in (/)q, and in any other quan- 
tity in the rest of this subsection, denotes that it is eval- 
uated at T = 0. 

The solutions of eqs. (|146p and (|147p are 



,(1) 



(t,a-): 



2GeffA/e 



27GcffMe 



Vr 



2 
r 



(148) 
(149) 



where Mq = (f>o J dr'x' p (t, x') /(f). The effective Newton 
constant Goft and the post-Newtonian parameter 7 are 
defined as 



G, 



off 



G/ M. 

cjiO \ Mq 
Mq - My 

Mq+Mv' 



(150) 
(151) 



where My = k^Vo / d^x' [Vo/(j)o - V((l))/ ^]. 

As stated in different words in ( Olmol . l2005b[ l. if the 
Newtonian mass is defined as Mat = J (Px'p{t,x'), the 
requirement that a theory has a good Newtonian limit 
is that GeffM0 equals GM^, where N denotes Newto- 
nian, and 7 ~ 1 to very high precision. Additionally, the 
second term on the right hand side of both eqs. (|148p 
and (|149p should be negligible, since it plays the role of 
a cosmological constant term. ft{T) should also be small 
and have a negligible dependence on T. 

Even though it is not impossible, as mentioned before, 
to prescribe / such that all of the abov e are satisfied fo r 
some range of densities within matter fSotirioul. 120063) . 
this does not seem possible over the wide range of den- 
sities relevant for the Solar System tests. As a matter 
of fact, is nothing but an algebraic function of T and, 
therefore, of the density (since (p is an algebraic func- 
tion of 72.). The presence of the ft{T) term in eqs. (|148p 
and (|149p signals an algebraic dependence of the post- 
Newtonian metric on the density. This direct dependence 
of the metric on the matter field is not only surprising but 
also seriously problematic. Besides the fact that it is ev- 
ident that the theory cannot have the proper Newtonian 
limit for all densities (the range of densities for which it 
will fail depends on the functional form of /), consider 
the following: What happens to the post-Newtonian met- 
ric if a very weak point source (approximated by a delta 
function) is taken into account as a perturbation? And 
will the post-Newtonian metric be continuous when going 
from the interior of a source to the exterior, as it should? 
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We will refrain from further analysis of these issues 
here, since evidence coming from considerations differ- 
ent than the post-Newtonian limit, which we will review 
shortly, will be of significant help. We will, therefore, 
return to this discussion in Sec. IVI.C.21 



B. Stability issues 

In principle, several kinds of instabilities need to be 
considered to make sure that /(i^ ) grav i ty is a viable 
alternative to G R (ICalcagni et all l2006t IChibal. 12005 
De FeUce et d\ . l200fil: ISokolowskiL l2007al lhllcl: IWan 



2005D 



The Dolgov-Kawasaki ( Dolgov and Kawasakil l2003al ) 
instability in the matter sector, specific to metric f{R) 
gravity, imposes restrictions on the functional form of 
/, and is discussed below. More generally, it is be- 
lieved that a stable ground state, the existence of which 
is necessary in a gravitational theory, should be highly 
symmetric, such as the de Sitter, or Minkowski, or 
perhaps the Einstein static space. Instabilities of de 

Sitter space in the gravity sector have been found in 
iT^ — — — n'fcTTTr^ rrr-: ^ ^ - ■ i r^TT;^. 



6 



Barrow and Hervik, 2006a; 'Dolgov and Pelliccia', '2006'; 



sec 



Faraoni^ 2004b , c , 2005a ; Faraoni and Nadeau 2005 ^ 
also (jBarrow and Ottewill . Il983l: iMuller et al\ . iiOQOf) for 
pre-1998 discussions], while stability in first loop quan- 
tization of f{R) gravity and w ith respect to black 



hole nucleation was studi e d in ( Cognola et 



20081: ICognola and Zerbinj 120061: 



Paul and Pau 



2005 



2006 . 



20051 ). The linear stability of de Sitter space with re- 



spect to homogeneous perturbations in generalized theo- 
ries of the form f (R, R^^Ri^" R^^^ gR'"'"''^) was stud- 
ied in ( Cognola and Zerbini l2008f) . The stability of 
the Einstein static space in metric f{R) gravity with 
respect to homogeneou s perturbations was studied in 
( Boehmer et aLl . [2007bl ). while stability of this space with 
respect to inhomogeneous isotropic perturbations was es- 
tablished, with a gau ge-invariant formalism and under 
certain conditions, in ( Goswami et all\200^ . 



1. Ricci stability in the metric formalism 

In the metric formalism, Dolgov and Kawasaki dis- 
covered an instability in the prototype model f{R) = 
R — fi'^/R (now called "Dolgov-Kawasaki", or "Ricci 
scalar" or "matter" instability) , which manifests itself on 
an extremely shor t time scale and i s sufficient to rule out 
this model (.Dolgov and Kawasakil [2003a). Their re sult 
was confirmed in ( Noiiri and Odintsovl[2003al[2004bl ). in 
which it was also shown that adding to this specific f{R) 
an R^ term re moves this instab i lity. T he instability was 
rediscovered in ([Baghram et all |2007[ ) for a specific form 
of the function f{R). The analysis of this instability is 
generalized to arbitrary f(R) theories in the metric for- 
malism in the following way ([Faraoni [2006a} ). 

We parametrize the deviations from Einstein gravity 



fiR)=R + e^iR), 



(152) 



where e is a small parameter with the dimensions of a 
mass squared and (p is arranged to be dimensionless (in 
the example f = R — /i^/i?, one has e = /i^, = —fjL^/R, 
and fi-Ho^ 10-33 eV). 

By using the trace equation ^ 

inf\K) + f{R)R-2f{K)^KT, (153) 

and evaluating □/', 



(ey^-l) „^ 2^ 

(154) 

We assume that (p" ^ 0: if tp" = on an interval then the 
theory reduces to GR. Isolated zeros of ip" , at which the 
theory is "instantaneously GR" , are in principle possible 
but will not be considered here. 

Consider a small region of spacetime in the weak-field 
regime and approximate locally the metric and the cur- 
vature by 

gf^i^ ^ Vp-i^ + hp.^, R = -kT + Ri, (155) 

where 77^^ is the Minkowski metric and \Ri/kT\ <^ 1. 
This inequality excludes the case of conformally invariant 
matter with T = 0, a situation considered later. Equa- 
tion (|155p yields, to first order in 

, 2k Lp'" ■ ■ 2k tp'" ^ 
Ri - W^Ri ^ TRi + — ^ VT • VRi 



ip" 



+3;7U"''J — ' 

(156) 

where V and V'^ are the gradient and Laplacian in 
Euclidean three-dimensional space, respectively, and an 
overdot denotes differentiation with respect to time. The 
function ip and its derivatives are now evaluated at i? = 
— kT. The coefRcient of Ri in the fifth term on the left 
hand side is the square of an effective mass and is domi- 
nated by the term (3e ip") ^ due to the extremely small 
value of e needed for these theories to reproduce the cor- 
rect cosmological dynamics. Then, the scalar mode Ri of 
the f{R) theory is stable if ip" = f" > 0, and unstable if 
this effective mass is negative, i.e., if ip" — f" < 0. The 
time scale for this instability to manifest is estimated to 



be of the order of the inverse effective mass 
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in the example e(p{R) = —fi'^/R ([Dolgov and Kawasakil 
l2003a[) . The small value of ip" gives a large effective mass 
and is responsible for the small time scale over which the 
instability develops. 

Let us consider now matter with vanishing trace T of 
the stress-energy tensor. In this case eq. (jl56p becomes 



Ri 



ip" 



Ri -V^Ri 



ip" 



(vs,)= 



iip" 



Zip" 



(157) 
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Again, the effective mass term is ~ {3eip")~ , which has 
the sign of /" and the previous stability criterion is re- 
covered. The stability condition f"{R) > is useful to 
veto f{R) gravity models. 

When /" < 0, the instability of these theories can be 
interpreted, following eq. (|156p . as an instability in the 
gravity sector. Equivalently, since it appears inside mat- 
ter when R starts deviating from T [see eq. (jl55p ]. it can 
be seen as a matter instability [this is th e interpretation 
taken in ( Dolgov and Kawasakil l2003al) ]. Whether the 
instability arises in the gravity or matter sector seems to 
be a matter of interpretation. 

The instability of stars made of any type of matter 
in theories with /" < and sufficiently small is con- 
firmed, with a differ ent approach (a generalized varia- 
tional principle) in ( Seifertj . |2007[ ). in which the time 



scale for instability found by Dolgov and Kawasaki in 
the 1/R model is also recovered. The stability condition 
/" > is recovered in studies of cosmological perturba- 
tions (jSawicki and Hul . l2007l ). 

The stability condition f"{R) > 0, expressing the 
fact that the scalar degree of freedom is not a ghost , 
can be given a simple physical interpretation ( Faraoni 
l2007bl ). Assume that the effective gravitational coupling 
Geff{R) = G/ f'{R) is positive; then, if Ge// increases 
with the curvature, i.e., 



dG. 



dR 



-f"{R)G 



>o. 



(158) 



at large curvature the effect of gravity becomes stronger, 
and since R itself generates larger and larger curvature 
via eq. (|153|) . the effect of which becomes stronger and 
stronger because of an increased Gef f{R), a positive feed- 
back mechanism acts to destabilize the theory. There is 
no stable ground state if a small curvature grows and 
grows without limit and the system runs away. If in- 
stead the effective gravitational coupling decreases when 
R increases, which is achieved when f"{R) > 0, a neg- 
ative feedback mechanism operates which compensates 
for the increase in R and there is no running away of 
the solutions. These considerations have to be inverted 
if /' < 0, which can only happen if the effective energy 
density pe/f also becomes negative. This is not a phys- 
ically meaningful situation because the effective gravi- 
tational coupling becomes negative and the tensor field 
and the scalar field of metri c f{R) gravity become ghosts 
(jNunez and Solganiki [200I . 

GR, with f"{R) = and Gefj = constant, is the bor- 
derline case between the two behaviours corresponding to 
stability (/" > 0) and instability (/" < 0), respectively. 

Remarkably, besides the Dolgov-Kawasaki instability 
which manifest itself in the linearized version of equa- 



tion (|153p . there are also recent claims that R can be 
driven to infinity due to strong non-linear effects related 
to the same equation ( Applebv and Battvd . 12008 : Frolov . 
'2008; 'Tsuiikawal, l2008l ). More specifically, in (|Tsuiilffl!wa , 
i2008. ) an oscillating mode is found as a solution to the 
perturbed version of eq. (|153p . This mode appears to 
dominate over the matter-induced mode as one goes back 
into the past and, therefore , it can violate the stability 
conditions. In (|Frolovll2008l) . eq. was studied, with 

the use of a convenient variable redefinition but without 
resorting to any perturbative approach. It was found that 
there exists a singularity at a finite field value and energy 
level. The strongly non-linear character of the equation 
allows R to easily reac h the singularity in the pres ence of 
matter. As noticed in ("Appl ebv and Battvd . 120081) . since 
when it comes to cosmology the singularity lies in the 
past, it can in principle be avoided by choosing appropri- 
ate initial conditions and evolving forward in time. This, 
of course, might result in a hidden fine-tuning issue. 

All three studies mentioned consider models in which 
/(i?) includes, besides the linear term, only terms which 
become important at low curvatures. It is the form of 
the effective potential governing the motion of R, which 
depends on the functional form of f{R), that de termines 
how easy it is to drive R to infinity ( Frolovl . |2008| ) . There- 
fore, it seems interesting to study how the presence of 
terms which become important at large curvatures, such 
as positive powers of R, could affect these results. Fi- 
nally, it would be interesting to see in detail how these 
findings manifest themselves in the case of compact ob- 
jects, and whether there is any relation between this issue 
and the Dolgov-Kawasaki instability. 



2. Gauge-invariant stability of de Sitter space in the metric 
formalism 

One can consider the generalized gravity action 
/(</., i?) u;{cf,) 



S 



d^x^ 



v"0v„0-y(0) 



(159) 

incorporating both scalar-tensor gravity (if / (0, R) — 
4>{<p)R) and modified gravity (if the scalar field 4> is ab- 
sent and /rr 7^ 0). In a spatially flat FLRW universe 
the vacuum field equations assume the form 



3/' 



Rf f 



^^-^ (^r + F-Hr 



3H(P 



1 f duj 



2uj 



df 



V - 3Hf' ) ,(160) 
(161) 



=0,(162) 



22 Refs. l lMultamaki and Viliall2006al : lNoiirill2004l'l hinted towards 
the stabihty criterion, but did not fully derive it because a de- 
composition in orders of was not performed. 



where /' = df/d(j), F = df /dR, and an overdot denotes 
differentiation with respect to t. We choose {H,(j)) as 
dynamical variables; then, the stationary points of the 
dynamical system (|160p - (|162l) are de Sitter spaces with 
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constant scalar field {Ho,4>o). The conditions for these 
de Sitter solutions to exist are 



6H^ /o - /o + 2Fo = 0, 
df ^dV 



0, 



(163) 
(164) 



where = /'(^o,-Ro), /o = /(0o,-Ro), Vo = V (M, 
and Rq — 12Hq. The phase space is a curved two- 
dimensional surface embedded in a three-dimensional 



space (|de Souza and Far aoni. 2007). 

Inhomogeneous perturbations of de Sitter space have 
been studied using the covari a nt and gauge-in v ariant 
formalism of (iBardeenl . Il980l: lEllis and Brml. Il989 : 
Ellis et a/ J, '1990', '1989") in a version pr ovided by (|Hwane . 
1990allbl . 11997. ,1998; .Hwang and Nohl . Il996l ) for general- 
ized gravity. The metric perturbations are defined by 

500 = -a' (1 + 2Ar) , go, = -a^BY,, (165) 
g,, = a" [h,j (1 + 2HlY) + 2Ht Y,j] . (166) 



Here the Y are scalar spherical harmonics, hij is the 
three-dimensional metric of the FLRW background, Vi is 
the covariant derivative of hij , and k is the eigenvalue of 
V,;V*F ~ —k^Y. The Yi and Yij are vector and tensor 
harmonics satisfying 



1 ^ 



1 



1 



y. = -- v,r, r,, = -2 V, v,y + -y h,, , (i67) 



fc2 

respectively. The Bardeen gauge-invariant potentials 



B — — Ht 
k 



a / a ■ 
^a = A+-[B--Ht 



a 



the EUis-Bruni variable 



B 



(168) 



(169) 



(170) 



and analogous gauge-invariant variables A/, and 
Ai? sat isfy f irst o rder equations given in ([Hwand . 
Il990al lbl. Il997l . Il998t iHwang and Nohi r il996D . which sim- 
plify sign ificantl y in the de Sitter backgrou i id (g p, c/jq) 
(|Faraonil l2004bl l2005at iFaraoni and Nadeaul . |2005[ ). 

To first order and in the absence of ordinary ma t- 
ter, vector perturbations do not ap pear (lHwand . fl990allbl . 
ll997lll998tlHwang and Nohl . ll996D . and de Sitter space is 
always stable with respect to first order tensor perturba- 
tions. Focusing on scalar perturbations, modified grav- 
ity corresponds to </) = 1 and / = f{R) with f"{R) ^ 
in (I159p . The gauge- invariant perturbations (from 
which one easily obtains and Ai?) satisfy 



H 



3Ho^H + ( - 4i?o' + ^ ) = 0, (171) 



(|Faraonil l2004bl l2005at iFaraoni and Nadeaul . l2005l) . 
where the term k'^ jc? can be dropped at late times and 
for long wavelength modes. Linear stability ensues if the 
coefficient of is non-negative, i.e. (using eq. (|163p ). 



/o/o 



> 0. 



(173) 



The only term containing the comoving wave vector k 
in eq. (|171l) becomes negligible at late times and/or for 
zero-momentum modes and thus the spatial dependence 
effectively disappears. In fact, eq. (|173|) coincides with 
the stability condition that can be obtained by a straight- 
forward homogeneous perturbation analysis of eqs. (|160p 
and (|161l) . As a result, in the stability analysis of de 
Sitter space in modified gravity, inhomogeneous pertur- 
bations can be ignored and the study can be restricted to 
the simpler homogeneous perturbations, which are free of 
the notorious gauge-dependence problems. This result, 
which could not be reached a 'priori but relies on the 
inhomogeneous perturbation analysis, holds only for de 
Sitter spaces and not for different attractor (e.g., power- 
law) solutions that may be present in the phase space. 
The stability condition (|173p is equivalent to the condi- 
tion that the scalar field potential in the Einstein frame 
of the equivalent Brans-Dicke theory has a minimum at 
the configu ration ident i fied by the de Sitter space of cur- 
vature Rq (jSokolowskill20Q7b[ ). 

As an example, let us consider the prototype model 



fiR) = R 



R 



(174) 



The background de Sitter space has Rq — 12H'q — fi^ 
and the stability condition (|173p is never satisfied: this 
de Sitter solution is always unstable. An improvement is 
obtained by adding a quadratic correction to this model: 



f{R) = R- ^+aR\ 



(175) 



Then, the condition for the existence of a de Sitter so- 
lution is again Rq = v^/i^, while the stability condi- 



tion (|173p is satisfied if a > , in agreement with 

(|Noiiri and Odintsovl l2003ai l2004bl) who use indepen- 
dent methods. 

Different definitions of stability lead to different, al- 
beit close, stability criteria for de Sitter space [see 
( Cognola et aLl . [20051 2008) for the sen iiclassical stability 
of modified gravity, ( Bertolainl Il987l ) for scalar-tensor 
gravity, and (|Seifertl . 120071 ^ for a variational approach ap- 
plicable to various alternative gravities]. 



The generalizatio n of the condi tion II173II to D spacetime dimen- 
sions, derived in l|Radoij . l2007l ) for homogeneous perturbations, 



IS 



(D-2) {f^y-Dfof- 



> 0. 



(172) 
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3. Ricci stability in the Palatini formalism 



C. The Cauchy problem 



For Palatini f{R) gravity the field equations 
and ([T9|) are of second order and the trace equation ((20|) 
is 



/'(7^)7^-2/(7^) = Kr, 



(176) 



where TZ is the Ricci scalar of the non-metric connection 
r{^^ (and not that of the metric connection {{^^} of g^jy). 
Contrary to the metric case, eq. (I176P is not an evolution 
equation for TZ; it is not even a differential equation, but 
rather an algebraic equation in TZ once the function / (TZ) 
is specified. This is also the case in GR, in which the 
Einstein field equations are of second order and taking 
their trace yields R = —kT. Accordingly, the scalar field 
(j) of the equivalent ojq = —3/2 Brans-Dicke theory is not 
dynamical. Therefore, the Dolgov-Kawasaki instability 
can not occur in Palatini /(i?) gravity (|Sotirioui . i2007ail . 



4. Ghost fields 

Ghosts (massive states of negative norm that cause 
apparent lack of unitarity) appear easily in higher order 
gravities. A viable theory should be ghost-free: the 
presence of ghosts in f(R, R,,u R^" . R,,„„_aR'"'"^) 
gravi ty has been studied in (iBuchbinder et al. 



19921: ICodello and Percaccil 120061: iDe Fehce 
De Felice and Hindmarshl . 120071: IStelld. Il976 . 



StromingeiT 



Vilkoviskv 



119841: 
I1992D . 



Utivama and DeWittl . 



1977 



1962t 



Due to the Gauss-Bonnet iden- 
tity, if the initial action is linear in Rfj.ua^R'^'^"^ , 
one can reduce the theory under consideration to^^ 
f{R,Rf^,^R'^'^) which, in general, contains a massive 
spin 2 ghost field in addition to the usual massless gravi- 
ton and the massiv e scala r . f(R) theories haveno ghosts 



(IBuchbinder et all Il992 : 



1978L Il977l: IStromingeT 



Ferraris et al 
1984 



1988; Stelle 



Utivama and DeWitt 



1962i IVilkoviskvl 119921). and the stability condi- 



tion f"(R) ^ of (jPolgov and Kawasakil l2003al: 
iFaraonil l2006a ) essentially amounts to guarantee that 
the scalaron is not a ghost. Theories of the kind 
f{R,Rf,^R^"',Rf,^p^R^"'P'') in general are plagued by 
ghosts [this is the case, for example, of conformal 
gravity, as noticed long before t h e 199 8 discovery of 
the cosmic acceleration ( Riegertj . Il984l )]. but models 
with only / {R, R^ - AR^^RP" + R^upaR^"'"') terms in 
the action have been claimed to be gh ost- free (|Comellil . 
120051: iNavarro and Van Acolevenll2006l) . 



Furthermore, Ruv W^'^ can be ex pressed in terms of in a 
FLRW background ll Wandj . [l993) . 



A physical theory must have predictive power and, to 
this extent, a well-posed initial value problem is a re- 
quired feature. GR satisfies this requirem ent for most 
reasonable forms of matter ( Waldl . Il984[) . The well- 



posedness of the Cauchy problem for f{R) gravity is an 
open issue. Using harmonic coordinates, Noakes showed 
that theories with action 



S 



1 

2k 



aR^,RP" + PR )+Sm (177) 



in the metric form alism ha v e a w ell posed initial value 
problem in vacuo (iNoakesL I1983D . By using the dy- 
namical equivalence with the scalar-tensor theory (j54p 
when f"{R) ^ 0, the well-posedness of the Cauchy 
problem can be reduced to the analogous problem for 
Brans-Dicke gravity with cjq = (metric formalism) or 
LjQ — —3/2 (Palatini formalism). The fact that the initial 
value problem is well-posed wa s demonstrated for par- 
ticular scalar-tensor theories in (ICocke and Cohenl . ll968l: 
INoakesL IT983h and a genera l analysis has recentl y been 
presented in ( Salgadol 120061 : ISalgado et all l2008l ). This 
work, however, does not cover the wq = 0, —3/2 cases. 

A system of 3 -I- 1 equations of motion is said to be 
well-formulated if it can be rewritten as a system of equa- 
tions which are of only first order in both time and space 
derivatives. When this set can be put in the full first 
order form 



(178) 



2007t where u collectively denotes the fundamental variables 
hij,Kij, etc. introduced below, M* is called the charac- 
teristic matrix of the system, and S (u) describes source 
terms and contains only the fundamental variables but 
not their derivatives. The initial value formulation is 
well-posed if the system of partial differential equations 
is symmetric hyperbolic [i.e., the matrices are sym- 
metric) and strongly hyperbolic if SiM^ has a real set of 
eigenvalues and a complete set of eigenvectors for any 
1-form Sj, a nd obeys some boundedness conditions [see 
(|Solinl . l2006l) ]. 

The Cauchy problem for metric f{R) gravity is well- 
formulated and is well-posed in vacuo and with matter, 
as shown below. For Palatini f{R) gravity, instead, the 
Cauchy problem is unlikely to be well-formulated or well- 
posed unless the trace of the matter energy-momentum 
tensor is constant, due to the presence of higher deriva- 
tives of the matter fields in the field equations and to the 
impossibility of eliminating them (see below). 

A systematic covariant approach to scalar-tensor the- 
ories of the form 



2k 2 



d"<l)do,(l> - W{(l)) 



S 



M 



(179) 

is due to ( Salgadol l2006l ). who showed that the Cauchy 
problem of these theories is well-posed in the absence of 



32 



matter and well-formulated otherwise. With the excep- 
tion of ujq — —3/2, as we will see later, most of Salgado's 
results can be extended to the more general action 



2k 2 



d''(f>do,(p - W{<j>) 



-S 



M, 



(180) 

which contains the additional coupling function lo{4>) 

i which is diff erent from the Brans-Dic ke parameter ujq) 
Lanahan-Tremblav and Faraoni l20Q7t ). 

The field equations, after setting k = 1 for this section, 
are 



Guy 



1 



(181) 



a;n</-+yi?-VK'(0) + yV"(/)Vc.0 = O, (182) 

where a prime denotes differentiation with respect to (j). 
Eq. (jl8ip can be cast in the form of the effective Einstein 



equatio n Ggi, = Tgu 



tensor (jSakadol [2006h 

V('/>) V '^'^ '"^ J' 



with the effective stress-energy 

(183) 



and 



-f^'(0)(V^V,</>- 5^,00), (184) 
;(0) (^V^0V,0 - ig^,V'3(/.V^0^ 

-M^(0)5m- (185) 
The trace of the effective Einstein equations yields 



T. 



3W0 
2Vj 



{uj + 3^") 



2V''1^(0) 



(186) 



The 3 + 1 Arnowitt-Deser-Misner (ADM) formulation of 
the theory proceeds by introducing lapse, shift, extrin- 
sic curvature, and gradients of (R cula . .1998 ; Salgado, 
120061 : IWaldl . Il984) . Assume that a time function t ex- 
ists such that the spacetime {M,g^i,) admits a foliation 
with hypersurfaces Ej of constant t with unit timelike 
normal n". The 3-metric and projection operator on Et 
are h^j^i, — g^^i, + n^Ui, and /i"^, respectively. Moreover, 



rt^n^ = — 1, 



hapn^ — hai3 , 



(187) 



The metric is then 

ds^ = - (Ar2 „ AT'AT.) ~ 2N,dtdx^ + h^jdx'dx^ (188) 

(i,j = 1,2,3), where iV > 0, = -A^V^i, A^" = 
-~h°' pt^ is the shift vector, while obeys f^V at — 1 and 
t° = -AT" TVn" so that A^ = -n^t" and N°'na = 0. 
The extrinsic curvature of is 



(189) 



and the 3D covariant derivative of ft-a^ on is defined 
by 

^T"-^,... = . . . h'^p, . . . h\VfT^^-5,... 

(190) 

for any 3-tensor ^^^T'^i - , with D^h^,^ = 0. The spa- 
tial gradient of the scalar is = D^(f) (where 
denotes the covariant derivative of h^,j), while its mo- 
mentum is n = Cn(j) — n'^^u4> Slid 

K,, = -V,7i, = --^ I ^ + D,N, + D.Ni] , (191) 



2N\ dt 

n^j^idt^ + N^Q,), (192) 
dtQi + N'diQi + QAN' = A (A^n) . (193) 

The effective stress-energy tensor T^^^^'^^ is decomposed 
as 



a, 

where 



1 



^a/ = ^ i^o^P + -^""/s + Jp'^o. + En^np) , (194) 



S^, ^ K-^h.'T^l^f' = i (^S^ + ^S^) ,(195) 



- -f£;(« -f . (197) 



Its trace is T^^fi") ^S-E, where 8 = 8''^. The Gauss- 
Codazzi equations then yiel d the Hamiltonian constraint 
(ISakadol [2OO6I: IWaldl . fl98l 



(198) 



(199) 



^^^R + K^ - KijK'^ = 2E, 
the vector constraint 

DiK\~D,K^J,, 
and the dynamical equations 

dtK'j + N^dilCj + K\d.jN^ - K^jdiN' + D'DjN 
-^^'>R'jN - NKK'j = ^[{S- E) 5] - 2S'j] , (200) 

where K = K\. The trace of eq. (^00]) yields 

dtK + N^diK+^^'^AN-NKijK'^ = ^ {8 + E) , (201) 
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where *^^^A = D'^Di. Our purpose is to eventually 
eliminate all second derivatives. which is present 

in eqs. (jl95p - ()197p can actually be eliminated using 
eq. (fTSS)) . provided that w ^ -3(?/'')V(2V')- 

To be more precise, a direct calculation yields the ip- 
and 0-quantities of eqs. (|195|) - (|197p 

£;W = ^j' {D^'Q,, + KIi)+i,"Q'', (202) 

jw = -^'(x2g^ + i?an)-v"nQ„, (203) 



where (ISakadol l2006h 



/:„n - IVK - Q^'D^ (IniV) ~ D^Q^" = -D(j) 

1 [ V''t(™) 



3(V>')^ 



2-0 



-^-(o. + SV- )y 



(214) 



[/ia/3 (Q' - tf) - Q,Q^] , (204) 
where = Q^Qv, while 



S'(^) = {D^Q"" +KI{- 300) + V" (3n2 - 2Q2) , 

(205) 

and 



i;W) = |(n2 + Q2)+p^(^)^ (206) 

j^"^) = -c^ng^, (207) 

5^5 = ^QaQ/j - Kp (Q' - tf) + Vl^(0)] ,(208) 

(209) 

The Hamiltonian and the vector constraints become 



7/; 



DiK\ - D,K + \ [i^' (K.'Q, + An) 



In vacuo, the initial data {hij,Kij,(f>,Qi,Il) on an ini- 
tial hypersurface Sq obey (|210p . (|21ip . and Qi = Di(j), 
DiQj = DjQi. In the presence of matter, the variables 

E("^\ j!r\ S^^^ must also be assigned on the initial hy- 
persurface Sq. Fixing a gauge corresponds to specifying 
the lapse and the shift vector. The system (I210p - (|213p 
contains only first-order derivatives in both space and 
time once the d'Alembertian is written in terms of 
(/), V^(/)Vp0, ip, and its derivatives by means of eqs. (|186p 
or pi4p . As mentioned earlier, this can be done whenever 
w# — 3(^')^/(2?/'). As already pointed out in ( Salgadd . 
|2006| ) for the specific case w = 1, and can be now gen- 
eralized for any w 7^ — 3(i/'')^/(2'0), the reduction to a 
first-order system shows that the Cauchy problem is well- 
posed in vacuo and well-formulated in the presence of 
reasonable matter. 

Let us now consider the results specific to f{R) gravity. 
Recall that Brans-Dicke theory, which is of interest for 
us due to its equivalence with /(i?) gravity, corresponds 
(e'^"'^ + W{(f>)\ , (210) to w(0) = a;o/0, with ^(0) = 0, and W 2V. This 

yields the constraints 



+ (w -f V") ng.] = 



J 



(m) 



(211) 



respectively, and the dynamical equation pOOp is 

+^ [r iQ' - tf) + 21^(0) + ^'□(/.j 



y((/))l , (215) 



Am) 



-\ [kIQi An + ^ng, ) = (216) 



and the dynamical equations 



24 



(212) 



dtK'j + N^diK'j + K\djN^ - K/diN' + D'DjN 

N 
24)' 



-5] {2V{^) + 



with trace 

dtK + N^diK +(3) AiV - NK,.K'^ 



N 



Nloq 



{D^'Q^ + UK) 



+^[^V-(2c. + 30n2] 



_ N 
~ 24 

dtK + N^diK +^^"1 AN - NK^K''^ 

-— (i^'^g. + ni^)- 



(217) 



N 
2^ 



(-2H/(0) - 3V''n</> + 5'"^ + £;'™^) (213) 



iV 
24> 



-2F(0) - 3n</> -f -f 



(218) 
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with 



UJQ 



-2V{cj))+cl)V'{cl)) + '^ 



(219) 

The condition oj ^ —3{ip')'^/{2tp), which needs to be sat- 
isfied in order to for one to be able to use eq. (|186p in 
order to eliminate Dcj) can be written in the Brans-Dicke 
theory notation as ujq ^ —3/2. One could of course had 
guessed that by looking at eq. (|219|) . Therefore, met- 
ric f{R) gravity, which is equivalent to wq = Brans- 
Dicke gravity, has a well-formulated Cauchy problem in 
general and is well-p osed in vacuo. Further work by 
( Salgado et all |2008^ established the well-posedness of 
the Cauchy problem for scalar-tensor gravity with uj = 1 
in the presence of matter; this can be translated into the 
well-posedness of metric f{R) gravity with matter along 
the lines established above. 

How about Palatini /(i?) gravity, which, correspond- 
ing to = — 3/2, is exactly the case that the constraint 
ujQ ^ —3/2 excludes? Actually, for this value of the 
Brans-Dicke parameter, eq. ([M)) . and consequently at 
eq. (|219l) . include no derivatives of cj). Therefore, one can 
actually solve algebraically for (p. [The same could be 
done using eq. (jl86p in the more general case where uj is 
a function of 4>, when uj = —3{ip')'^/{2ip).] We will not 
consider cases for which eq. (I69p has no roots or when it 
is identically satisfied in vacuo. These cases lead to theo- 
ries for which, in the Palatini /(i?) formulation, eq. (pij) 
has no roots or when it is identically satisfied in vacuo 
respectively. As already mentioned in Sec. III.Bl the first 
case leads to inconsistent field equ ations, and the second 
to a con formally i nvarian t theory ()Ferraris et all Il992f) , 
see also (ISotirioul . [2006bh for a discussion. 

Now, in vacuo one can easily show that the solutions 
of eq. (|69l) or (I219p will be of the form cj) — constant. 
Therefore, all derivatives of (fi vanish and one conclude 
that ujQ — —3/2 Brans-Dicke theory or Palatini f{R) 
gravity have a well-formulated and well-posed Cauchy 
problem. This could h ave been expected, as noticed in 
( Olmo and Singhl . |2009[ ). considering that Palatini f{R) 
gravity reduces to GR with a cosmological constant in 
vacuo. 

In the presence of matter, things are more complicated. 
The solutions of eq. ([S^ or (I219p will give as a function 
of T, the trace of the stress-energy tensor. This can still 
be used to replace cj) in all equations but it will lead to 
terms such as DT. Therefore, for the Cauchy problem 
to be well-formulated in the presence of matter, one does 



This has been missed in l lLanahan-Tremblav and Faraonill2007t) . 
where it is claimed that the Cauchy problem is not well-posed 
because the constraint loq —3/2 does not allow for the use 
of eq. I I219I I in order to eliminate 0<p. Note also that in the 
absence of a potential (there is no corresponding Palatini f{R) 
gravity) uiq = —3/2 Brans-Dicke theor y does not ha ve a well- 
posed Cauchy problem, as noticed in llNoakesl . Il983l ). because 
this theory is conformally invariant and ip is indeterminate. 



not only have to assume that the matter is "reasonable" , 
in the sense that the matter fields satisfy a quasilinear, 
diag onal, second order hyperbolic system of equations 
[see ( Waldl . Il984[) ]. but also to require that the matter 
field equations are such that allow to express all deriva- 
tives of T present in eq. ([STS]) to ([2T8| for = -3/2 
in terms of only first derivatives of the matter fields. It 
seems highly implausible that this requirement can be 
fulfilled for generic matter fields. This seems to imply 
that Wo = —3/2 Brans-Dicke theory and Palatini f{R) 
gravity are unlikely to have a well-formulated Cauchy 
problem in the presence of matter fields. However, more 
precise conclusions can only be drawn if specific matter 
fields are considered on a case by case basis. The com- 
plication arising from the appearance of derivatives of T, 
and consequently higher derivatives of the matter fields 
in the equations, and seem to be critical on whether the 
Cauchy problem can be well-formulated in the presence 
of matter will be better understood in Sec. I VI. C. 21 



VI. CONFRONTATION WITH PARTICLE PHYSICS AND 
ASTROPHYSICS 

A. Metric f{R) gravity as dark matter 

Although most recent motivation for f{R) gravity 
originates from the need to find alternatives to the myste- 
rious dark energy at cosmological scales, several authors 
adopt the same perspective at galactic and cluster scales, 

■ as a substitute for dark matter 



(CaDOzziello et al. 2004 


2005a 


2006al 2007al 


2006b 


120074 lorio and Ruggiero, 


2007a 


b: Jhingan et al. 


. 2008 


'Martins and Salucci], 


2007 


;jNoiiri and Odintsovl, 


2008a 



Saffari and Sobouti, 120071 : [Zliao and Li, 2003). Given 
the equivalence between f{R) and scalar-tensor gravity, 
these efforts resemble previo us attempts to m odel dark 
matter using scalar fields ( Alcubierre et all 2002a^ ; 
Bernal and Guzmanl. J2006a,b.c; Bernal and Matosl. 12005 
i2007a i h;~ iMatos and Guzman , 
2OOII : iMatos and Urena-Lopez 



Cervantes-Cota et 



2000; 



Matos et al. 



2OO1I. I 2OO7': Rodrigucz-Meza and Cervantes-Cotal 12004 : 
Rodrigucz-Mcza et al . 2005. , .20071^ 

Most works concentrate on models of the form 
f{R) = i?". A t heory of this form with n = 1 — 
a/2 was studied in ( Mendoza and Rosas- Gtievaral [2007t 



ISaffari and Soboutil . 2007() bv using sphericallv svmmet- 
ric solutions to approximate galaxies. The fit to galaxy 
samples yields 



a = (3.07 ±0.18) X lO"'^ 



M 



0.494 



(220) 



where M is the mass appearing in the spherically sym- 
metric metric (the mass of the galaxy) . Notice that hav- 
ing a depending on the mass of each individual galaxy 
straightforwardly implies that one cannot fit the data for 
all galactic masses with the same choice of f{R). This 
make the whole approach highly implausible. 
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( Capozziello et all l2004 l2005al l2006aL l2007a§ ) com- 
puted weak-field limit corrections to the Newtonian 
galactic potential and the resulting rotation curves; 
when matc hed to galaxy samples, a best fit yields 
n ~ 1.7. ( Martins and Salucci |2007() performed a 
fit u sing two broad e r samp les, finding n ~ 2.2 [see 
also ( Boehmer et all l2007al ) for a variation of this ap- 



proach focusing on the constant velocity tails of the 
rotation curves]. All these values of the parameter n 
are in vio lent contrast wit h the bounds obta i ned by 
(iBarrow and Cli fton. 200(| IClifton and Barrowl l2005al 
12006 : IZakharo\^e^ all [2000) and have been shown to vio- 
late also the current constraints on t he precession of peri- 
helia o f several Solar System planets ( lorio and Ruggierol 
l2007al lb[). In addition, the consideration of vacuum met- 
rics used in these works in order to model the gravita- 
tional field of galaxies is highly questionable. 

The potential obtained in the weak- field limit of f{R) 
gravity can affect other aspects of galactic dynamics as 
well: the scattering probability of an intruder star and 
the relaxation time of a stellar sy stem were studied by 
( Hadiimichef and Kokubunl . Il997 ). originally motivated 
by quadratic corrections to the Einstein-Hilbert action. 



B. Palatini f{R) gravity and the conflict with the Standard 
Model 

One very important and unexpected shortcoming of 
Palatini f{R) gravity is that it appears to be in con- 
fiict with the Standard Model of particle physics, in the 
sense that it introduces non-perturbative corrections to 
the matter action (or the field equations) and strong cou- 
plings between gravity and matter in the local frame and 
at low energies. The reason why we call this shortcom- 
ing unexpected is that, judging by the form of the ac- 
tion (I13p . Palatini f{R) gravity is, as we mentioned, a 
metric theory of gravity in the sense that matter is only 
coupled minimally to the metric. Therefore, the stress 
energy tensor is divergence- free with respect t o the metric 
covariant derivative, the metric postulates ( Willi . Il98l[ ) 
are fulfilled, the theory apparently satisfies the Einstein 
Equivalence principle, and the matter action should triv- 
ially reduce to that of Special Relativity locally. 

Let us see how this conflict comes about. This issue 

2,) using Dirac 



was first pointed out in (|Flanagan. . 20041 
particles for the matt er action as a n exam ple, and later 
on studied again in ( Iglesias et all [2003) by assuming 
that the niatter action is that of the Higgs field [see 
also (lOlmd . I2008D ]. Both calculations use the equiva- 
lent Brans-Dicke theory and are performed in the Ein- 
stein frame. Al though the us e of the Einstein frame has 
been criticized (IVollickl . l200l .^6. this frame is equivalent 
to the Jordan frame and both are perfectly suitable for 



performing calculations ( Flanaganl 2004a) [see also the 
relevant discussion in Sec. lIIII and ( Faraoni and Nadeaul . 
120071 : ISotiriou et aLl . [2007h ]. 

Nevertheless, since test particles are supposed to fol- 
low geodesies of the Jordan frame metric, it is this met- 
ric which becomes approximately fiat in the laboratory 
reference frame. Therefore, when the calculations are 
performed in the Einstein frame they are less transpar- 
ent since the actual effects could be confused with frame 
effects and vice-versa. Consequently, for simplicity and 
clarity, we prese nt the calculation in th e Jordan frame, 
as it appears in ( Barausse et aU[2008b[ ). We begin from 
the action (pS)) . which is the Jordan frame equivalent of 
Palatini f{R) gravity, and we take matter to be repre- 
sented by a scalar field H {e.g., the Higgs boson), the 
action of which reads 



Sm — — 



d.Hd^H^^H'] (221) 



(in units in which G — c — 1). As an example, we choos e 
f{R) = 7^ - fi'^/TZ (jCarroU et all \2004 IVoUickl . I2003D . 
For this choice of /, the potential is V{(j)) = 2/i^(0 — 
iy/2_ To go to the local frame, we want to expand the 
action to second order around vacuum. The vacuum of 
the action with (|22ip as a matter action is H = 0, 
(/) = 4/3 [using eq. (|69l) ] and gf^i, ~ rj^i, {fi^ ^ A acts as 
an effective cosmological constant, so its contribution in 
the local frame can be safely neglected). 

However, when one tries to use a perturbative ex- 
pansion for (f), things stop being straightforward: 4> 
is algebraically related to the matter fields as is obvi- 
ous from eq. (|69)) . Therefore, one gets Scj) ^ T / ^ 
rn^SH^ / {h? at energies lower than the Higgs mass 
(mn ^ 100 ~ 1000 GeV). Replacing this expression in 
the action (pS)) perturbed to second order, one immedi- 
ately obtains that the effective action for the Higgs scalar 



5" 



;ffccti^ 




(222) 



at energies E ^ mn- Taking into account the fact that 



K ^ Hi where H^^ 



4000 Mpc is the Hubble 
radius and 5H ^ mn because E <C toh, it is not difficult 
to estimate the order of magnitude of the corrections: 
at an energy E = 10"'^ eV (corresponding to the length 



° Note that in the case of l|Flanaganl . [2004bh in which fermions are 
used as the matter fields, one could decide to couple the indepen- 



dent connection to them by allowing it to enter the matter action 
and define the covariant derivative (which would be equivalent 
to assuming that the spin connecti on is as independent variable 
in a tetra d formalism ) , as noted in l lVoUickl l2005f) . Although the 
results of l|Flanaganl .f2004b') would cease to hold in this case, this 
can not be considered a problem: clearly in this case we would 
be talki ng about a different theory, namely metric-affine f{R) 
gravity l lSotiriou and Liberatil [20071J) . 
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scale L = h/E = 2 x 10 ^ m), the first correction is 
of the order m^5H^/{n^h^) - {H'^ /Xii)^{mu/Mp)^ > 
1, where Ah = ^/toh 2 x lO^^^ - 2 x 10"^* m 
is the Compton wavelength of the Higgs and Mp — 
{hc^/Gy/^ = 1.2 x 10^9 GeV is the Planck mass. The 
second correction is of the order m\{d5HY / [ji^'h?) ^ 
{H^^/XxH)HmH/Mp)^H^^/L)^ > 1. Clearly, having 
such non-perturbative corrections to the local frame mat- 
ter action is unacceptable. 

An alternative way to see the same problem would be 
to replace (50 ~ rn^SH^ / {hP ^j?) directly in (|66l) . Then 
the coupling of matter to gravity is described by the in- 
teraction Lagrangian 



d^Sg 



ft3 



L 



(223) 



This clearly exhibits the fact that gravity becomes non- 
perturbative at microscopic scales. 

It is obvious that the algebraic dependence of (p on 
the matter fields stands at the root of this problem. We 
have still not given any explanation for the "paradox" 
of seeing such a behavior in a theory which apparently 
satisfies the metric postulates both in the f{R) and the 
Brans-Dicke representation. However, this will become 
clear in Sec. IVLC.21 



C. Exact solutions and relevant constraints 

1. Vacuum and non-vacuum exact solutions 

Let us now turn our attention to exact solutions start- 
ing from metric f{R) gravity. We have already mentioned 
in Sec. lII.Al that. as it can be seen easily from the form of 
the field equations © , the maximally symmetric vacuum 
solution will be either Minkowski spacetime, if i? = is a 
root of eq. or de Sitter and anti-de Sitter spacetime, 
depending on the sign of the root of the same equation. 
Things are slightly more complicated for vacuum solu- 
tions with less symmetry: by using eq. ([6]) it is easy to 
verify that any vacuum solution (i?^!/ — h.g^^, T^^ — 0) 
of Einstein's theory with a (possibly vanishing) cosmolog- 
ical constant, including black hole solutions, is a solution 
of metric /(i?) gravity (except for pathological cases for 
which eq. © has no roots). However, the converse is not 
true. 

For example, when spherical symmetry is imposed, the 
Schwarzschild metric is a solution of metric f{R) grav- 
ity if i? = in vacuum. If R is constant in vacuo, 
then Schwarzschild-(anti-)de Sitter spacetime is a solu- 
tion. As we have already mentioned though, th e Jebsen- 
Birkhoff theorem (IWaldl . 11984 IWeinberd . Il972h does not 
hold in metric f{R) gravity [unless, of course one wishes 
to impose further c onditio ns, such as that R is constant 
( Capozziello et al\ . l2007d( )]. Therefore, other solutions 



can exists as well. An interesting finding is that the cos- 
mic no- hair theorem valid in GR and in pure f{K) gravity 
is not valid, in general, in theories of the form 



S ■ 



I3R,..R^'' - 2A) 



for which exact anisotropic 
continue to inflate anis otropic 
found (iBarrow and Hervip . 



(|Kluske and Schmidt! . Il996 



2006al 




iMaedal 


1988h 



ever, isotropizatio n during inflatio n occurs 
/ ((/), R) models (jMaeda et all Il989[ ) 

In addition to the exact cosmological 
explored for the purpose of explaii iing 



(224) 

solutions that 
y have been 
[see also 
] . How- 
mixed 



solutions 
the cur- 



rent c osmic acceleration [see^ e.g.. (lAbdalla et~ai. 



2005 



2007 



iBarrow and Cliftonl. 
IClifton and Barrowl 



120061: [Ciliftonl. l2006a. 



(ICapozziello and De Felice 



I2005U 
20M 



20061) and 



Capozziello et all 



l2007cl: iModak am 200'5rivakili 120081) for an approach 
based on Noether symmetries; see (jCarloni et dl I2006D 
for bouncing solutions and the conditions that they 
satisfy], exact spherically symmetric solutions of metric 
f{R) gravity have been explored in the literature, with 
most recent studies being motivated by the need to un- 
derstand the weak-field limit of cosmologically-motivated 
theories. 

Regarding non-vacuum solutions, the most common 
matter source is a perfect fluid. Fl uid dynamics in metric 
fjR) gravity was stud ie d in (iMaartens and Tavloi 



1994 iMohseni Sadiadil. l2007t iRippl et all Il996 : 



Tavlor and Maartensl . ll995l V Sp herically symmetric 



solutions were found in ([ Bronnik ov and Chernakova , 
2005a!'t?.'3: Bustelo and Barraco, 2007; Capozziello et ai, 
2007d: Clifton, ,2006a.b,: .Migncmi and Wihshire , 1992 ; 



Multamaki and Vilial . I2006U 



IWhittl . [19851) 

We regret not being able to present these solutions 
extensively here due to space limitations and refer the 
reader to the literature for more details. 

Stability issues for spherica lly symmetric solutions 
were discussed in (jSeifertl . (20071 ). In the theory 



[R - aR^ - 13R^^R'"'] 



(225) 



where a, /3, and e are constants and x is the Gauss- 
Bonnet invariant, the Schwarzschild metric is a solution 
and th e stabi l ity of Schwarzschild black holes was stud- 
ied in ( Whittl . Il985l) . Surprisingly, it was found that the 
massive ghost graviton ( "poltergeist" ) present in this the- 
ory stabilizes small niass black holes against quantu m in- 
stabilities [see also (jMvers and Simonl . [l988i Il989l) ]. In 
the case /? = e = 0, which reduces the theory (|225|) to 
a quadrat i c f{R ) gravity, the stability criterion found 
m (IWhittl . Il985l) reduces to a < 0, which corresponds 
again to f"{R) > 0. For a = we recover GR, in 
which black holes are stable classically (but not quantum- 
mechanically, due to Hawking radiation and their nega- 
tive specific heat, a feature that persists in f{R) gravity). 
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so the classical stability condition for Schwarzschild black 
holes can be enunciated as f"{R) > 0. 

Let us now turn our attention to Palatini f{R) grav- 
ity. In this case things are simpler in vacuo: as we saw 
in Sec. III.Bl the theory reduces in this case (or more 
precisely even for matter fields with T =const., where 
T is the trace of the stress energy tensor) to GR with 
a cosmological co nstant, which might as well be zero 
for some models (iBarraco et al\ . Il999l : iBorowiec et al\ . 
119981 : iFerraris et all ll992L 11991 . Therefore, it is quite 
straightforward that Palatini f{R) gravity will have 
the same vacuum solutions as GR with a cos mologi- 
cal constant. A l so, th e Jebsen-BirkhofF theorem (|Waldl . 
119841: IWeinberd. I1972D is valid in the Palatini formalism 
( Barausse et alV I2008al lbll3: iKainulainen et ai l. l2007bD . 

Cosmol ogical solutions in quadrat ic gravity were ob- 
tained in (jShahid-Salessl . 119901 Il99l[ ). Spherically sym- 
metric interior solutions in the Palatini formalism can 
be found by using the generalization of the Tolman- 
Oppenheimer-Volkoff equa tion valid for these theo - 
ries, which was found in ( Barraco and Hamitv*, '200 
iBustelo and Barrac 

E l2007t IKainulainen et al, 2007 
Indeed, such solutions have been found and matched 
with the unique ex terior (ant i-)de Sitter solution 
(iBarraco and Hamitv . 1998, 200 ^ Bustelo and Barracol 
120071: IKainulainen etaU . i2007allbf l Nevertheless, a 
matching between exterior and interior which can lead 
to a sensible solution throughout spacetime is not always 
feasible and this seems to ha ve serious conse quences for 
the viability of f{R) gravity ( Barausse et aZ.I . i2008a.b.c). 
This is discussed extensively in the next section. 

Let us close this section with some remarks on black 
hole solutions. As is clear from the above discussion, all 
black hole solutions of GR (with a cosmological constant) 
will also be solutions of f{R) i n both the metric and 
the Palatini forma l ism [s ee also (jBarausse and Sotirioul . 
120081: iPsaltis et aLl . 1200 8")]. However, in the Palatini for- 
malism they will constitute the complete set of black 
hole solutions of the theory, whereas in the metric for- 
malism other black hole solutions can exist in princi- 
ple, as the Jebsen-BirkhofF theorem does not hold. For 
a discussion on black hole e ntropy in f(R ) gravity see 
(IJacobson et a/.l . Il994 fT995t IVollickll2007D .^^ 



2. Surface singularities and tlie incompleteness of Palatini f{R) 
gravity 



In Secs. lV.A.3(IV.CI and lVLB| we aheady spotted three 
serious shortcomings of Palatini /(i?) gravity, namely 
the algebraic dependence of the post-Newtonian metric 
on the density, the complications with the initial value 
problem in the presence of matter, and a conflict with 



particle physics. In this section we will study static 
spherically symmetric interior solutions and their match- 
ing to the unique exterior with the same symmetries, 
the Schwarzschi l d-de S itter solution, along the lines of 
( Barausse et all . l2008al lblR). As we will see, the three 
problems mentioned earlier are actually very much re- 
lated and stem from a very specific characteristic of Pala- 
tini /(i?) gravity, which the discussion of this section will 
help us pin down. 

A common way of arriving to a full description of a 
spacetime which includes matter is to solve separately 
the field equations inside and outside the sources, and 
then match the interior and exterior solutions using ap- 
propriate jun ction conditi ons [called Israel junction con- 
dition in GR ( Israeli [i 966 *)]. This is what we are going to 
attempt here. We already know the exterior solution so, 
for the moment, let us focus on the interior. Since we as- 
sume that the metric is static and spherically symmetric, 
we can write it in the form 



(226) 



We can then replace this metric in the field equations 
of Palatini f{R) gravity, preferably in eq. (P5| which 
is the simplest of all the possible reformulations. As- 
suming also a perfect fluid description for matter with 
Tfj.1, = {p + P) Ufj^Uu + Pg^v, where p is the energy den- 
sity, P is the pressure, and u'^ is the fluid 4-velocity, 
and representing d/ dr with a prime, one arrives at the 
equations 
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where F = df/dR. To determine an interior 
solution we need a generalization of the Tolman- 
Oppenheimer-Volkoff (TOV) hydrostatic equilibrium 
equati on. This has been derived for Palatini f (R) grav- 
ity in (iBarraco and Haniitvl. l2000l : iBustelo and Barracol 
l2007l IKainulainen et al\ . l2007bj r Defining mtot(r) = 
r(l — e~^)/2 and using Euler's equation 



P' 



A' 



(P + P): 



(231) 



See also l lEling et all , |2006| ) for a derivation of the field equations 
of metric f{R) gravity based on thermodynamical arguments ap- 
plied to local Rindler horizons. 



In this section we modify our standard notation for economy 
and a prime denotes differentiation with respect to the radial 
coordinate instead of differentiation with respect to the argument 
of the function. 
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one can use eqs. (I227P and (I228|) to arrive to the gener- 
alised TOV equations: 



P' = - 



1 



ip + P) 



1 + 7 r(r - 2mtot) 



(232) 
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We now have three differential equations, namely 
(|23ip . (|232p and (|233p . and four unknown functions, 
namely >1, mtot (or B), P, and p. The missing piece 
is the information about the microphysics of the mat- 
ter configuration under investigation. In the case of a 
perfect fluid, this is effectively given by an equation of 
state (EOS). A one-parameter EOS relates the pressure 
directly to the energy density, i.e., P — P{p). A simple 
form of such an EOS is a polytropic equation of state 
P = kpQ, where po is the rest-mass density and k and 
r are constants. This is the case that we will consider 
here. Note that the rest-mass density can be expressed 
in terms of the energy density p and the internal energy 
U as Po = p — U . Assuming an adiabatic transformation 
and using the first law of thermodynamics, one can ex- 
press the internal energy in terms of the pressure, i.e., 
U = P/(r - 1). Therefore, the polytropic EOS can be 
rewritten as 



P 



i/r 



P 



r- 1' 



(234) 



giving a direct link between P and p. 

Without specifying the interior solution, we can al- 
ready examine the appropriate matching conditions 
needed. One needs continuity of the metric and of A' 
on the surface of the matter configuration {A is given by 
a second order differential equation). Since we know that 
the exterior solution is unique and it is the Schwarzschild- 
de Sitter solution with a cosmological constant equal to 
TZo/4, where TZq is the vacuum value of TZ (see Sec. III. Bp . 
we can directly write for the exterior 

exp(-B(r)) = iexp{A{r)) = 1 - 2™/?- - TZor^/U, 

(235) 

where £ and m are integration constants to be fixed by 
requiring continuity of the metric coefficients across the 
surface, which is implicitly defined by r = rout where 
P = p — 0. Using the definition of mtot{r) this gives, in 
the exterior. 



rutotir) = m + ^"^o- 



(236) 



On the other hand, based on the exterior solution, one 
gets on the surface 



A'(rout) 



Jlo - 12m) 



r-out CR-orl^t - 12fout + 24m) 



(237) 



Assuming that, approaching the surface from the inte- 
rior, A and mtot indeed take the correct values required 
for the matching, it can be shown that continuity of A' 
across the surface req uires -F'(rout) = for r — >• r~^^ 



(IBarausse et aLl.KOOSal). Ad ditionallv. if this is the case 



then ( Barausse all 2008a[ ) 



™tot(w) = VST^ — , (238) 
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where 



dF , ^ dF dp , , , , 

Let us now examine the behavior of m^ot ^-t the sur- 
face for different values of the polytropic index F. For 
1 < F < 3/2, C = dC/dPP' cx dC/dP {P + p) ^ a.t 
the surface so that the expression (|238p is finite and it 
gives continuity of m(ot across the surface [cf. eq. (|236|) ]. 
However, for 3/2 < F < 2, C — >■ 00 as the sur- 
face is approached, provided that dF /dTZ{TZo) ^ and 
dTZ/dT{To) ^ (note that these conditions are satisfied 
by generic forms of /(T?.), i.e., whenever an TZ^ term or 
a term inversely proportional to TZ is present). There- 
fore, even though mtot remains finite (as can be shown 
using the fact that P' = at the surface), the diver- 
gence of m^ot drives to infinity the Riemann tensor of 
the metric, R^uaX, and curvature invariants, such as R 
or R^"'^^ RfiuaXi as can easily be checked. Clearly, such 
a singular behaviour is bound to give rise to unphysical 
phenomena, such as infinite tidal forces at the surface 
{cf. the geodesic deviation equation) which would de- 
stroy anything present there. We are, therefore, forced 
to conclude that no physically relevant solution exists for 
any polytropic EOS with 3/2 < F < 2. 

The following points about the result just presented 
should be stressed: 

• The sufficient condition for the singularity to oc- 
cur is that a polytropic EOS with 3/2 < F < 2 
should adequately describe just the outer layer of 
the matter configuration (and not necessarily the 
whole configuration). 

• In practice, there is no dependence of the result 
on the functional form of f{TZ) [a few unrealisti c 
exceptions can be found in ( Barausse et a/.l . [2008a[) ] 
so what is revealed is a generic aspect of Palatini 
f{R) gravity as a class of theories. 

• The singularities discussed are not coordinate, but 
true singularities, as can be easily verified by check- 
ing that curvature invariants diverge. 



' This fact seems to have been missed in l|Barraco and Hamitvl . 
120001) . 
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• The only assumptions made regard the EOS and 
the symmetries. Thus, the result applies to all 
regimes ranging from Newtonian to strong gravity. 

Let us now interpret these results. Obviously, one 
could object to the use of the polytropic EOS. Even 
though it is extensively used for simple stellar models, 
it is not a very realistic description for stellar configu- 
rations. However, one does not necessarily have to refer 
to stars in order to check whether the issue discussed 
here leaves an observable signature. Consider two very 
well known matter configurations which are exactly de- 
scribes by a polytropic EOS: a monoatomic isentropic gas 
and a degenerate non-relativistic electron gas. For both 
of those cases F = 5/3, which is well within the range 
for which the singularities appear. Additionally, both of 
these configuration can be very well described even with 
Newtonian gravity. Yet, Palatini f{R) gravity fails to 
provide a reasonable description. Therefore, one could 
think of such matter configurations as gedanken experi- 
ments whic h reveal that Pa l atini f (R) gravity is at best 
incomplete ( Barausse et all [20Q8al lblla) . 

On the other hand, the use of the polytropic EOS re- 
quires a perfect fluid approximation for the description 
of matter. One may, therefore, wish to question whether 
the length scale on which the tidal forces become impor- 
tant is larger than the lengt h scale for which t he flui d 
approximation breaks down ( Kainulainen et "all . l2007a[) . 
However, quantitative results fo r tidal forces have been 
given m (|Barausse et am2008bD . and it has been shown 
that the length scales at which the tidal forces become 
relevant are indeed larger than it would be required for 
the fluid approximation to break down. The observable 
consequences on stellar configurations have also been dis- 
cussed there. To this, one could also add that a theory 
which requires a full description of the microscopic struc- 
ture of the system in order to provide a macroscopic de- 
scription of the dynamics is not very appealing anyway. 

In any case, it should be stated that the problem dis- 
cussed is not specific to the polytropic EOS. The use of 
the latter only simplifies the calculation and allows an 
analytic approach. The root of the problem actually lies 
with the differential structure of Palatini f{R) gravity. 

Consider the field equations in the form (|28l) : it is 
not difficult to notice that these are second order partial 
differential equations in the metric. However, since / is a 
function of TZ, and TZ is an algebraic function of T due to 
eq. (|20)) . the right hand side of eq. (|28t includes second 
derivatives of T. Now, T, being the trace of the stress 
energy tensor, will include up to first order derivatives of 
the matter fields (assuming that the matter action has 
to lead to second order field equations when varied with 
respect to the matter fields). Consequently, eq. (1^5)) can 
be up to third order in the matter fields! 

In GR and most of its alternatives, the field equations 
are only of first order in the matter fields. This guar- 
antees that gravity is a cumulative effect: the metric is 
generated by an integral over the matter sources and, 
therefore, any discontinuities (or even singularities) in 



the latter and their derivatives, which are allowed, will 
not become discontinuities o r singularities of the me tric, 
which are not allowed [see (iBarausse et all l2008bD for 
a detailed discussion]. This characteristic is not present 
in Palatini /(i?) gravity and creates an algebraic depen- 
dence of the metric on the matter fields. 

The polytropic description not only does not cause this 
problem but, as a matter of fact, it makes it less acute 
than it is, simply because in the fluid approximation the 
stress-energy tensor does not include derivatives of the 
matter fields and effectively "smoothes out" the matter 
distribution. Actually, the fact that the metric is ex- 
tremely sensitive to rapid changes of the matter field has 
been exhibited also in the i nterior of stars describ ed by 
realistic tabulated EOSs in ([Barausse et "aZl. l20Q8aD . 

One should not be puzzled by the fact that this awk- 
ward differential structure of Palatini f{R) gravity is not 
manifest in the f{Tl) formulation of the theory [and the 
field eqs. and p^ ]. We have already mentioned that 
the independent connection is actually an auxiliary field 
and the presence of auxiliary fields can always be mis- 
leading when it comes to the dynamics. In fact, it just 
takes a closer look to realize that the Palatini f{Tl) ac- 
tion does not contain any derivatives of the metric and 
is of only first order in the derivatives of the connection. 
Now, given that the connection turns out to be an auxil- 
iary field and can be algebraically related to derivatives 
of the matter and of the metric, it no longer comes as a 
surprise that the outcome is a theory with higher differ- 
ential order in the matter than the metric. 

By now, the fact that the post-Newtonian metric turns 
out to be algebraically dependent on the density, as dis- 
cussed in Sec. IV.A.3[ should no longer sound surprising: 
it is merely a manifestation of the problem discussed here 
in the weak field regime. The fact that it is unlikely that 
the Cauchy problem will be well-formulated in the pres- 
ence of matter also originates from the same feature of 
Palatini f{R) gravity, as already mentioned. Similarly, 
the fact that a theory which manifestly satisfies the met- 
ric postulates and, therefore, is expected to satisfy the 
Equivalence Principle, actually exhibits unexpected phe- 
nomenology in local non-gravitational experiments and 
conflicts with the Standard Model, as shown in Sec. I VI. 51 
also ceases to be a puzzle: the algebraic dependence of 
the connection on the derivatives of matter fields (as the 
former is an auxiliary field) makes the matter enter the 
gravitational action through the "back door". This in- 
troduces strong couplings between matter and gravity 
and self-interactions of the matter fields which manifest 
themselves in the local frame. Alternatively, if one com- 
pletely eliminates the connection (or the scalar field in 
the equivalent Brans-Dicke representation) at the level 
of the action, or attempts to write down an action which 
leads to the field eqs. (^5)) directly through metric varia- 
tion, then this action would have to include higher order 
derivatives of the matter field and self-interactions in the 
matter sector. In this sense, the fjlZ) representa tion 
is simply misleading [see also ( Sotiriou et all . l2007l) for 
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a general discussion of representation issues in gravita- 
tional theories! . 



D. Gravitational waves in f{R) gravity 

By now it is clear that the metric tensor of f{R) grav- 
ity contains, in addition to the usual massless spin 2 
graviton, a massive scalar that shows up in gravitational 
waves in the metric version of these theories (in the Pala- 
tini version, this scalar is not dynamical and does not 
propagate). A scalar gravit ational wave mode is familiar 
from scalar-tensor gravity ( Willi . Il98ll) . to which /(i?) 
gravity is equivalent. Because this scalar is massive, 
it propagates at a speed lower than the speed of light 
and of massless tensor modes and is, in principle, de- 
tectable in the arrival times of signals from an exploding 
supernova when gravitational wave detectors are suffi- 
ciently sensitive [this possibility has been pointed out as a 
discr iminator between Ten s or- Ve ctor- Scalar theories and 
GR teahva and Woodardl I2007D ]. This massive scalar 
mode is longitudinal and is of dipole nature to lowest or- 
der (ICordal l20073 IWilll . Il98ll) . The study of its genera- 
tion, propagation, and detection falls within the purview 
of scalar-tensor gravity (Will: 1981, ). The propagation of 
gravitationa l waves in the specific model f ( R) = was 
studied in ( Mendoza and Rosas-Guevaral |2Q07[ ) where 
the massive scalar mode is, however, missed. 

The generation of gravitational waves specifically in 
f{R) gravity has not received much attention in the 
literature. Even though the fact that the black hole 
solutions of GR will also be solutions of metric f{R) 
gravity (without the converse being true) implies that 
determining the geometry around a black hole is un- 
likel y to provide eviden ce for such modifications of grav- 
ity (Psaltis et all l2008l ). solutions describing perturbed 
black holes do behave differently and could, therefore, 
leave a detectable imprint on gravitational wave radia- 
tion (Barauss e and Sotiriou, 2008). Note the analogy to 
the fact that cosmological FLRW solutions are shared by 
most gravitational theories, but cosmological perturba- 
tions reveal more about the underlying theory of gravity 
than the exact solutions themselves. Additionally, grav- 
itational radiation from binary systems would probably 
be more revealing than that coming from perturbed black 
holes when it comes to modified gravity. 

The detection of gravitational waves generated in the 
theorie s .fjR) = 1/R [alrea d y rule d out by Solar Systein 
data (Barrow and_ Cliftonl . 120061: IClifton and Barrowl 
l2005a. 2006)1 and f{R}_=R±_aJf were studied in 
(|Cordal [20073^ and (jCordal [2007bl ). respectively. 

The study of cosmological gravitational waves in f{R) 
gravity is perhaps more promising than that of astro- 
physically generated waves. The stochastic gravitational 
wave bac kground produced in the earl^Miniverse was ana- 
lyzed in ()Capozziello et a^.l . 2007b. 2003). The authors of 
this last reference consider the model f{R) = R^~^^ and 
derive an evolution equation for the metric perturbations 



hij = h{t) e ' Cij in a background FLRW universe with 
scale factor a{t) — aa(^j^ : 

h+ ~ ^'^U + fc^ao ( -)\h^ 0. (240) 

This can be solved in terms of Bessel functions; plots of 
these wave amplitudes are reported in ( Capozziello et all . 
,2008) for various values of the para meter 6, but the limit 
< <5 < 7.2 • 10-19 obtained b y (iBarrow and Chftonl . 
120061 : IClifto n and BarrowL[2005all2006D leaves little hope 
of detecting /(i?) effects in the gravitational wave back- 
ground; 

( Ananda et d] . l2008l ) give a covariant and gauge- 
invariant description of gravitational waves in a per- 
turbed FLRW universe filled with a barotropic perfect 
fluid in the toy model f{R) — i?". The perturbation 
equations are solved (again, in terms of Bessel functions 
of the first and second kind) in the approximation of 
scales much larger or much smaller than the Hubble ra- 
dius H^^, finding a high sensitivity of the tensor modes 
evolution to the value of the parameter n. In particular, 
a tensor mode is found that grows during the radiation- 
dominated era, with potential implications for detectabil- 
ity in advanced space interferometers. This study, and 
others of this kind expected to appear in future literature, 
are in the spirit of discriminating between dark energy 
and dark gravity, or even between different f{R) theories 
(if this class is taken seriously) , when gravitational wave 
observations will be available: as already remarked, this 
is not possible by considering only unperturbed FLRW 
solutions. 



VII. SUMMARY AND CONCLUSIONS 
A. Summary 

While we have presented f{R) gravity as a class of toy 
theories, various authors elevate modified gravity, in one 
or the other of its incarnations corresponding to specific 
choices of the function f{R), to the role of a fully realistic 
model to be compared in detail with cosmological obser- 
vations, and to be distinguished from other models. A 
large fraction of the works in the literature is actually de- 
voted to specific models corresponding to definite choices 
of the function /(i?), and to specific parametrizations. 

Besides the power law and power series of R models 
which we have already mentioned extensively, some other 



as \n(XR) (iNoiiri and Odintsovl l2004bl: Perez Bcrehaffa 
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volved functions of R, such as f{R) = R—a(R — Ai) ™-|- 
b (R — A2)" with n,m,a,b > (jNoiiri and Odintsovl . 
I2003al) . Some models have actually been tailored to 
pass all or most of the know'n cons traints, such as the 
one proposed in ( Starobinskvl |2007[ ) where f{R) = R + 
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Ai?o[(l + R'^/Rly - 1] with n,A > and i?o being of 
the order of Hq. Here we have tried to avoid considering 
specific models and we have attempted to collect gen- 
eral, model-independent results, with the viewpoint that 
these theories are to be seen more as toy theories than 
definitive and realistic models. 

We are now ready to summarize the main results on 
f{R) gravity. On the theoretical side, we have explored 
all three versions of f{R) gravity: metric. Palatini and 
metric-affine. Several issues concerning dynamics, de- 
grees of freedom, matter couplings, etc. have been ex- 
tensively discussed. The dynamical equivalence between 
both metric/Palatini f{R) gravity and Brans-Dicke the- 
ory has been, and continues to be, a very useful tool to 
study these theories given some knowledge of the aspects 
of interest in scalar-tensor gravity. At the same time, the 
study of f{R) gravity itself has provided new insight in 
the two previously unexplored cases of Brans-Dicke the- 
ory with Wo = and wq = —3/2. We have also considered 
most of the applications of f{R) gravity to both cosmol- 
ogy and astrophysics. Finally, we have explored a large 
number of possible ways to constrain /(i?) theories and 
check their viability. In fact, many avatars of f{K) have 
been shown to be subject to potentially fatal troubles, 
such as a grossly incorrect post-Newtonian limit, short 
time scale instabilities, the absence of a matter era, con- 
flict with particle physics or astrophysics, etc. 

To avoid repetition, we will not attempt to summa- 
rize here all of the theoretical issues, the applications or 
the constraints discussed. This, besides being redundant, 
would not be very helpful to the reader, as, in most cases, 
the insight gained cannot be summarized in a sentence or 
two. Specifically, some of the constraints that have been 
derived in the literature are not model or parametrization 
independent (and the usefulness of some parametriza- 
tions is questionable). This does not allow for them to 
be expressed in a straightforward manner through simple 
mathematical equations applicable directly to a general 
fimction f{R). Particular examples of such constraints 
are those coming from cosmology (background evolution, 
perturbations, etc.). 

However, we have encountered cases in which clear-cut 
viability criteria are indeed easy to derive. We would, 
therefore, like to make a specific mention of those. A 
brief list of quick-and-easy-to-use results is: 

• In metric f{R) gravity, the Dolgov-Kawasaki in- 
stability is avoided if and only if f"{R) > 0. The 
stability condition of de Sitter space is expressed 
by eq. ([T75)) . 

• Metric f{R) gravity might pass the weak-field limit 
test and at the same time constitute an alternative 
to dark energy only if the chameleon mechanism 
is effective — this restricts the possible forms of the 
function /(i?) in a way that can not be specified by 
a simple formula. 

• Palatini f{R) gravity suffered multiple deaths, 
due to the differential structure of its field equa- 



tions. These conclusions are essentially model- 
independent. (However, this theory could poten- 
tially be fixed by adding extra terms quadratic in 
the Ricci and/or Riemann tensors, which would 
raise the order of the equations.) 

• Metric-affine gravity as an extension of the Palatini 
formalism is not sufficiently developed yet. At the 
moment of writing, it is not clear whether it suffers 
or not of the same problems that afflict the Palatini 
formalism. 

Of course, as already mentioned, the situation is often 
more involved and cannot be summarized with a quick 
recipe. We invite the reader to consult the previous sec- 
tions and especially the references that they contain. 



B. Extensions and new perspectives on f{R) gravity 

We have treated f{R) gravity here as a toy theory 
and, as stated in the Introduction, one of its merits is 
its relative simplicity. However, we have seen a number 
of viability issues related to such theories. One obvious 
way to address this issue is to generalize the action even 
further in order to avoid these problems, at the cost of 
increased complexity. Several extensions of f{R) gravity 
exist. Analyzing them in detail goes beyond the scope of 
this review, but let us make a brief mention of the most 
straightforward of them. 

We have already discussed the possibility of having 
higher order curvature invariants, such as R^^R^^ , in 
the action. In fact, from a dimensional analysis perspec- 
tive, the terms R^ and R^^R^" should appear at the 
same order. However, theories of this sort seem to be 
bur dened with wha t is called the Ostrogradski instabil- 
ity (|Woodardl . [2007h . Ostrogradski's theorem states that 
there is a linear instability in the Hamiltonians associated 
with Lagrangians which depend upon higher than first 
order derivatives in such a way that the de pendence can- 
not b e eliminated by partial integration ()Ostrogradskii 
llSSOf ). f{R) gravity seems to be the only ca se that man- 
ages to avoid this theorem ( Woodardi . I200l1) and it obvi- 
ously does not seem very appealing to extend it in a way 
that will spoil this.'^" 

The alert reader has probably noticed that the above 
holds true only for metric f{R) gravity. In Palatini f{R) 
gravity (and metric-affine f{R) gravity), as it was men- 
tioned earlier, one could add more dynamics to the action 
without having to worry about making it second order 
in the fields. Recall that, in practice, the independent 
connection is an auxiliary field. For instance, the term 
Tlfiu'Tl^'^ still contains only first derivatives of the con- 
nection. In fact, since we have traced the root of some of 



However, one could consider adding a function of the 
Gauss-Bonnet invariant Q = B? - 4_R„yfl'^'" - RaPn^R""^^" 



l lCoenoIa et all |2006| : iNoiiri and Odintsovl. boOSt) . 
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the most crucial viability issues of Palatini f{R) gravity 
to the lack of dynamics in the gravity sector, such gener- 
alizations could actually help by promoting the connec- 
tion from the rol e of an auxiliary field t o that of a truly 
dynamical field (|Barausse et all 2008b ). Such g eneral- 
izations have been considered in ( Li et 'aLl. l2007al i. 

Another extension of metric f{R) gravity that ap- 
peared recently is that in which the action includes also 
an explicit coupling between R and the matter fields . 



In ("Berto lami aLl.boOTtlBertolami and Paramosl[2007t 
iBoch nicr et all 2008[) the following action was consid- 
ered: 



S = 




A/2(i?)]ir, 



(241) 



where L„i is the matter Lagrangian and /i^2 are (a pri- 
ori arbitrary) functions of the Ricci curvature R. Since 
the matter is not minimally coupled to R, such theories 
will not lead to energy conservation and will generically 
exhibit a violation of the Equivalence Principle (which 
could potentially be controlled by the parameter A). 

The m otivation for con si dering such an action spelled 
out in (jBertolami et al. I, I2007D was that the non- 
conservation of energy could lead to extra forces, which in 
turn might give rise to phenomenology similar to Mod- 
ified Newtonian Dynamics (MOND) (jMilgroml . 1 19831 ) 
on galactic scales. Other variant s of this action have 
also been considered elsewhere: in iri and Odintsovl 
l2004al ). as an alternative to dark energy by setting 
/i (i?) = R and keeping only the nonminimal coupling of 
matter to the Ricci curvature; in (IDolgo v and Kawasaki, 
l2003bl: iMukoh^ ama and Randall l2004t) . where the idea 
of making the kinetic term of a (minimally coupled) 
scalar field dependent on the curvature, while keeping 
fi{R) — R was exploited in a ttempts to cure the cos- 
mological constant problem; in ( Bertolami and Paramosi 
120071 ) the consequences of such a theory for stellar equi- 
librium were studied; finally, generalized constraints in 
order to av oid the instabili ty discussed in Sec. IV. BTT] were 
derived in ( FaraoniL l2007al ) . The viability of theories with 
such couplings between R and matter is still under inves- 
tigation. However, the case in which bot h fi and /2 are 
linear has been shown to be non- viable (Sotirio ul. I2OO8I) 
and, for the more general case of the action (I241II. se- 
rious doubts have been expressed ( Sotiriou and Faraonil 
I2OO8D on whether extra forces are indeed present in galac- 
tic environments and, therefore, whether this theory can 
really account for the MON D-like phenomenology t hat 
initially motivated its use in ( Bertolami et all |2007^ as 
a substitute for dark matter. 

One could also consider extensions of f{R) gravity in 
which extra fields appearing in the action couple to differ- 
ent curvature invariants. A simple example with a scalar 
field is the action (11591). which is sq i netini es dubbed ex- 
tended quintessence (|Perrotta et all . l2000l ). similarly to 
the extended inflation realized with Brans-Dickc theory. 
However, such generalizations lie beyond the scope of this 



Finally, it is worth mentioning a different perspective 
on f{R) gravity. It is common in the literature that we re- 
viewed here to treat f{R) gravity as an exact theory: the 
generalized action is used to derive field equations, the so- 
lutions of which describe the exact dynamics of the grav- 
itational field (in spite of the fact that the action might 
be only an approximation an d the th eory merely a toy 
theory). A different approach (jBel and Zia, 19 85: Simon . 
1990 ) which was recently revived in (jDePeo and Psaltis . 
20071 ) is that of treating metric f{R) as an effective field 
theory. That is, to assume that the extra terms are an 
artifact of some expansion of which we are considering 
only the leading order terms. Now, when we consider a 
correction to the usual Einstein-Hilbert term, this cor- 
rection has to be suppressed by some coefficient. This 
approach assumes that this coefficient controls the order 
of the expansion and, therefore, the field equations and 
their solutions are only to be trusted to the order with 
which that coefficient appears in the action (higher or- 
ders are to be discarded) . Such an approach is based on 
two assumptions: first, some power (or function) of the 
coefficient of the correction considered should be present 
in all terms of the expansion; second, the extra degrees 
of freedom (which manifest themselves as higher order 
derivatives in metric f{R) gravity) are actually an arti- 
fact of the expansion (and there would be a cancellation 
if all orders where considered). This way, one can do 
away with these extra degrees of freedom just by proper 
power counting. Since many of the viability issues trou- 
bling higher order actions are related to the presence of 
such degrees of freedom {e.g., classical instabilities), re- 
moving these degre es of freedom could certa inly allevi- 
ate many problems ( DeDeo and Psalti^l2007t ). However, 
the assumptions on which this approach is based should 
not be underestimated either. For instance, early results 
that showed renormalization of higher order actions were 
based on an exact treatmen t, i.e., it is fo urth order grav- 
ity that is renormalizable ( Stelld . [Ht^). Even though, 
from one hand, the effective field theory approach seems 
very reasonable (these actions are regarded as low energy 
limits of a more fundamental theory anyway), there is no 
guarantee that extra degrees of freedom should indeed 
not be present in a non-perturbative regime. 



C. Concluding remarks 

Our goal was to present a comprehensive but still thor- 
ough review of f{R) gravity in order to provide a starting 
point for the reader less experienced in this field and a 
reference guide for the expert. However, even though we 
have attempted to cover all angles, no review can replace 
an actual study of the literature itself. It seems inevitable 
that certain aspects of f{R) might have been omitted, or 
analyzed less than rigorously and, therefore, the reader 
is urged to resort to the original sources. 

Although many shortcomings of f{R) gravity have 
been presented which may reduce the initial enthusiasm 
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with which one might have approached this field, the fact 
that such theories arc mostly considered as toy theories 
should not be missed. The fast progress in this field, 
especially in the last five years, is probably obvious by 
now. And very useful lessons, which have helped signif- 
icantly in the understanding of (classical) gravity, have 
been learned in the study of f{R) gravity. In this sense, 
the statement made in the Introduction that f{R) grav- 
ity is a very useful toy theory seems to be fully justified. 
Remarkably, there are still unexplored aspects of f{R) 
theories or their extensions, such as those mentioned in 
the previous section, which can turn out to be fruitful. 
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